Displaying 1 result from an estimated 1 matches for "data_hn".
Did you mean:
data_ac
2006 Jan 12
0
question for mshapiro test
...R code I used and the result.
library(mvnormtest,lib.loc="~/mshapiro")
library(MASS)
n=50
dim=2
ntrial = 1000
x<-matrix(1:(dim*n),ncol=dim)
count = 0
for( trial in 1:ntrial)
{
x<-mvrnorm(n,rep(0,dim),diag(1,dim,dim))
data_hn<-x
p<-mshapiro.test(t(data_hn))$p.value
if( p<= 0.05) count <- count+1
}
print(count)
[1] 117
Can you help me out? Thank you very much.
Liang