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
>
> ______________________________________________
> [email protected] 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 |
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
