#dset<-read.table('MardiaMV25.dat.txt', na.string='-99')
#dset<-data.matrix(dset)
#n<-dim(dset)[1]
#p<-dim(dset)[2]
#miss_pattern<-rsem.pattern(n,p,dset)
#misinfo<-miss_pattern$misinfo
#V_forana<-c(1,2,4,5)
#em_results<-rsem.emmusig(dset,misinfo)
#hmu1<-em_results$mu
#hsigma1<-em_results$sigma
#rsem.Ascov(x, hmu1, hsigma1)
Run the code above in your browser using DataLab