R-beta: Result_abcnon

[Douglas Steele]

In message <35B7E998.5BC6 at kku1.kku.ac.th>, Nikom Thanomsieng
<nikom at kku1.kku.ac.th> writes
>Dear, Sir.
>       I run R program (rw0613b) under Windows 95  and used Functions for the
>Book 
>"An Introduction to the Bootstrap", that help in html files to compute
>nonparametric
>abc confidence intervals for the mean. I write command following
>"
>> library(bootstra)
>> library(help = bootstra)
>Read 8 items
>abcnon         Nonparametric ABC Confidence Limits 
>abcpar         Parametric ABC Confidence Limits  
>bcanon         Nonparametric BCa Confidence Limits 
>bootpred       Bootstrap Estimates of Prediction Error 
>bootstrap      Non-parametric Bootstrapping 
>boott          Bootstrap-t Confidence Limits 
>crossval       K-fold Cross-Validation 
>jackknife      Jackknife Estimation 
>> x <- rnorm(10)
>> x
> [1]  0.1808591  0.4303385  1.4218043  0.7222164 -0.9713840  0.6121021
> [7]  0.0555533 -0.9721635  0.5111223 -2.0564528
>> theta <- function(p,x) {sum(p*x)/sum(p)}
>> results <- abcnon(x, theta)
>> 
>"
>        The command that I wrote can not show the output of abc
>intervals for
>the mean.
>       This command correct?
What do you get when you just type results (or results$limits).  For
example I get...
> results
$limits
     alpha         abc        stan
[1,] 0.025 -0.30521640 -0.25558396
[2,] 0.050 -0.22529440 -0.18902905
[3,] 0.100 -0.13625620 -0.11229551
[4,] 0.160 -0.06816337 -0.05165805
[5,] 0.840  0.35371418  0.36842367
[6,] 0.900  0.40806807  0.42906113
[7,] 0.950  0.47470505  0.50579467
[8,] 0.975  0.53062840  0.57234958

If you are not getting anything at all, have you got the library
properly installed ?  (Why bootstra instead of bootstrap since you are
running under Win95 so no 8 character truncation required and the
origional binaries are named bootstrap ?)




-- 
Douglas Steele
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._