Estou tentando determinar se as lagartas que comem uma dieta natural (monkeyflower) são mais resistentes aos predadores (formigas) do que as lagartas que comem uma dieta artificial (uma mistura de gérmen de trigo e vitaminas). Fiz um estudo experimental com uma amostra pequena (20 lagartas; 10 por dieta). Eu pesava cada lagarta antes do experimento. Ofereci um par de lagartas (uma por dieta) a um grupo de formigas por um período de cinco minutos e contei o número de vezes que cada lagarta foi rejeitada. Repeti esse processo dez vezes.
É assim que meus dados são (A = dieta artificial, N = dieta natural):
Trial A_Weight N_Weight A_Rejections N_Rejections
1 0.0496 0.1857 0 1
2 0.0324 0.1112 0 2
3 0.0291 0.3011 0 2
4 0.0247 0.2066 0 3
5 0.0394 0.1448 3 1
6 0.0641 0.0838 1 3
7 0.0360 0.1963 0 2
8 0.0243 0.145 0 3
9 0.0682 0.1519 0 3
10 0.0225 0.1571 1 0
Estou tentando executar uma ANOVA em R. É assim que meu código se parece (0 = Dieta artificial, 1 = Dieta natural; todos os vetores são organizados com dados para as dez lagartas artificiais da dieta primeiro, seguidos pelos dados da dieta dez natural lagartas):
diet <- factor (rep (c (0, 1), each = 10)
rejections <- c(0,0,0,0,3,1,0,0,0,1,1,2,2,3,1,3,2,3,3,0)
weight <- c(0.0496,0.0324,0.0291,0.0247,0.0394,0.0641,0.036,0.0243,0.0682,0.0225,0.1857,0.1112,0.3011,0.2066,0.1448,0.0838,0.1963,0.145,0.1519,0.1571)
all.data <- data.frame(Diet=diet, Rejections = rejections, Weight = weight)
fit.all <- lm(Rejections ~ Diet * Weight, all.data)
anova(fit.all)
E são estes os meus resultados:
Analysis of Variance Table
Response: Rejections
Df Sum Sq Mean Sq F value Pr(>F)
Diet 1 11.2500 11.2500 9.8044 0.006444 **
Weight 1 0.0661 0.0661 0.0576 0.813432
Diet:Weight 1 0.0748 0.0748 0.0652 0.801678
Residuals 16 18.3591 1.1474
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Minhas perguntas são:
- A ANOVA é apropriada aqui? Sei que o pequeno tamanho da amostra seria um problema em qualquer teste estatístico; este é apenas um estudo de avaliação no qual eu gostaria de exibir estatísticas para uma apresentação em classe. Eu pretendo refazer este estudo com um tamanho de amostra maior.
- Eu inseri meus dados no R corretamente?
- Isso está me dizendo que a dieta é significativa, mas o peso não?