Displaying 3 results from an estimated 3 matches for "thetastar".
Did you mean:
theta_star
2004 Apr 09
1
bootstrap function coefficients
...s:
C0
k 0.7101
# Second, I try to use bootstrap to get the 10 estimates of k and C0
> library(bootstrap)
> theta<-function(tabMika){coef(eval(DegSA.nls$call))}
> bootSolA.nls<-bootstrap(tabMika,10,theta)
Warning message:
multi-argument returns are deprecated in: return(thetastar, func.thetastar,
jack.boot.val, jack.boot.se,
> bootSolA.nls
$thetastar
[,1] [,2] [,3] [,4] [,5] [,6]
C0 25.68210358 25.68210358 25.68210358 25.68210358 25.68210358 25.68210358
k 0.08735615 0.08735615 0.08735615 0.08735615 0.08735615 0.0873...
2004 Mar 12
0
Basic questions on nls and bootstrap
...e of course but 10 is retained for this example) of the parameters C0
and t. I use this way:
> library(bootstrap)
> theta<-function(tabMika){coef(eval(DegSA.nls$call))}
> bootSolA.nls<-bootstrap(tabMika,10,theta)
Warning message:
multi-argument returns are deprecated in: return(thetastar, func.thetastar,
jack.boot.val, jack.boot.se,
> bootSolA.nls
$thetastar
[,1] [,2] [,3] [,4] [,5] [,6]
C0 25.68210358 25.68210358 25.68210358 25.68210358 25.68210358 25.68210358
k 0.08735615 0.08735615 0.08735615 0.08735615 0.08735615 0.087...
2004 Apr 27
4
Problems raised to 1/3 power and NaN
...(x)
Sb <- sqrt((sum((x-mean(x))^2)/nb))
Kb <- ((sum((x-mean(x))^3)/(nb*psd^3)))
W <- (bmean-mean)/Sb
Q <- W + (Kb * W^2)/3 + (Kb^2)*(W^3)/27 + Kb/(6*nb)
}
#Bootsrap 10000 values for Q
Q <- bootstrap(x,10000,Qstat,pmean)
#Determine the quantile of Q to use
Qa <- quantile(Q$thetastar, 0.05, na.rm = TRUE, names = FALSE)
#Calulate Wq*******Here is the code tha does not work***************
Wq <- (3/k)*((1+k*(Qa-(k/(6*n))))^(1/3)-1)
#Halls bootstrap
UCL <- pmean-(Wq*psd)
UCL
}
Michael J. Bock, PhD.
ARCADIS
24 Preble St. Suite 100
Portland, ME 04101
207.828.0046
fax 2...