[R] Bootstrap and Jackknife Bias using Survey Package
Ramón Casero Cañas
8-T at gmx.net
Tue Apr 11 13:25:10 CEST 2006
Carlos Creva Singano (M2004078) wrote:
>
> 1. How to compute Bootstrap and Jackknife Bias of estimates, like mean?
Have you had a look at packages "boot" and "bootstrap"? E.g. you can
compute the bias and s.e. of an estimate theta using bootstrap
library("boot")
a <- boot( data, theta, R=1000 )
where boot is a function of the boot package, and theta is a function
that you can define yourself. In your case, it is not necessary and you
can use mean instead of theta.
> 2. How to see each replicate estimate of the parameter?
In the example above, they can be found in the boot object "a" in:
a$t
--
Ramón Casero Cañas
http://www.robots.ox.ac.uk/~rcasero/wiki
http://www.robots.ox.ac.uk/~rcasero/blog
More information about the R-help
mailing list