通过这个图 可以判定数据符合多元正态分布吗

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

如图 ，从一个5列的数据里，通过自己写的qqmultinorm函数得到的图，图中是经验分布函数的100个p值与[0,1]均匀分布的经验分布函数的比较。从这张图可以看出此数据集是多变量正态吗？

我已经用shapiro检验了，其中4个符合正态分布，1个不符合。。


贴一下这个函数的代码：
> qqmultinorm <- function(x,nproj=50,scale=T,plot.qq=F,plot.edf=F) {
+     p <- ncol(x)
+     if(scale) x <- scale(x,center=T,scale=T)
+     if(plot.qq) devAskNewPage(ask = T)
+     pvals <- NULL
+     for (i in 1:nproj) {
+         a <- rnorm(p)
+         a <- a/sqrt(sum(a^2))
+         y <- as.vector(x%*%a)
+         if (plot.qq) qqnorm(y)
+         r <- shapiro.test(y)
+         pvals <- c(pvals,r$p.value)
+     }
+     if (plot.edf) {
+         plot(ecdf(pvals),xlab="p-values",ylab="probability",
+              main=" ")
+         abline(0,1,lwd=2)
+     }
+     pvals
+ }
> r = qqmultinorm(cbind(mec,vec,alg,ana,sta),nproj=100,plot.edf = T)

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：多元正态分布 正态分布 shapiro multi 分布函数 正态分布

try test using: "nortest" package in R.  Compare with shapiro.test() result.