Displaying 1 result from an estimated 1 matches for "coefficientsboot".
2008 May 04
1
Residual resampling for non linear reg model
...ad.table("http://www.creem.st-and.ac.uk/len/classes/mt3607/data/manhours_surgical.dat", header=TRUE)#attach dataattach(manpower)B<-999#number of data pointsn<-dim(manpower)[1]#alpha level to use for the confidence limitsalpha<-0.05#matrix that's going to contain the bootstrap coefficientsboot.coef<-matrix(NA, B+1, 2)#fit the initial modelinit.values<-c(alpha=20,beta=0)model<-nls(manhours~cases/(alpha+(beta*cases)), start=init.values, trace=TRUE)pred<-predict(model)resid<-resid(model)#do the bootstrapfor (i in 1:B){ #resample the residuals resample.index<-sample(1:n,n...