我想在r中创建一个boxplot,但我很难将boxplot分组为第三个变量。
我希望在一个y轴(治疗结果调查问卷上得分)(称为HAMD)上有箱形图 - )
周0-5(称为周)
两个治疗组(称为治疗组)均为两次。
我用
制作了两个箱图boxplot(MP $ HAMD~MP $周)
和六个
boxplot(MP $ HAMD~MP $ Treatment)
但现在我想要12个箱图,每个治疗每周一次。 我怎么能在r?
中这样做非常感谢你。
答案 0 :(得分:3)
您可以使用interaction
功能和基本图形执行此操作:
boxplot( HAMD ~ interaction(treatment,week), data=MP )
boxplot( HAMD ~ interaction(week,treatment), data=MP )
boxplot( HAMD ~ interaction(week,treatment), data=MP,
at= c(1:6, 8:13) )
以下是使用网格包的一个选项:
library(lattice)
bwplot( HAMD ~ week|treatment, data=MP )
bwplot( HAMD ~ treatment|week, data=MP )
bwplot( HAMD ~ treatment|week, data=MP, layout=c(6,1) )
使用ggplot2包的选项:
library(ggplot2)
p <- qplot(interaction(treatment,week), HAMD, data=MP, geom="boxplot")
p
p + aes(fill=week)