Efecto global del tratamiento a partir de un modelo mixto en metafor

Question

Esto es la continuación de mi post anterior ( Prueba del efecto de los moderadores en el paquete metafor ) Gracias@Wolfgang por la aclaración. Ahora he utilizado el modelo: A<-rma.mv(yi, vi, mods = ~ factor(suelo)+factor(cultivo)+ ppt, random=~1|id/año-estudio, method = "REML", data = a) [Ahora también he añadido la precipitación (ppt); como variable continua]

Mis preguntas son:

1. ¿Cómo se extrae el efecto global del tratamiento de este modelo? (He extraído el efecto global eliminando todos los moderadores, pero me pregunto si podría obtener el efecto global de este modelo incluyendo los moderadores).

2. También el intercepto incluye tanto un nivel de suelo como un nivel de cultivo como referencia, estoy confundido sobre cómo comparar otros niveles de suelo con el intercepto ya que el intercepto también está actuando como referencia para el factor de cultivo al mismo tiempo. ¿Hay alguna manera de separar estos dos? Cuando elimino el cultivo del modelo, el valor del intercepto cambia y se vuelve significativo (en este caso, ¿significa que ese nivel de suelo es significativamente diferente de cero?) Quiero mirar los valores medios de cada uno de los niveles de los factores y si difieren significativamente de cero. ¿Cómo puedo hacerlo?

3. He nombrado las identificaciones individuales de cada estudio en el conjunto de datos, pero un par de estudios se realizan en realidad en los mismos sitios de campo. Entonces, ¿debería el "id" de mi modelo ser el estudio o el emplazamiento (y, por tanto, tengo que renombrar los ids en el conjunto de datos según los emplazamientos en lugar de según los estudios)?

4. El número total de estudios en mi conjunto de datos es de 30, lo que parece estar en el límite. Entonces, si el efecto de cualquier nivel de factor resulta no significativo, ¿significa que realmente no es significativo o que se debe a la falta de datos suficientes? En otras palabras, ¿cómo podemos calcular la potencia para sacar conclusiones? (No significación = falta de potencia frente a no significación real)? También estoy probando otra variable de respuesta para la que el número de estudios es sólo de 8. ¿Se trata de la misma pregunta o debo realizar el análisis para esta variable?

Les agradecería mucho que me orientaran al respecto. Gracias.

Answer 1

P1 - Si muestra un efecto del moderador, entonces no hay un efecto global. Sería como estimar la altura media de los humanos cuando se ha demostrado que los hombres son más altos que las mujeres.

P2: el coeficiente para el cultivo se aplica independientemente del nivel del suelo y viceversa, a menos que se ajuste una interacción.

P3: con el tipo de cifras que cita, sospecho que si intenta ajustar efectos aleatorios a nivel de centro y a nivel de estudio, uno de ellos se basará en pocos datos y será bastante inestable, por lo que me decantaría por un efecto de estudio.

P4 - Con 30 observaciones no querría utilizar más de 2 grados de libertad para mis moderadores o quizás 3 si me presionan.