Caros colegas, Estou realizando uma análise de variânica do peso vazio inicial de ovelhas em função da idade gestacional e manejo. Após rodar a análise a análise de variância tipo 3 não indicou haver diferenças para o nível do manejo (0,05) entretanto ao utilizar a função glht (multcomp) a análise mostrou haver diferenças entre os pesos. Em qual análise confiar? segue CMR<br>
<br>library(multcomp)<br>library(agricolae)<br>library(doBy)<br>library(phia)<br>library(car)<br><br>analise<-lm(PVINICIAL~Manej,data=GANHO1)<br>Anova(analise,type="II")<br>anova(analise)<br>summary(glht(analise,linfct=mcp(Manej="Tukey",interaction_average=TRUE,covariate_average=TRUE)))<br>
<br><span class="" style="border-collapse:separate;color:rgb(0,0,0);font-family:DotumChe;font-size:13px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:15px;text-indent:0px;text-transform:none;white-space:pre-wrap;word-spacing:0px;background-color:rgb(225,226,229)"><pre tabindex="0" class="" style="font-family:DotumChe;font-size:10pt!important;outline-style:none;border-style:none;white-space:pre-wrap!important;margin:0px;line-height:1.2">
<span class="" style="color:rgb(127,0,85)">dput(GANHO1)
</span>structure(list(Animal = c(25L, 38L, 42L, 50L, 53L, 18L, 22L, 
34L, 41L, 44L, 48L, 2L, 12L, 14L, 39L, 46L, 4L, 5L, 21L, 26L, 
49L, 3L, 32L, 37L, 43L, 47L, 8L, 15L, 30L, 40L, 1L, 17L, 23L, 
28L, 7L, 9L, 19L, 35L), Gest = structure(c(1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 
4L), .Label = c("0", "100", "130", "140"), class = "factor"), 
    Manej = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 
    2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 
    1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), .Label = c("1", 
    "2"), class = "factor"), PVINICIAL = c(36.5, 37, 25.6, 30.3, 
    30.1, 27.3, 31, 25.1, 32.5, 29.6, 35, 40, 34.3, 33.8, 33.7, 
    31.4, 40, 32.3, 34.5, 29.5, 39.7, 38, 33.2, 34.5, 34.8, 31.3, 
    29.2, 32, 30.5, 27.5, 37, 32, 36.5, 40.6, 35.1, 29.8, 36.2, 
    28.6)), .Names = c("Animal", "Gest", "Manej", "PVINICIAL"
), class = "data.frame", row.names = c(NA, -38L))</pre></span><br><br>