Hi R-users, I've found the error.... In these rows: # Model (a) testtemp <- lm.bp(doctorco~sex+age+income, prescrib~sex+age+income, data=bootdata) betafound <- c(testtemp$beta,testtemp$beta3) results[i,] <- betafound
betafound must to be equal to: betafound <- c(testtemp$beta1,testtemp$beta2,testtemp$beta3) ... Thanks a lot to everyone... -- View this message in context: http://r.789695.n4.nabble.com/Bootstraps-standard-error-tp3313322p3317675.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
