Hello,
I want to draw a histogram of the mean of sample observations drawn from
multivariate t distribution. I am getting the following error corresponding to
the code I used. Though I am getting the graph, but I am curious to know the
warning message.
Warning messages:
1: In if (freq) x$counts else { :
the condition has length > 1 and only the first element will be used
2: In if (!freq) "Density" else "Frequency" :
the condition has length > 1 and only the first element will be used
-----------------------------------------------------------------------------------
library(mvtnorm);
s12<-c(1:1000);
var21<-lapply(s12,function(x){
rs<-rmvt(10, sigma=diag(8), df=1);
my<-mean(rs);
sy<-sqrt(var(rs))
return(cbind(my,sy))
});
data1<-do.call(rbind,var21);
dataMat<-data.frame(data1);
W<-dataMat$my;
hist(W,breaks=20,probability=T)
-----------------------------------------------------------------------------------
Thanks.
Shant
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
