Displaying 1 result from an estimated 1 matches for "lnpri25".
2012 Apr 30
3
95% confidence interval of the coefficients from a bootstrap analysis
...riation of my variables to
replacement.
I am trying to obtain the coefficients 95% confidence interval from the
bootstrap procedure.
Here is my script for the bootstrap:
N = length (data_Pb[,1])
B = 10000
stor.r2 = rep(0,B)
stor.r2 = rep(0,B)
stor.inter = rep(0,B)
stor.Ind5 = rep(0,B)
stor.LNPRI25 = rep(0,B)
stor.NPRI10 = rep(0,B)
stor.Mine = rep(0,B)
for (i in 1:B){
idx = sample(1:N, replace=T)
newdata = data_Pb[idx,]
L_NPRI_25k <- log(newdata$NPRI_25k+1)
data_Pb.boot = lm(newdata$Log_Level ~
newdata$Ind_5k + newdata$MineP_50k +
newdata$NPRI_10k + L_NPRI_25k )
stor.r2[i...