search for: seunghee

Hi R-users, I have a problem for updating the estimates of correlation coefficient in simulation loop. I want to get the matrix of correlation coefficients (matrix, name: est) from geese by using loop(500 times) . I used following code to update, nsim<-500 est<-matrix(ncol=2, nrow=nsim) for(i in 1:nsim){ fit <- geese(x ~ trt, id=subject, data=data_gee, family=binomial,

Hi, I have a problem in programming for bootstrapping. I don't know why it show the error message. Please see my code below: #'st' is my original dataset. #functions of 'fml.mlogl','pcopula.fam4','ltd','invltd' are already defined boot.OR<-function(data,i) { E=data[i,] ml1<-glm(c_VAsex90_bf ~ trt,family=binomial,data=E) ml2<-glm(c_VAsex90_bm ~

Hi All, I would like to do double boostrapping to estimate 95% CI coverage. So, I can only need to estimate 95% confidence interval from each bootstrapped sample. Since we don't have a closed form of 95% CI, in order to get 95% CI for each sample, we need to use bootstrapping. For outer bootstrapping, I used boot function in library(boot). How can I add double bootstrap inside the function.