我在这方面工作了一天，并希望你们能对这个奇怪的事情做出暗示 .
y0，y1和y2通过相同的方法独立生成 .
它们通过相同的方法分成20组 .
然而，manova说它们有很大的不同？为什么？
Manova测试（存储在变量s中）的摘要说：
Pr（> F）值小于2.2e-16 .
y0 <- runif(100, 0, 1)
y1 <- runif(100, 0, 1)
y2 <- runif(100, 0, 1)
y0 <- c(y0, runif(100, 0, 10) )
y1 <- c(y1, runif(100, 0, 10) )
y2 <- c(y2, runif(100, 0, 10) )
y0=as.numeric(unlist(y0))
y1=as.numeric(unlist(y1))
y2=as.numeric(unlist(y2))
b=10
a=length(y0)/b
g=rep(1:a,rep(b,a))
m1 <- manova(cbind(y0, y1, y2) ~ g)
s=summary(m1, test = "Wilks")
a = s$stats
a = a[11]
s
a
摘要如下：
Df Wilks approx F num Df den Df Pr(>F)
g 1 0.37069 110.91 3 196 < 2.2e-16 ***
Residuals 198
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
1 回答
我真的不确定为什么，但我开始运行你的代码，前三行似乎揭示了为什么你会得到不同的结果 . 如果您运行该部分代码然后只是询问
你会看到所有三个对象都有不同的元素 . 事实上，所有元素都是不同的 . 至于为什么会这样，我不确定是否会以同样的方式定义它们，但它们肯定是不同的 . 将它们绘制出来并看一看 .
希望这可以帮助