Displaying 1 result from an estimated 1 matches for "nlmfunc1".
Did you mean:
lmfunc
2011 Mar 19
0
Problems using NLS in conjunction with non-parameteric bootstrapping
...5)
mydata <- data.frame(x,y)
names(mydata) <- c("x","y") # variable names
%x=Trap success; proportion of total trap nights where at least one deer was
photographed
%y=Square root transformation of total number of deer photographed
#asymptotic equation y = (a*x)/(b+x)
NlMFunc1 <- function(x, a, b) {(a*x)/(b+x)}
NlMFunc1 <- nls(y ~ NlMFunc1(x, alpha, beta),
data = RPOvTS, start = list(alpha = 10, beta = 1), control = list(maxiter =
10000, minFactor=1/20048),
trace = T)
summary(NlMFunc1)
logLik(NlMFunc1)
#Generate non-parametric bootstrap estimates of alpha and b...