Pessoal; Necessito da grande ajuda de vcs.... Tem um delineamento em esquema fatorial, mas os dados estão desbalaceados. Nesse caso, um dos fatores não está completo.... Eu tenho 4 clones, 3 locais e 2 espaçamentos, mas para um dos clones eu não tenho todos os locais. Como posso proceder com essa estatística no R? O Walmes já me deu alguns toques sobre algo parecido a algum tempo.... Conto com vcs! Muito obrigado! --- Em sáb, 12/2/11, walmes zeviani <walmeszeviani@yahoo.com.br> escreveu: De: walmes zeviani <walmeszeviani@yahoo.com.br> Assunto: Res: Res: [R_STAT] Dados não balanceados! Para: R_STAT@yahoogrupos.com.br Data: Sábado, 12 de Fevereiro de 2011, 13:21   Thiago, Como o Ivan disse, a coisa do teste dos termos do modelo ou anova, pode ser empregada a car::Anova(..., type="III"). Mas os cuidados com desbalanceamentos atravessam a fronteira da anova. Os testes de médias (ou demais inferências) não são mais baseados nas médias amostrais, uma vez que essas são viesadas devido a perda de ortogonalidade causada pelo desbalanceamento. O que deve ser obtido são as médias ajustadas que requerem manipulação matricial do vetor de solução do sistema linear, ou coef(modelo.lm). Antes de tudo isso, deve se saber se o seu desbalanceamento não comprometeu a estimação dos efeitos, pois existem duas "intensidades" de desbalanceamento: * apenas perda de repetições: aqui toda combinação entre níveis dos fatores tem pelo menos uma repetição com observação, e.g. fator A (A1,A2,A3) e B (B1,B2) todas as 6 combinações de A com B possuem pelo menos uma observação. Aqui todos os efeitos podem ser estimados embora com erros padrões diferentes em função do desbalanceamento. * perda de caselas: aqui pelo menos uma combinação dos níveis dos tratamentos não foi observada, e.g. fator A (A1,A2,A3) e B (B1,B2) alguma das 6 combinações foi completamente não observada (e.g. A1B1). Nesse caso não é possível estimar alguns efeitos justamente por não haver as observações experimentais correspondentes para estimá-los. No caso de dois fatores sem interação entre eles tem um exemplo na página de Rídiculas do LEG que fiz ontem para outro fim. Quem sabe possamos fazer uma Rídicula com o seu problema. http://www.leg.ufpr.br/doku.php/ridiculas À disposição. Walmes.   ============================================================ Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br / @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes ============================================================ De: Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Para: R_STAT@yahoogrupos.com.br Enviadas: Sábado, 12 de Fevereiro de 2011 10:38:53 Assunto: Res: [R_STAT] Dados não balanceados!   Bom dia Thiago!! Nestes casos é aconselhável vc utilizar a soma de quadrados do tipo III. Para obter resultados semelhantes a alguns softwares consagrados vc pode  utilizar a reparametrização do tipo soma zero. Para isso: options(contrasts = c("contr.sum", "contr.helmert")) library(car) modelo <- lm( y ~ A*B) Anova(modelo, type="III") Abraço. (S,f,P) Allaman   M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman@yahoo.com.br Tel: (35)3826-6608/9925-6428 De: Thiago De paula protásio <depaulaprotasio@yahoo.com.br> Para: R_STAT@yahoogrupos.com.br Enviadas: Sábado, 12 de Fevereiro de 2011 10:20:28 Assunto: [R_STAT] Dados não balanceados!   Olá Pessoal; Estou precisando de algumas orientações para rodar um um fatorial com dados desbalaceados. Espero que alguém posso me ajudar com alguma função ou material. Muito obrigado! Bom final de semana a todos!!