小编MYa*_*208的帖子

tidyverse:绑定相同维度的列表元素

使用reduce(bind_cols),可以组合相同维度的列表元素。但是,我想知道如何从可能具有不同维度元素的列表中仅组合相同维度(可能以某种方式指定维度)的元素。

library(tidyverse)

df1 <- data.frame(A1 = 1:10, A2 = 10:1)
df2 <- data.frame(B = 11:30)
df3 <- data.frame(C = 31:40)

ls1 <- list(df1, df3)
ls1

[[1]]
   A1 A2
1   1 10
2   2  9
3   3  8
4   4  7
5   5  6
6   6  5
7   7  4
8   8  3
9   9  2
10 10  1

[[2]]
    C
1  31
2  32
3  33
4  34
5  35
6  36
7  37
8  38
9  39 …
Run Code Online (Sandbox Code Playgroud)

reduce r dplyr purrr tidyverse

16
推荐指数
4
解决办法
564
查看次数

ggplot2:标签表达式中的字体样式

更改字体面和大小非常有效,而无需在标签中进行表达.这是代码,我有问题更改标签表达式中的字体面.

p <- ggplot(data = mtcars, aes(x=wt, y=mpg)) + 
    geom_point() + 
    labs(x="Weight", y=expression(paste("mpg (  ", m^{-2}, ")"))) + 
    theme(axis.title.x = element_text(family="serif", face="bold", 
        size=12, angle=00, hjust=0.54, vjust=0)) + 
    theme(axis.title.y = element_text(family="serif", face="bold", 
        size=12, angle=90, vjust=0.25))

p
Run Code Online (Sandbox Code Playgroud)

产量

任何帮助将受到高度赞赏.谢谢

r ggplot2

15
推荐指数
1
解决办法
2万
查看次数

在ggplot2中将facet标签更改为数学公式

我想知道如何将facet标签更改为数学公式ggplot2.

d <- ggplot(diamonds, aes(carat, price, fill = ..density..)) +
  xlim(0, 2) + stat_binhex(na.rm = TRUE) + opts(aspect.ratio = 1)
d + facet_wrap(~ color, ncol = 4)
Run Code Online (Sandbox Code Playgroud)

在此输入图像描述

例如,我想将facet标签从Dto 更改为Y[1],其中1是下标.在此先感谢您的帮助.

我找到了这个答案,但它对我不起作用.我正在使用R 2.15.1ggplot2 0.9.1.

r ggplot2

15
推荐指数
4
解决办法
5530
查看次数

从R中的lsmeans包中提取lsmobj的元素

我想知道如何lsmobjlsmeans包中获取元素R.

require(lsmeans)
fiber.lm <- lm(strength ~ diameter + machine, data = fiber)
fiber.lsm <- lsmeans (fiber.lm, "machine")
fiber.lsm
machine   lsmean        SE df lower.CL upper.CL
 A       40.38241 0.7236252 11 38.78972 41.97510
 B       41.41922 0.7444169 11 39.78077 43.05767
 C       38.79836 0.7878785 11 37.06426 40.53247

Confidence level used: 0.95 


str(fiber.lsm)
'lsmobj' object with variables:
    machine = A, B, C
Run Code Online (Sandbox Code Playgroud)

我想要fiber.lsm的提取物lsmeansSE列.任何帮助将受到高度赞赏.在此先感谢您的帮助.

r extraction lsmeans

15
推荐指数
1
解决办法
5213
查看次数

lme4中的多元线性混合模型

我想知道如何拟合多元线性混合模型lme4.我使用以下代码拟合单变量线性混合模型:

library(lme4)
lmer.m1 <- lmer(Y1~A*B+(1|Block)+(1|Block:A), data=Data)
summary(lmer.m1)
anova(lmer.m1)

lmer.m2 <- lmer(Y2~A*B+(1|Block)+(1|Block:A), data=Data)
summary(lmer.m2)
anova(lmer.m2)
Run Code Online (Sandbox Code Playgroud)

我想知道如何拟合多元线性混合模型lme4.数据如下:

Block A B    Y1    Y2
 1    1 1 135.8 121.6
 1    1 2 149.4 142.5
 1    1 3 155.4 145.0
 1    2 1 105.9 106.6
 1    2 2 112.9 119.2
 1    2 3 121.6 126.7
 2    1 1 121.9 133.5
 2    1 2 136.5 146.1
 2    1 3 145.8 154.0
 2    2 1 102.1 116.0
 2    2 2 112.0 121.3
 2 …
Run Code Online (Sandbox Code Playgroud)

r lme4

14
推荐指数
2
解决办法
1万
查看次数

将小平面标题注释为小平面上的条带

