[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