[R] problem in loop

2009-09-02 Thread Seunghee Baek
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

[R] double bootstrap

2009-07-01 Thread Seunghee Baek
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 bootstra

[R] 95% confidence interval coverage

2009-06-26 Thread Seunghee Baek
Hi R-users, I would like to compute 95%CI coverage for each bootstrapped sample. But, I don't know how to get 95% CI of estimate for each bootstrapped sample. I could only see the estimate of each bootstrapped sample. Thanks, Becky __ R-help@r-project.

[R] a difficulty in boot package

2009-06-19 Thread Seunghee Baek
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

[R] resampling the entire row

2009-06-19 Thread Seunghee Baek
Hi, For bootstrapping method, I would like to resample the entire row instead of one column. What should I do? Thanks, Becky __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R