resamp.std {agce} | R Documentation |
Bootstrap a new data set from the original data set by keeping the general structure of the groups. Then compute the standard deviation based on the bootstrap estimates.
resamp.std(factors, Y, X, B=1000)
factors |
A vector (character or numeric) containing the group ids. |
Y |
The matrix of measurements. |
X |
The design matrix. |
B |
A numeric value, the number of bootstraps. |
The matrix of standard deviation.
Raphael Gottardo
data(Bt20) ###Create the design matrix X<-design.matrix(Bt20[,1]) ###Convert as numeric the measurement Y<-as.matrix(Bt20[,3:dim(Bt20)[2]]) ### Number of groups p<-dim(X)[2] ### Compute the standard deviation associated with the mean std<-resamp.std(Bt20[,1],Y,X,B=1000)