0

我为此工作了一天，希望你们能对这个奇怪的事情给出提示。

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
4

1 回答 1

0

我真的不知道为什么，但我开始运行你的代码，前三行似乎揭示了为什么你会得到不同的结果。如果你运行那部分代码然后简单地问

y0
y1
y2

您将看到所有三个对象都有不同的元素。事实上，所有元素都是不同的。至于为什么会这样，我不确定您是否以相同的方式定义它们，但它们肯定是不同的。把它们画出来看看。

希望这可以帮助

于 2015-10-23T00:28:59.250 回答