Displaying 1 result from an estimated 1 matches for "iffshelf".
2012 Apr 30
0
Extracting coefficients values with bootstrap
...categorical variable
(GeoRock) with 10 categories.
Here is my script:
N = length (data_As[,1])
B = 10000
stor.r2 = rep(0,B)
stor.inter = rep(0,B)
stor.Rds = rep(0,B)
stor.Bimod = rep(0,B)
stor.grano = rep(0,B)
stor.mafic_intru = rep(0,B)
stor.mafic_vol = rep(0,B)
stor.nonmarine = rep(0,B)
stor.iffshelf = rep(0,B)
stor.sedi = rep(0,B)
stor.vol = rep(0,B)
for (i in 1:B){
idx = sample(1:N, replace=T)
newdata = data_As[idx,]
LogRds <- log(newdata$Rds_25k+1)
GeoRock <- factor(newdata$GeoRock)
data_As.boot = lm(newdata$Log_Level ~
LogRds + GeoRock )
stor.r2[i] = summary(data_As.boot...