search for: yearsmarri

...Affairs) # function to obtain regression weights bs <- function(formula, data, indices) { d <- data[indices,] fit <- tobit(formula, data=Affairs) return(coef(fit)) } # bootstrapping with 1000 replications results <- boot(data=Affairs, statistic=bs, R=1000, formula=affairs~age+yearsmarried+religiousness+occupation+rating) # view results results That gives me the right coefficients, but the SEs are all zero, so I cannot get 95% confidence intervals... Any help is welcome... Ravi -- View this message in context: http://r.789695.n4.nabble.com/Bootstrapped-Tobit-regression-g...