11

我的旧代码如下所示：

library(ggplot2)
gp<-ggplot(NULL,aes(x=Income))
gp<-gp+geom_density(data=dat$Male,color="blue")
gp<-gp+geom_density(data=dat$Female,color="green")
gp<-gp+geom_density(data=dat$Alien,color="red")
plot(gp) #Works

现在我已经开始使用优秀的 data.table 库（而不是 data.frame）：

library(data.table)
cols<-c("blue","green","red")
gp<-ggplot(NULL,aes(x=Income))
dat[, list(gp+geom_density(data=.SD, color=cols[.GRP])), by=Gender]
#I even tried
dat[, list(gp<-gp+geom_density(data=.SD, color=cols[.GRP])), by=Gender]
plot(gp) #Error: No layers in plot

我不确定出了什么问题，但似乎我在J()中运行的代码在外部范围内没有被识别。

如何以 data.table 惯用的方式实现这一点？

4

1 回答 1

12

ggplot2应该以与长格式 data.frames 相同的方式与长格式 data.tables 一起使用：

library(data.table)
set.seed(42)

dat <- rbind(data.table(gender="male",value=rnorm(1e4)),
             data.table(gender="female",value=rnorm(1e4,2,1))
             )

library(ggplot2)
p1 <- ggplot(dat,aes(x=value,color=gender)) + geom_density()
print(p1)

不要将宽格式 data.frames（或 data.tables）提供给 ggplot2。

如果你有很多组，绘图会很慢，但由于它的内在魔力ggplot2，没有什么data.table能真正帮助你（直到 Hadley 以某种方式实现它）。您可以尝试计算外部的密度ggplot2，但这只会对您有所帮助：

set.seed(42)
dat2 <- data.table(gender=as.factor(1:5000),value=rnorm(1e7))
plotdat <- dat2[,list(x_den=density(value)$x,y_den=density(value)$y),by=gender]
p2 <- ggplot(plotdat,aes(x=x_den,y=y_den,color=gender)) + geom_line()
print(p2) #this needs some CPU time

当然，如果你有很多组，你可能会做错类型的情节。

于 2013-03-20T16:10:37.023 回答