[R-meta] Questions about phylogenetic meta-analisys

[Viechtbauer Wolfgang (SP)]

Dear Rafael,

First of all, 'sezf' sounds like these are standard errors. The second argument of rma.mv() (and rma()/rma.uni()) is for the *sampling variances* (the square of the standard errors). So, you should use: rma.mv(zf, sezf^2, ...).

This aside, I cannot diagnose what the problem is without a reproducible example or at least more information. Could you at least post:

unique(h1$species)
class(phylo)
dimnames(phylo)

Best,
Wolfgang

-----Original Message-----
From: Rafael Rios [mailto:biorafaelrm at gmail.com] 
Sent: Wednesday, 27 September, 2017 20:49
To: r-sig-meta-analysis at r-project.org; Viechtbauer Wolfgang (SP)
Subject: Questions about phylogenetic meta-analisys

Dear Dr. Wolfgang Viechtbauer and other meta-analysists,

I am trying to conduct a phylogenetic meta-analysis using function rma.mv with the following format:

meta1=rma.mv(zf, sezf, mods= ~phylum, 
             random = list (~1|study, ~1|species, ~1|field), 
             R=list(species=phylo),
             data = h1)

Where I have a moderator (phylum) and three random variables (study, species and if it was an observational or experimental study). We included phylogenetic correlation matrix (phylo) in argument R. However, I found the following error message: "Error in R[[j]][s.levels[[j]], s.levels[[j]]] : subscript out of bounds". I also tryied to use a matrix with phylogenetic distances, but I found the same error. I hope that this e-mail finds you well.

Best wishes,

Rafael R. Moura.
scientia amabilis

Doutorando da Pós-graduação em Ecologia e Conservação de Recursos Naturais
Universidade Federal de Uberlândia, Uberlândia, MG, Brasil

ORCID: http://orcid.org/0000-0002-7911-4734
Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4244908A8
Research Gate: https://www.researchgate.net/profile/Rafael_Rios_Moura2