我想在一个刻面的情节中添加一个方面标题作为条带ggplot2.我的MWE抛出一个错误.任何帮助将受到高度赞赏.谢谢

library(ggplot2)
library(gtable)
p <- ggplot(mtcars, aes(mpg, wt)) + geom_point()
p <- p + facet_grid(. ~ cyl)


# get gtable object
Plot1 <- ggplot_gtable(ggplot_build(p))
# add label for top strip
Plot1 <- gtable_add_rows(Plot1, Plot1$heights[[3]], 2)
Plot1 <- gtable_add_grob(Plot1,
  list(rectGrob(gp = gpar(col = NA, fill = gray(0.5))),
  textGrob("Cyl", gp = gpar(col = gray(1)))),
  3, 4, 3, 10, name = paste(runif(2)))

# add margins
Plot1 <- gtable_add_rows(Plot1, unit(1/8, "line"), 2)

# draw it
grid.newpage()
print(grid.draw(Plot1))
Run Code Online (Sandbox Code Playgroud)

r facet ggplot2 gtable

14
推荐指数
1
解决办法
2348
查看次数

用ggpmisc显示nls模型的方程

Rggpmisc可用于显示lm模型和poly模型的方程ggplot2(参见此处参考).想知道如何使用nls模型方程结果.以下是我的MWE.ggplot2ggmisc

library(ggpmisc)
args <- list(formula = y ~ k * e ^ x,
             start = list(k = 1, e = 2))
ggplot(mtcars, aes(wt, mpg)) +
  geom_point() +
  stat_fit_augment(method = "nls",
                   method.args = args)
Run Code Online (Sandbox Code Playgroud)

r nls ggplot2 lm ggpmisc

14
推荐指数
2
解决办法
469
查看次数

ggplot2中来自glm和stat_smooth的逻辑回归的预测值是不同的

我正在努力制作这个逻辑回归图ggplot2.

df <- structure(list(y = c(2L, 7L, 776L, 19L, 12L, 26L, 7L, 12L, 8L,
24L, 20L, 16L, 12L, 10L, 23L, 20L, 16L, 12L, 18L, 22L, 23L, 22L,
13L, 7L, 20L, 12L, 13L, 11L, 11L, 14L, 10L, 8L, 10L, 11L, 5L,
5L, 1L, 2L, 1L, 1L, 0L, 0L, 0L), n = c(3L, 7L, 789L, 20L, 14L,
27L, 7L, 13L, 9L, 29L, 22L, 17L, 14L, 11L, 30L, 21L, 19L, 14L,
22L, 29L, 28L, 28L, 19L, 10L, 27L, 22L, 18L, 18L, …
Run Code Online (Sandbox Code Playgroud)

r ggplot2

13
推荐指数
2
解决办法
1万
查看次数

使用knitr为R代码添加行号

我不知道是否有把任何功能line numbersknitr.Rnw.我发现了这个讨论和这些文件(这里这里),但找不到放置行号的方法.任何帮助将受到高度赞赏.谢谢

latex r sweave knitr rnw

13
推荐指数
1
解决办法
2510
查看次数

带有月份和日期的时间序列日期在ggplot2中与facet一起绘制多年

我想有两个month,并day在使用的时间序列图的x轴facet的年ggplot2.我的MWE如下:

set.seed(12345)
Date <- seq(as.Date("2010/1/1"), as.Date("2014/1/1"), "week")
Y <- rnorm(n=length(Date), mean=100, sd=1)
df <- data.frame(Date, Y)

df$Year <- format(df$Date, "%Y")
df$Month <- format(df$Date, "%b")
df$Day <- format(df$Date, "%d")

df$MonthDay <- format(df$Date, "%d-%b")


p <- ggplot(data=df, mapping=aes(x=MonthDay, y=Y, shape=Year, color=Year)) + geom_point() +geom_line(aes(group = 1))
p <- p + facet_grid(facets = Year ~ ., margins = FALSE) + theme_bw()
print(p)
Run Code Online (Sandbox Code Playgroud)

在此输入图像描述

我尝试使用以下命令控制x轴标签

p + scale_y_continuous() + scale_x_date(labels = date_format("%d-%b"))
Run Code Online (Sandbox Code Playgroud)

但它会抛出以下错误消息.

Error: Invalid input: date_trans …
Run Code Online (Sandbox Code Playgroud)

r time-series ggplot2

12
推荐指数
2
解决办法
1万
查看次数

标签 统计

r ×10

ggplot2 ×6

dplyr ×1

extraction ×1

facet ×1

ggpmisc ×1

gtable ×1

knitr ×1

latex ×1

lm ×1

lme4 ×1

lsmeans ×1

nls ×1

purrr ×1

reduce ×1

rnw ×1

sweave ×1

tidyverse ×1

time-series ×1