我有这个数据:
a <- c(12, 14, 16, 12, 14, 16, 12, 14, 16)
b <- c(107.13, 83.01, 164.01, 113.01, 85.48, 86.08, 112.96, 82.18, 51.92)
a
是具有3个级别12、14和16的因子,b
是我的响应变量
当我尝试运行和绘制pairs.anova
时
pairw.anova(y=a, x=b, conf.level = 0.95, method = "tukey", MSE = NULL, df.err = NULL)
我得到这个错误:
dims[产品0]与对象[1]的长度不匹配
你知道为什么吗?
谢谢你的帮助。
也许您可以交换分段的值,因为响应变量(定量(可能是b
,而包含组的分类向量可能是a
。
你必须确保a
是一个因素:
library(asbio)
pairw.anova(y=b, x=factor(a), conf.level = 0.95,
method = "tukey", MSE = NULL, df.err = NULL)
95% Tukey-Kramer confidence intervals
Diff Lower Upper Decision Adj. p-value
mu12-mu14 27.47667 -55.80207 110.7554 FTR H0 0.596586
mu12-mu16 10.36333 -72.9154 93.64207 FTR H0 0.923774
mu14-mu16 -17.11333 -100.39207 66.1654 FTR H0 0.809506