[R-sig-ME] Logistic exposure glmer polynomial covariate issue
Ben Bolker
bbo|ker @end|ng |rom gm@||@com
Mon May 17 21:06:10 CEST 2021
I worked on this a little more, posted at
https://github.com/bbolker/mixedmodels-misc/blob/master/notes/rsm_logexp.R
I tried the power-logistic link and got it to work with linear (but
not quadratic) models: it does appear to do a bit better on this data
set, although with the small sample I wouldn't take that too seriously.
(I left the predictor as G250, not log(1+G250).)
I would certainly try some of this stuff out on your full data set
(but be careful not to fish/snoop too much ...)
On 5/17/21 1:43 PM, Kirsti Carr wrote:
> Hi Ben,
>
> Thank you for your thoughts. I will look into cloglog and your other
> suggestions. I should note that the data I provided was not the full
> dataset (after reading the recommendations for posting I decided to provide
> a minimum dataset that still produced the error). Actual effective sample
> size/ number of failures is 23. Still small -- so perhaps taking out the
> cubic and quadratic models is best? I can also provide full dataset if
> necessary.
>
> Thanks so much,
> Kirsti
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models using r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
More information about the R-sig-mixed-models
mailing list