Displaying 1 result from an estimated 1 matches for "theparams".
Did you mean:
dhparams
2011 Oct 18
2
Non-linear maximization function in R
...heData) {
x <- theData[ ,1]
y <- theData[ ,2]
logitReg <- glm(y ~ x + I(x^2), family=binomial)
b0 <- logitReg$coefficients[1]
b1 <- logitReg$coefficients[2]
b2 <- logitReg$coefficients[3]
opt <- (-b1)/(2*b2)
tol <- 1 / sqrt(-2*b2)
pmax<-1/(1+exp(b1^2/(4*b2)-b0))
theParams <- cbind(opt, tol, pmax)
theParams
}
>theData<-cbind(Sample,Aequipecten)
>tBLparamsForOneSpecies(theData)
> tBLparamsForOneSpecies(theData)
opt tol pmax
x 0.6337473 0.1468823 0.2433407
#Nor do to the coefficients look similar to p either.
> glm(Aequipecten~...