[R-meta] residual error term in rma.mv

Daniel Noble d@niel@w@@noble @ending from gm@il@com
Thu May 31 00:57:05 CEST 2018


Hi Clint,

You could estimate this by including an observation-level random effect
(~1|obs) after adding an observation-level identifier to each row in your
data frame (i.e., 1:nrow(data).

Cheers,
Dan

–––
Dr. Daniel Noble | ARC DECRA Fellow
Level 5 West, Biological Sciences Building (E26)
Ecology & Evolution Research Centre (E&ERC)
School of Biological, Earth and Environmental Sciences (BEES)
*The University of New South Wales*
Sydney, NSW 2052
AUSTRALIA

T : +61 430 290 053
E : daniel.noble using unsw.edu.au <daniel.noble using mq.edu.au>
W: www.nobledan.com
Github: https://github.com/daniel1noble

On Thu, May 31, 2018 at 3:44 AM, Kelly, Clint Dale <kelly.clint using uqam.ca>
wrote:

>
> A manuscript reviewer recently asked my colleagues and I:
>
> "Were the author’s careful to ensure the residual error term was included
> in their models. rma.mv does not, by default include this term, and it
> needs to be specified.”
>
> It is not clear to me to what they refer and how to “specific this term”.
> Any guidance or insight would be much appreciated.
>
> Here is my model in which species identity (Species), study identity
> (Code) and phylogeny are entered as random factors:
>
> model.1<-rma.mv(yi,vi,random = list(~1|Species,
>                                                             ~1|Code,
>
> ~1|Species.phylo),
>                         R=list(Species.phylo=Species.phyl), struct="UN",
>                         data=meta.data)
>
> Thanks,
> Clint
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-meta-analysis mailing list
> R-sig-meta-analysis using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis
>

	[[alternative HTML version deleted]]



More information about the R-sig-meta-analysis mailing list