Hi Nathan, sorry about the delay in responding.
The test of whether all of the bootstrap values are equal is in the wrong place. I will correct it in the next release. Meanwhile, I think the simplest solution is to fix(boot.ci) and move the line if (const(t, min(1e-08, mean(t)/1e+06))) { print(paste("All values of t are equal to ", mean(t), "\n Cannot calculate confidence intervals")) return(NULL) } from its current position to directly after the line t <- t[fins] boot.ci should then work properly for missing statistic values. Angelo On Sat, 4 Dec 2004, Nathan Leon Pace, MD, MStat wrote: > Hi, > > I using the boot package 1.2-20 on R 2.0.1. > > My statistics function estimates 6 parameters. > > In a small percentage of resampled data sets my statistics function > doesn't produce an estimate for one parameter and the boot function > stops with an error. > > I can write an ifelse(exists('parameter.estimate'), parameter.estimate, > NA) statement within the statistic function to substitute an NA for the > missing estimate value. > > However, the boot.ci function to generate CIs from the boot object > won't accept NAs. > > My problem is writing code to impute a numeric value for the missing > estimate. ifelse won't generate a numeric value if the test is mode > logical. > > Any suggestions? > > Nathan > > > Nathan Leon Pace, MD, MStat Work:[EMAIL PROTECTED] > Department of Anesthesiology Home:[EMAIL PROTECTED] > University of Utah Work:801.581.6393 > Salt Lake City, Utah Home:801.467.2925 > Fax:801.581.4367 > Cell:801.558.3987 > > ______________________________________________ > R-help@stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html > -- ------------------------------------------------------------------ | Angelo J. Canty Email: [EMAIL PROTECTED] | | Mathematics and Statistics Phone: (905) 525-9140 x 27079 | | McMaster University Fax : (905) 522-0935 | | 1280 Main St. W. | | Hamilton ON L8S 4K1 | ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html