[R] bcanon bug???

[Bruno Giordano]

Hallo,
I'm using bcanon to get bootstrap confidence intervals.
I need to modify the default confidence levels
(alpha=c(0.025, 0.05, 0.1, 0.16, 0.84, 0.9, 0.95, 0.975)),

adding two more extreme values:

alpha=c((0.05/x),0.025, 0.05, 0.1, 0.16, 0.84, 0.9, 0.95, 0.975,1-(0.05/x)))

trying to adapt for x/2 comparisons with the bootstraped statistic.

The function seems to work ok, but in some cases the bootstrap estimates for
the
highest confidence level are obviously wrong:


0.05/x    -1.33318711781158
0.025    -1.09386961626344
0.05    -0.868768454819694
0.1    -0.652649120212949
0.16    0.0453058935697044
0.84    0.104645679223495
0.9    0.213385388775905
0.95    0.368838243723381
0.975    0.615868703042345
1-0.05/x    -1.33318711781158

where the estimate for the highest confidence level is simply the copy of
that for the lowest.

Did someone experience a similar problem?
I guess I should do some debugging on the bcanon code.

Thanks
    Bruno

~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Bruno L. Giordano - Ph. D. student
Dipartimento di Psicologia Generale
Via Venezia 8 - 35131 Padova, Italy

currently hosted by:

Equipe Perception et Cognition Musicales
Ircam-CNRS (UMR 9912)
1 place Igor-Stravinsky
F-75004 Paris, France