Displaying 2 results from an estimated 2 matches for "toy2".
Did you mean:
toy
2011 Feb 24
1
parallel bootstrap linear model on multicore mac (re-post)
...the bootstrapping linear
models example in Venables and Ripley. Please excuse my rather terrible
programming skills. I am always open to suggestions. Below the example I
describe what methods I have tried.
library(boot)
library(lme4)
dat<-read.table("http://www2.fiu.edu/~adick/downloads/toy2.dat", header = T)
nboot<-1000 # number of bootstraps
attach(dat)
x<-dat[,2] # IV number 1
y<-dat[,4] # DV
z<-dat[,3] # IV number 2
subj<-dat[,1] # random factor
boot.fun<-function(data,i) { # function to resample residuals
d<-data
d$y<- d$fit...
2011 Feb 23
0
parallel bootstrap linear model on multicore mac
...suggestions. Below the example I
describe what methods I have tried (which is the reason for all of the
unnecessary library calls).
>library(snow)
library(snowfall)
library(boot)
library(lme4)
library(multicore)
library(rlecuyer)
dat<-read.table("http://www2.fiu.edu/~adick/downloads/toy2.dat", header = T)
nboot<-1000
attach(dat)
x<-dat[,2]
y<-dat[,4]
z<-dat[,3]
subj<-dat[,1]
m2.ph.fun<-function(data,i) {
d<-data
d$y<- d$fitted+d$res[i]
as.numeric(coef(update(m2.fit,data=d))[1][[1]][1,c(1:4)])
}
m2.fit&...