##       id trial state                  loc            cult pairs ctrl_dismean dis_pres2 dis_pres3
## 1 2012_1     1    MT          Campo_Verde         TMG_803    36         42.5      High      High
## 2 2012_2     2    MT          Pedra_Preta         TMG_803    36         36.3      High    Medium
## 3 2012_3     3    MT   Primavera_do_Leste         TMG_803    36         41.3      High      High
## 4 2012_4     4    MS      Chapadao_do_Sul           5G830    36         30.6      High    Medium
## 5 2012_6     5    MS Sao_Gabriel_do_Oeste BMX_Potencia_RR    36         27.0       Low    Medium
## 6 2012_7     6    MT          Deciolandia         TMG_803    45         44.8      High      High
##   Ry_pres2 Ry_pres3 Ry_max        cor          z         Vz
## 1     1Low     1Low   1.18 -0.1805764 -0.1825785 0.03030303
## 2     1Low  2Medium   1.22 -0.5247035 -0.5828082 0.03030303
## 3    2High  2Medium   1.32 -0.3490060 -0.3643115 0.03030303
## 4     1Low     1Low   1.11 -0.1865706 -0.1887817 0.03030303
## 5     1Low     1Low   1.16 -0.1866246 -0.1888376 0.03030303
## 6     1Low     1Low   1.11 -0.4933235 -0.5404434 0.02380952

Fitted models

Random effects

## 
## Random-Effects Model (k = 26; tau^2 estimator: ML)
## 
##   logLik  deviance       AIC       BIC      AICc  
##  -3.7911   52.3380   11.5822   14.0984   12.1040  
## 
## tau^2 (estimated amount of total heterogeneity): 0.0504 (SE = 0.0219)
## tau (square root of estimated tau^2 value):      0.2245
## I^2 (total heterogeneity / total variability):   64.03%
## H^2 (total variability / sampling variability):  2.78
## 
## Test for Heterogeneity: 
## Q(df = 25) = 72.7184, p-val < .0001
## 
## Model Results:
## 
## estimate       se     zval     pval    ci.lb    ci.ub 
##  -0.5681   0.0551 -10.3122   <.0001  -0.6760  -0.4601

Fixed effects

## 
## Fixed-Effects Model (k = 26)
## 
## Test for Heterogeneity: 
## Q(df = 25) = 72.7184, p-val < .0001
## 
## Model Results:
## 
## estimate       se     zval     pval    ci.lb    ci.ub 
##  -0.5709   0.0330 -17.3081   <.0001  -0.6356  -0.5063

Including a moderator variable (qualitative) (Ry_max)

## 
## Mixed-Effects Model (k = 26; tau^2 estimator: ML)
## 
##   logLik  deviance       AIC       BIC      AICc  
##  -0.2090   45.1739    6.4181   10.1924    7.5090  
## 
## tau^2 (estimated amount of residual heterogeneity):     0.0313 (SE = 0.0166)
## tau (square root of estimated tau^2 value):             0.1769
## I^2 (residual heterogeneity / unaccounted variability): 52.49%
## H^2 (unaccounted variability / sampling variability):   2.10
## 
## Test for Residual Heterogeneity: 
## QE(df = 24) = 54.7823, p-val = 0.0003
## 
## Test of Moderators (coefficient(s) 1,2): 
## QM(df = 2) = 148.9380, p-val < .0001
## 
## Model Results:
## 
##                estimate      se      zval    pval    ci.lb    ci.ub
## Ry_pres21Low    -0.4281  0.0684   -6.2577  <.0001  -0.5622  -0.2940
## Ry_pres22High   -0.7041  0.0672  -10.4775  <.0001  -0.8358  -0.5724

zest zest.inf zest.sup r.zRy_pres2 r.zRy_pres2_inf r.zRy_pres2_sup
Ry_pres21Low -0.43 -0.56 -0.29 -0.40 -0.51 -0.29
Ry_pres22High -0.70 -0.84 -0.57 -0.61 -0.68 -0.52
## 
## Random-Effects Model (k = 13; tau^2 estimator: ML)
## 
## tau^2 (estimated amount of total heterogeneity): 0.0311 (SE = 0.0238)
## tau (square root of estimated tau^2 value):      0.1762
## I^2 (total heterogeneity / total variability):   51.34%
## H^2 (total variability / sampling variability):  2.06
## 
## Test for Heterogeneity: 
## Q(df = 12) = 26.7055, p-val = 0.0085
## 
## Model Results:
## 
## estimate       se     zval     pval    ci.lb    ci.ub 
##  -0.4281   0.0683  -6.2692   <.0001  -0.5619  -0.2942
## 
## Random-Effects Model (k = 13; tau^2 estimator: ML)
## 
## tau^2 (estimated amount of total heterogeneity): 0.0315 (SE = 0.0231)
## tau (square root of estimated tau^2 value):      0.1775
## I^2 (total heterogeneity / total variability):   53.58%
## H^2 (total variability / sampling variability):  2.15
## 
## Test for Heterogeneity: 
## Q(df = 12) = 28.0768, p-val = 0.0054
## 
## Model Results:
## 
## estimate       se     zval     pval    ci.lb    ci.ub 
##  -0.7041   0.0673 -10.4594   <.0001  -0.8361  -0.5722