[R] bca ci's and NaN's in boot.out
Thomas W Blackwell
tblackw at umich.edu
Fri Apr 25 00:04:48 CEST 2003
I think that calling bca.ci() with argument type="reg"
in addition to the arguments you've used already MIGHT do
This argument will be passed to empinf() and it ought
to guarantee that the method empinf.reg() is used. That
method looks as though it will tolerate NAs. I've never
used any of this stuff, so I don't know. I'm just looking
at the function definitions for bca.ci(), empinf() and
empinf.reg(), inf.jack(), usual.jack(), positive.jack().
These last four functions are called by empinf() under
For more details, I encourage you to look at the function
definitions yourself, and try to trace exactly which ones
will be called for the particular set of arguments you are
using. To see a function definition, type the name of the
function at the command prompt without any parentheses ()
after the function name.
It's quite possible that you are testing the package with
a combination of arguments that has never been used before,
and that some default setting isn't quite doing its job.
I'm afraid I will leave the more detailed investigation to
HTH - tom blackwell - u michigan medical school - ann arbor -
On Thu, 24 Apr 2003, Katalin Csillery wrote:
> Dear All,
> I am trying to use the bca.ci function on a boot.out object
> which consists a few NaN's and I want to ignore those NaN's,
> and get a ci only for the "normal" values.
> boot.out$t has R number of values for 3000 different statistics,
> so when I use boot.ci(boot.out, index=i) and i happens to be a
> column in boot.out$t with some NaN's in there I get an error
> message. I tried na.omit for the problematic columns, but bca.ci
> needs a boot.out object so I cannot feed in a single column.
> The number of NaN's varies, and for some statistics
> there isn't any.
> Interestingly the standard errors that the boot() calculates
> and returns there is an automatic NA exclusion.
> (A side question, in S+ the se's are saved in boot.out$statistics$SE, but
> in R boot.out I couldn't find them. Are they saved somewhere in boot.out?)
> Any idea would be greatly appreaciated!
> Katalin Csillery
> Division of Biological Sciences
> University of Montana, Missoula MT 59801
> Phone: 406 243 6106, E-mail: csillery at selway.umt.edu
More information about the R-help