[R] Conditional Logistic Regression with Multilevel Data [using clogistic() in Epi and lmer() in lme4???]

[Wolfgang Viechtbauer]

Ouch, that was a brain fart on my part!

But indeed, since we can fit a 'single-site model' using glmer() (as shown), we should be able to just add 'site' as another random effect with (1 | site / id).

Best,
Wolfgang

On July 3, 2019 1:26:54 PM GMT+02:00, Marc Schwartz <marc_schwartz using me.com> wrote:
>Wolfgang,
>
>Perhaps I am missing something here, but with clogit(), or Heather's
>initial attempt at using clogistic(), as both have 'strata' arguments,
>wouldn't the 'strata' argument need to be able to accept a nested
>structure, along the lines of strata(id / site)?
>
>On the page that you link to below, I don't see an example of that type
>of nesting structure.
>
>The mixed effects models would generally support a nested random
>effects specification, if that makes sense here.
>
>Regards,
>
>Marc Schwartz
>
>
>> On Jul 3, 2019, at 4:41 AM, Viechtbauer, Wolfgang (SP)
><wolfgang.viechtbauer using maastrichtuniversity.nl> wrote:
>> 
>> Dear Heather,
>> 
>> You can do this with clogit() from the 'survival' package, using
>strata(). See:
>> 
>>
>http://www.metafor-project.org/doku.php/tips:clogit_paired_binary_data
>> 
>> You can also use lme4::glmer(), glmmML::glmmML(), or
>metafor::rma.glmm() for that. [1]
>> 
>> [1] Rerunning this example, I noticed a slight issue with the
>structuring of the output from rma.glmm() -- if you want to run that
>example, install the 'devel' version of metafor as described here:
>> 
>> https://wviechtb.github.io/metafor/#installation
>> 
>> Best,
>> Wolfgang