Displaying 1 result from an estimated 1 matches for "mafic_vol".
2012 Apr 30
0
Extracting coefficients values with bootstrap
...ks
There is basically 1 variable (LogRds_25k) and 1 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 ~
L...