Displaying 1 result from an estimated 1 matches for "lithologyb".
Did you mean:
lithology
2009 Apr 30
0
Categorical variable in a custom nonlin function with gnm
...el. I have previously done something similar with nonlinear least squares models (nls) and it worked very well, however gnm doesn't seem to like these dummy variables, and doesn't estimate the parameters.
It works if I put the dummy variables into a linear model using gnm, like so:
(where LithologyB is a series of ones and zeros, representing the occurrence of Lithology B, etc.)
> summary(gnm(response~LithologyB+LithologyC+LithologyD+LithologyE+LithologyF, family=binomial))
Linear predictor - using glm.fit
Call:
gnm(formula = response ~ LithologyB + LithologyC + LithologyD + Lithology...