[R] boot package
Angelo Canty
canty at math.mcmaster.ca
Thu Jan 6 16:34:41 CET 2005
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:n.l.pace at utah.edu
> Department of Anesthesiology Home:nlpaces at comcast.net
> 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 at 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: cantya at mcmaster.ca |
| Mathematics and Statistics Phone: (905) 525-9140 x 27079 |
| McMaster University Fax : (905) 522-0935 |
| 1280 Main St. W. |
| Hamilton ON L8S 4K1 |
More information about the R-help
mailing list