[R-sig-ME] Mixes effect with polynomial term

Thierry Onkelinx thierry.onkelinx at inbo.be
Tue Jan 24 10:45:26 CET 2017


Dear Ratana,

Make a plot of the fitted value using only the fixed effect. That will be
easier to interpret that the raw coefficients.

Best regards,

ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
Forest
team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
Kliniekstraat 25
1070 Anderlecht
Belgium

To call in the statistician after the experiment is done may be no more
than asking him to perform a post-mortem examination: he may be able to say
what the experiment died of. ~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data. ~ Roger Brinner
The combination of some data and an aching desire for an answer does not
ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey

2017-01-24 10:26 GMT+01:00 Pin chanratana <chanratana.pin op gmail.com>:

> Hello,
>
> I'm fitting data of Lesser adjutant that were camera-trapped at waterholes
> during the dry season.
> My final model consist of depth of waterholes and the depth power 2. Can
> anyone
> help me in interpreting the result shown below.
>
> glmmadmb ( LA ~ depth + I(depth^2) + offset(log(trap)) + (1 | obs), data =
> ndata4, family = "nbinom")
>
> Coefficients:
>             Estimate Std. Error z value Pr(>|z|)
> (Intercept)  -3.0644     0.2397  -12.78  < 2e-16 ***
> depth         1.5807     0.2633    6.00  1.9e-09 ***
> I(depth^2)   -0.3598     0.0737   -4.89  1.0e-06 ***
>
> Number of observations: total=294, obs=9
> Random effect variance(s):
> Group=obs
>             Variance StdDev
> (Intercept)   0.3244 0.5696
>
> Negative binomial dispersion parameter: 0.34381 (std. err.: 0.064328)
>
>
> Thanks,
>
> Ratana
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models op r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>

	[[alternative HTML version deleted]]



More information about the R-sig-mixed-models mailing list