search for: yboot

Dear R users, I want to conduct a small simulation study and I have to use the lme function in a loop to save the restricted log likelihood. However, for one simulated data set the lme function gives this error Error en lme.formula(yboot ~ X[, -1], data = data.fr, random = Z.block) : nlminb problem, convergence error code = 1 message = singular convergence (7) and then, the simulation stops. I would like to skip this simulated data, and then continue with the loop, but I cannot find the way of doing this. I am using R2...

...ment to be the same (or eliminate this line altogether), I get reproducible results.=20=20 =20 I'm using R 2.2.0 =20 Thanks, -Lori =20 test.fun<-function(nsamp, mu,q, seed, nsim=3D2){ set.seed(seed)=20=20=20=20 mnsq<-rep(NA,nsim) for (i in 1:nsim){ y<-rnorm(nsamp,mean=3Dmu) yboot<-sample(y, size=3Dq, replace=3DT) mnsq[i]<-(mean(y)^2) } return(mnsq) } =20 test.fun(50,0,30,12345) test.fun(50,0,10,12345) =20 Lori E. Dodd, Ph.D. National Cancer Institute Biometric Research Branch 6130 Executive Blvd, MSC 7434 Rockville, MD 20892=20 (For overnight deliveries...

...lues that I can then use to do parametric bootstrap analysis and generate the confidence intervals? Something like this: n<-length(A) niter<-1000 y<-matrix(nrow=n,ncol=niter*2) for (i in 1:niter) { y[,I(i*2-1):I(i*2)]<-simulate(model)[,1] } rvalues<-numeric() for (i in 1:niter) { yboot<-cbind(y[,I(i*2-1)],y[,I(i*2)]) mboot<-lmer(y~A+B+(1|C/D)+(1|E),binomial) rvalues[i]<- attr(lme4::VarCorr(mboot)$E, "stddev")^2/(1*(pi^2)/3 + attr(lme4::VarCorr(mboot)$E, "stddev")^2 + attr(lme4::VarCorr(mboot)$C, "stddev")^2 + attr(lme4::VarCorr(mboot)$D, &...