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_schwa...@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.viechtba...@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 ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
