Crecimiento

Datos

La curva del peso en función de la edad de los pulpos muestra un crecimiento exponencial bajo las dos condiciones, control (naranja) y tratamiento (azul):

Análisis estadísticos

Por lo tanto, se realiza un modelo estadístico lineal sobre el logaritmo del peso, teniendo en cuenta el tratamiento (control / tratamiento), y la edad como variables predictoras. Además, se utiliza una variable aleatoria (cada juvenil), considerando que se midieron los mismos juveniles a las diferentes edades (cada curva representa un juvenil).

El mejor modelo fue seleccionado en base al criterio del mínimo AIC (criterio de información de Akaike) y además teniendo en cuenta que fueran homocedásticos y cumplieran con la normalidad de los supuestos. Esto fue verificado gráficamente y a través de la prueba de Shapiro-Wilks.

Gráfico de residuos estandarizados en función de valores predichos, con colores según el tratamiento. No se observan grandes desviaciones a lo esperado bajo homocedasticidad, pero es curioso que se ve un patrón ondulante en los residuos. Es como si el crecimiento no fuera regular y constante sino como en “pulsos”.

En cuanto a la normalidad, se observa un buen ajuste de los residuos, que además coincide con la prueba de Shapiro-Wilks (W=0.989, p=0.245).

Resultados

El modelo muestra una interacción significativa entre los tratamientos y la edad (Dias:Tipo), que implica que la tasa intrínseca de crecimiento es diferente de acuerdo a la carga parásitaria (control vs tratamiento).

## Analysis of Deviance Table (Type II Wald chisquare tests)
## 
## Response: log(Peso)
##               Chisq Df Pr(>Chisq)    
## Dias      6011.4802  1  < 2.2e-16 ***
## Tipo         0.5834  1      0.445    
## Dias:Tipo   42.0747  1  8.785e-11 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

En relación con la variabilidad entre juveniles, esta representa un 70% de la variación total observada, es decir que la mayor parte de la variabilidad en el peso está explicada por las diferencias entre juveniles, mientras que hay un 30% de variación aleatoria.

Este modelo sugiere que si bien ambos grupos de pulpos crecen exponencialmente, los que están sometidos al tratamiento lo hacen con una tasa intrínseca de crecimiento un 15% menor (azul; 0.0142 1/dias) que los controles (rojo; 0.0167 1/dias).

Supervivencia

Datos

Si se consiedera el día de última medición como la edad de los pulpos al morir, puede analizarse si el tipo de tratamiento influye en la mortalidad de los pulpos. En este gráfico se ve que la probabilidad de supervivencia es levemente mayor en los pulpos control, que en los pulpos tratamiento.

Análisis estadísticos

Para evaluar si existe una dependencia entre la superviencia y el tipo de tratamiento, se corrió un modelo lineal para compar la edad de los pulpos al nacer entre tratamientos.

No se observa heterocedasticidad:

Y en cuanto a la normalidad, también se observa un buen ajuste. Resultados prueba de Shapiro-Wilks (W=0.949, p=0.631).

Resultados

De acuerdo con este modelo, el tipo de tratamiento no afecta a la supervivencia (F=0.253, p=0.627), que es de 244 días en promedio para los dos grupos. Cabe destacar que quizás el n es un poco bajo para detectar diferencias significativas entre los grupos para esta variable.

## Anova Table (Type II tests)
## 
## Response: Dias
##           Sum Sq Df F value Pr(>F)
## Tipo         314  1  0.2534 0.6268
## Residuals  11152  9