sim.list2<-fun.sim(g=c(20,50,10),mean1=0.5,d=0,sds2=c(1,1,1,1),corr=c(0.1,0.5,0.5,0),
n.var=c(20,20))
datos2 <- as.data.frame(sim.list2$x)
datos2<-subset(datos2,select=-grp)
grupos <- sim.list2$grp
grupos<-factor(grupos,labels=c(1,2,3))
sde.method(mydata=datos2,groups=grupos,plt=FALSE)
Run the code above in your browser using DataLab