引用于:公众号“R语言和统计”,原标题为“这篇《新英格兰医学杂志》论文中的美丽的图到底是怎么画的?附R代码‘
> mydata<-data.frame(
+ ID=c(rep(1:10,2),rep(11:20,2)),
+ treatment=rep(c("v","p"),each=20),
+ time=rep(c("before","after","before1","after1"),each=10),
+values=c(rnorm(10,mean=100,sd=30),rnorm(10,mean=150,sd=30),rnorm(10,mean=100,sd=30),rnorm(10,mean=100,sd=30))
+ )
head(mydata,20)
tail(mydata,20)
ggplot(mydata,aes(time,values,fill=time))+
+ geom_point(shape=21