我正在测试 Iris 数据集中的多元正态性。我想比较 MVN 包下所有不同测试的结果(Mardia 的测试、Henze-Zikler、Royston、Doornik-Hansen 和 Energy 测试)并完成了以下代码:
library(MVN)
library(htmlTable)
attach(iris)
#Mardia's Test
Mar<-MVN::mvn(data=iris,subset="Species", mvnTest="mardia")
Mar$multivariateNormality
#Henze-Zirkler's Test
Hz<-MVN::mvn(data=iris,subset="Species",mvnTest="hz")
Hz$multivariateNormality
#Royston's Test
Roy<-MVN::mvn(data=iris,subset="Species", mvnTest="royston")
Roy$multivariateNormality
#Doornik-Hansen Test
Dh<-MVN::mvn(data=iris,subset="Species", mvnTest="dh")
Dh$multivariateNormality
我必须以一种易于在我的发现中呈现的方式获得结果。我在三个表格(每个物种一个)之后总结了上面输出的数据。基本上对于每个表,行是不同的测试,列是测试统计值、P 值和 mvn 结果。我尝试了以下
table <- rbind(Mar$multivariateNormality$setosa[1:2, ],Hz$multivariateNormality$setosa[1:2, ],Roy$multivariateNormality$setosa[1:2, ])
htmlTable(table, caption="Testing")
但是,我收到此错误
Error in match.names(clabs, names(xi)) :
names do not match previous names