我是R的新手,所以请原谅我的业余编码语言。我正在执行三向方差分析。 3个因子,每个因子有2个等级。我正在尝试使用ggplot绘制图形。我有执行2方差分析(2种水平2种治疗)的代码,想知道如何添加另一个因素。
我使用的代码是:
data.summary2 <- ddply(dat2,
.(SEDIMENT, NUTRIENTS ),
summarise,
mean = mean(Total_Invert_abundance),
se = se(Total_Invert_abundance))
此代码创建一个图表,显示四种治疗组合中每种组合的均值+-标准错误
ggplot(data.summary2, aes(x = SEDIMENT, y= mean, fill = NUTRIENTS)) +
geom_bar(position = "dodge", stat = "identity")+
geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
position=position_dodge(0.9), width = 0.2)
您可以看到泥沙和营养素这两个因素。我需要添加另一个名为Flow_Velocity的变量,该变量已经在其他箱形图中使用,并且在我的数据集中。因此,只需要知道如何将其添加到ggplot geom栏中即可。不确定是否有可能?
在x = area
中,我添加了不起作用的flow_velocity,然后在fill =区域中,我也添加了额外因子的名称,但它也不起作用
data.summary2 <- ddply(dat2,
.(SEDIMENT, NUTRIENTS ),
summarise,
mean = mean(Total_Invert_abundance),
se = se(Total_Invert_abundance))
此代码创建一个图形,显示每个图形的均值+-标准错误 四种治疗组合中的
ggplot(data.summary2, aes(x = SEDIMENT, y= mean, fill = NUTRIENTS)) +
geom_bar(position = "dodge", stat = "identity")+
geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
position = position_dodge(0.9), width = 0.2)
我希望能够制作相同的图,但是在运行三向方差分析之前添加其他因素作为探索性分析