[R] Stratified Bootstrap question
Tim Hesterberg
timh at insightful.com
Wed Mar 30 19:05:45 CEST 2005
Dear Qian,
You might try the S+Resample library, which has built-in support
for both sampling by subject and stratified sampling.
If you are a student, there is a free student version of S+.
See
www.insightful.com/downloads/libraries (S+Resample)
www.insightful.com/Hesterberg/bootstrap (has link to the student version)
For the missing values, consider the S+Missing library,
which offers multiple imputation. With S+, do
library(missing)
Tim Hesterberg
P.S. The combination of sampling by subject and stratified sampling
was terribly messy to program. If I'd known in advance how messy, I
never would have done it :-( But it is done now.
>Dear R users,
>
>I have a question regarding stratified bootstrap question and how to implement
>it using boot() in R's boot package.
>
>My dataset is a longitudinal dataset (3 measurements per person at year
>1, 4 and 7) composed of multiple clinic centers and multiple participants
>within each clinic. It has missing values.
>
>I want to do a bootstrap to find the standard errors and confidence
>intervals for my variance components. My model is a mixed model with
>random clinic and random participant within clinic.
>
>I thought two methods to do bootstrap:
>(1) bootstrap data; however, I have problem specifying the second
>parameter for my statistic function, shall I use indices, weight or
>frequency and how shall I relate to my dataset.
>(2) bootstrap residuals; however, the dataset has multiple measurements
>and missing values. I am wondering how to construct a new data frame
>containing the residuals and fitted values.
>
>Any ideas will be highly appreciated!
>Sincerely yours,
>Qian
========================================================
| Tim Hesterberg Research Scientist |
| timh at insightful.com Insightful Corp. |
| (206)802-2319 1700 Westlake Ave. N, Suite 500 |
| (206)283-8691 (fax) Seattle, WA 98109-3044, U.S.A. |
| www.insightful.com/Hesterberg |
========================================================
Download the S+Resample library from www.insightful.com/downloads/libraries
More information about the R-help
mailing list