5
$\begingroup$

I'm looking to run a linear mixed effect model using lme4, where my dependent variable one_syllable_words / total_words_generated is a proportion and my random effect (1 | participant_ID) reflects the longitudinal nature of the design. Independent, fixed effect variables of interest include age, group, timepoint, and interactions between them.

I've come across two main ways to deal with the proportional nature of the DV:

  1. Standard logistic regression / binomial GLM

    In my scenario, I envision the lme4 equation looking like this: 

    glmer(one_syllable_words / total_words_generated ~ age + group +
timepoint + age:timepoint + age:group + timepoint:group + (1 |
participant_ID), family = "binomial", weights =
total_words_generated, data = mydat)

  2. Beta regression

    I would apply a transformation to my DV (DV * (n - 1) + .5)/ n) so that it cannot equal 0 or 1. (There are a few instances where it equals zero, no instances where it equals one.)

I'm unclear whether logistic regression or beta regression is preferred in this example. My DV isn't a clear-cut case of successes and failures (unless we stretch the definition of "success"), so I'm worried logistic regression might not be appropriate. However, I'm having trouble getting a firm grasp on beta regression & all it entails. If beta regression is preferred:

  1. Why is it preferred?
  2. What is it doing "behind the scenes" to the data?
  3. How can it be applied in R?
Improve this question
$\endgroup$
8
  • 2
    $\begingroup$ +1. My very related Q from a couple of days ago: stats.stackexchange.com/questions/233366 (see my own answer there). $\endgroup$
    – amoeba
    Commented Sep 14, 2016 at 13:52
  • 1
    $\begingroup$ See also stats.stackexchange.com/questions/189115 and stats.stackexchange.com/questions/87956 for how to do the binomial model (which I think is the correct approach in your case) and related issues. Both threads are linked to in my answer that I mentioned above - I tried to make that one comprehensive. $\endgroup$
    – amoeba
    Commented Sep 14, 2016 at 14:08
  • 2
    $\begingroup$ (1) I think that the binomial model is the correct approach because your data is binomial: if there is some underlying probability of a word being one syllable then the number of one syllable words is binomially distributed. This number can be zero, and this cannot be modeled via beta (and your suggested transformation is just an arbitrary hack). (2) This is the same model! Just a different numerical optimizer that fits this model. I had to use it because the glmer function was not converging without it (I have no idea why), and I read that bobyqa is more reliable than the default. $\endgroup$
    – amoeba
    Commented Sep 14, 2016 at 14:27
  • 1
    $\begingroup$ No-no, these are different issues. I could not use my option #1 because I did not have fractions; my situation is different from yours. You have fractions and hence you have a binomial situation. It is still a good idea to model overdispersion though, see BenBolker's answer linked in my second comment here. You should add one more random term into your glmer call, that's it. $\endgroup$
    – amoeba
    Commented Sep 14, 2016 at 14:35
  • 1
    $\begingroup$ Yes, I meant (1|index). With only three observations per subject I am not sure this term will make a lot of difference, but you can still include it and see how it behaves. I think the idea is to include it as a safeguard. Check the variance of this random effect in the fitted model; if it's zero or close to zero then your random subject variance is already enough. Re quasibinomial - you cannot use it in lme4, it does not work. So it's not an option. $\endgroup$
    – amoeba
    Commented Sep 14, 2016 at 15:00

0

Reset to default

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.