Analise de experimentos em parcelas subdivididas

Ver o tópico anterior Ver o tópico seguinte Ir em baixo

Analise de experimentos em parcelas subdivididas

Mensagem  marciel lelis em Qua Jun 27, 2012 1:09 pm

Estou tentando fazer uma análise de um experimento em parcelas subdivididas, mas não consigo. Importo os dados escrevo o modelo pra poder rodar todos os dados de uma so vez mas tá dando erro. O banco de dados segue abaixo eu estou estudando doses de Zn em cinco clones.

Erro em if (projections) qr <- lmcall$qr <- TRUE :
argumento não é interpretável como lógico
Além disso: Mensagens de aviso perdidas:
In if (projections) qr <- lmcall$qr <- TRUE :
a condição tem comprimento > 1 e somente o primeiro elemento será usado


bloc Zn clon Prod enriaz CV REIT MS.PA MS.R MS.TOTAL
1 1 56 245 3.9375 83.3333 100.0000 66.6667 0.3852 0.1031 0.4883
2 1 56 429 3.6667 85.4167 89.5833 67.9487 0.2621 0.0638 0.3260
3 1 56 498 3.6250 83.3333 95.8333 69.0476 0.2974 0.0633 0.3607
4 1 56 523 4.2708 91.6667 97.9167 78.6111 0.3315 0.0748 0.4063
5 1 56 1213 4.1250 95.8333 97.9167 81.2500 0.3099 0.0935 0.4034
6 2 56 245 4.2708 87.5000 100.0000 72.9167 0.2881 0.0830 0.3711
7 2 56 429 3.5417 85.4167 97.9167 70.8333 0.2953 0.0756 0.3709
8 2 56 498 2.8750 91.6667 100.0000 66.6667 0.3319 0.0823 0.4141
9 2 56 523 3.7917 95.8333 97.9167 85.4167 0.3862 0.0884 0.4746
10 2 56 1213 3.9583 97.9167 100.0000 81.2500 0.3573 0.1209 0.4782
11 3 56 245 4.1667 91.6667 100.0000 66.6667 0.3356 0.1089 0.4445
12 3 56 429 4.1250 93.7500 100.0000 85.4167 0.3856 0.1011 0.4867
13 3 56 498 3.8750 91.6667 91.6667 68.7500 0.3638 0.0891 0.4528
14 3 56 523 3.3125 89.5833 97.9167 68.7500 0.3005 0.0944 0.3950
15 3 56 1213 4.4583 89.5833 95.8333 74.7222 0.3878 0.1227 0.5105
16 1 19 245 3.7292 87.5000 95.8333 77.0833 0.3470 0.0777 0.4247
17 1 19 429 3.3958 93.7500 95.8333 81.1111 0.2569 0.0591 0.3160
18 1 19 498 2.8958 95.8333 100.0000 66.6667 0.3092 0.0569 0.3661
19 1 19 523 3.3542 83.3333 91.6667 73.2143 0.4396 0.0823 0.5218
20 1 19 1213 3.7708 95.8333 100.0000 89.5833 0.3185 0.1106 0.4291
21 2 19 245 4.3542 93.7500 95.8333 75.0000 0.3635 0.0918 0.4553
22 2 19 429 4.0833 89.5833 93.7500 77.0833 0.3135 0.0977 0.4112
23 2 19 498 3.4792 97.9167 100.0000 87.5000 0.3688 0.0913 0.4601
24 2 19 523 3.5625 95.8333 97.9167 60.4167 0.3402 0.0716 0.4118
25 2 19 1213 4.6250 95.8333 97.9167 85.4167 0.3580 0.1395 0.4975
26 3 19 245 3.6875 83.3333 100.0000 68.7500 0.3000 0.0866 0.3866
27 3 19 429 3.9583 83.3333 87.5000 72.3214 0.3216 0.0773 0.3989
28 3 19 498 3.2917 93.7500 97.9167 60.4167 0.3210 0.0823 0.4033
29 3 19 523 2.6875 89.5833 97.9167 66.6667 0.3376 0.0708 0.4084
30 3 19 1213 4.6042 91.6667 100.0000 79.1667 0.3767 0.1406 0.5172
31 1 37 245 4.2917 83.3333 97.9167 68.7500 0.3858 0.0908 0.4766
32 1 37 429 3.8958 91.6667 95.8333 77.0833 0.3343 0.0758 0.4100
33 1 37 498 3.1250 91.6667 100.0000 81.2500 0.2952 0.0778 0.3730
34 1 37 523 3.4167 83.3333 89.5833 81.2500 0.3861 0.0758 0.4620
35 1 37 1213 4.1875 87.5000 89.5833 85.2778 0.3268 0.1051 0.4319
36 3 37 245 4.3542 91.6667 100.0000 64.5833 0.3506 0.0978 0.4484
37 3 37 429 3.6875 83.3333 97.9167 72.9167 0.3241 0.0846 0.4087
38 3 37 498 3.7292 100.0000 100.0000 77.0833 0.3485 0.0848 0.4333
39 3 37 523 3.2500 91.6667 95.8333 75.0000 0.3866 0.0894 0.4760
40 3 37 1213 4.2083 95.8333 95.8333 82.1429 0.3522 0.1143 0.4665
41 2 37 245 3.6875 89.5833 97.9167 59.7222 0.3567 0.0918 0.4485
42 2 37 429 3.1250 60.4167 75.0000 47.5962 0.2920 0.0886 0.3806
43 2 37 498 3.4167 89.5833 97.9167 85.4167 0.3146 0.0819 0.3965
44 2 37 523 3.1875 91.6667 93.7500 66.6667 0.3949 0.0928 0.4877
45 2 37 1213 3.5833 91.6667 95.8333 62.5000 0.3883 0.1441 0.5324
46 1 84 245 4.1042 87.5000 100.0000 68.7500 0.3671 0.0834 0.4505
47 1 84 429 3.6250 83.3333 89.5833 77.0833 0.3823 0.0815 0.4639
48 1 84 498 2.8125 93.7500 95.8333 80.9722 0.3241 0.0733 0.3974
49 1 84 523 3.0625 95.8333 97.9167 85.4167 0.3589 0.0668 0.4258
50 1 84 1213 4.1875 93.7500 100.0000 89.5833 0.3227 0.1069 0.4296
51 2 84 245 4.9583 83.3333 100.0000 64.5833 0.4215 0.1054 0.5268
52 2 84 429 3.4583 93.7500 97.9167 83.3333 0.3227 0.0898 0.4125
53 2 84 498 3.1250 91.6667 97.9167 66.6667 0.3385 0.0719 0.4104
54 2 84 523 3.6875 91.6667 100.0000 60.4167 0.3417 0.0732 0.4149
55 2 84 1213 4.3542 91.6667 97.9167 87.5000 0.3347 0.1398 0.4744
56 3 84 245 3.6875 56.2500 87.5000 35.4167 0.3847 0.1661 0.5508
57 3 84 429 3.3750 83.3333 100.0000 50.0000 0.2973 0.0688 0.3661
58 3 84 498 3.3542 91.6667 95.8333 61.8056 0.3192 0.0728 0.3920
59 3 84 523 2.5833 89.5833 100.0000 52.0833 0.3921 0.0884 0.4805
60 3 84 1213 3.6667 85.4167 97.9167 66.1111 0.3437 0.1118 0.4554
61 1 112 245 3.9167 89.5833 97.9167 77.0833 0.3487 0.0778 0.4265
62 1 112 429 3.6667 91.6667 95.8333 81.2500 0.3647 0.0841 0.4488
63 1 112 498 3.3125 93.7500 95.8333 80.6944 0.3700 0.0855 0.4555
64 1 112 523 3.3542 91.6667 95.8333 77.0833 0.3854 0.0794 0.4648
65 1 112 1213 4.3125 85.4167 97.9167 68.7500 0.3620 0.1121 0.4741
66 2 112 245 4.1458 95.8333 100.0000 50.0000 0.3750 0.0778 0.4528
67 2 112 429 4.2917 89.5833 91.6667 81.2500 0.3255 0.0788 0.4043
68 2 112 498 3.3958 95.8333 97.9167 70.8333 0.3546 0.0790 0.4336
69 2 112 523 3.6458 95.8333 97.9167 64.5833 0.3704 0.0830 0.4534
70 2 112 1213 4.4375 97.9167 100.0000 85.4167 0.3579 0.1244 0.4824
71 3 112 245 4.2500 91.6667 97.9167 62.5000 0.3454 0.1011 0.4465
72 3 112 429 3.4375 87.5000 93.7500 55.9722 0.2893 0.0661 0.3554
73 3 112 498 3.5625 93.7500 97.9167 74.7222 0.3308 0.0808 0.4116
74 3 112 523 3.0417 87.5000 95.8333 64.1667 0.3837 0.0884 0.4721
75 3 112 1213 3.8958 85.4167 93.7500 70.8333 0.3590 0.1275 0.4865

marciel lelis

Mensagens : 4
Data de inscrição : 27/06/2012

Ver perfil do usuário

Voltar ao Topo Ir em baixo

Re: Analise de experimentos em parcelas subdivididas

Mensagem  gustavoreis em Qui Jun 28, 2012 10:00 pm

Boa noite marciel!

Para que eu possa ajudá-lo preciso que você me informe os comandos que está utilizando.

att

Gustavo Mello Reis

gustavoreis
Admin

Mensagens : 70
Data de inscrição : 01/03/2011
Idade : 31

Ver perfil do usuário http://estatisticanor.forumeiros.com

Voltar ao Topo Ir em baixo

Re: Analise de experimentos em parcelas subdivididas

Mensagem  marciel lelis em Sex Jun 29, 2012 4:12 pm

Boa Tarde Gustavo!
Este foram os comandos que usei.

dados<-read.table("Zinco.txt",head=TRUE)
dados
dados <- data.frame(dados)
summary(dados)
variaveis <- dados[,-c(1:3)]

#ANOVA

modelo <- list()
for(i in 1:ncol (variaveis)){
modelo[[i]] <- aov(variaveis[[i]]~bloc+Zn+clon+Zn*clon+Error(bloc/Zn),dados, data=dados)
}
lapply(modelo,anova)

marciel lelis

Mensagens : 4
Data de inscrição : 27/06/2012

Ver perfil do usuário

Voltar ao Topo Ir em baixo

Re: Analise de experimentos em parcelas subdivididas

Mensagem  gustavoreis em Sab Jun 30, 2012 12:15 am

Boa noite Marciel!

Tente os seguintes comandos:

dados<-read.table("Zinco.txt",head=TRUE)
dados
dados <- data.frame(dados)
summary(dados)
variaveis <- dados[,-c(1:3)]

#ANOVA

modelo <- list()
for(i in 1:ncol (variaveis)){
modelo[[i]] <- aov(variaveis[[i]]~bloc+Zn+clon+Zn*clon+Error(bloc/Zn),dados)
}
lapply(modelo,summary)

Se continuar dando erro poste aqui, caso dê certo confirme.

att,

Gustavo Mello Reis

gustavoreis
Admin

Mensagens : 70
Data de inscrição : 01/03/2011
Idade : 31

Ver perfil do usuário http://estatisticanor.forumeiros.com

Voltar ao Topo Ir em baixo

Re: Analise de experimentos em parcelas subdivididas

Mensagem  marciel lelis em Seg Jul 02, 2012 9:39 am

Bom dia Gustavo!

Muito Obrigado. Os comandos que vc me passo deram certo.

Att.

Marciel.

marciel lelis

Mensagens : 4
Data de inscrição : 27/06/2012

Ver perfil do usuário

Voltar ao Topo Ir em baixo

Ainda sobre parcelas subdivididas...

Mensagem  Caroline.Rodrigues em Seg Mar 25, 2013 8:38 pm

Boa Noite pessoal! Também preciso de uma ajuda!

Estou tentando realizar a ANOVA de um experimento muito parecido com o que o Marciel trabalhou: DBC em parcelas subdivididas com 3 blocos, 5 níveis para parcela e 2 níveis para subparcela. Consigo executar os comandos no R porém o resultado da ANOVA não está correto. Acredito que seja algum erro de comando ou até mesmo no banco de dados que ainda não consegui identificar....
Os comandos que estou utilizando são:

dados<-read.table("Anova_241.txt",header=T)
> dados
Bloco A B VSCC
1 1 5 1 4.496
2 1 5 2 68.278
3 1 4 1 3.022
4 1 4 2 65.507
5 1 3 1 5.381
6 1 3 2 69.674
7 1 2 1 4.949
8 1 2 2 64.767
9 1 1 1 1.702
10 1 1 2 51.276
11 2 2 1 4.140
12 2 2 2 66.197
13 2 4 1 3.863
14 2 4 2 65.020
15 2 5 1 6.853
16 2 5 2 75.037
17 2 1 1 6.761
18 2 1 2 71.209
19 2 3 1 3.518
20 2 3 2 58.034
21 3 1 1 4.321
22 3 1 2 71.324
23 3 3 1 4.284
24 3 3 2 71.352
25 3 2 1 3.944
26 3 2 2 71.331
27 3 5 1 3.190
28 3 5 2 65.138
29 3 4 1 -0.104
30 3 4 2 46.184
> saida<-aov(VSCC~Bloco+A+B+A*B+Error(Bloco/A),dados)
> summary(saida)

Error: Bloco
Df Sum Sq Mean Sq
Bloco 1 0.1828 0.1828

Error: Bloco:A
Df Sum Sq Mean Sq
A 1 23.04 23.04

Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
A 1 146 146 4.635 0.0416 *
B 1 28214 28214 895.505 <2e-16 ***
A:B 1 1 1 0.022 0.8834
Residuals 24 756 32
---
Signif. codes: 0 ‘**

Caroline.Rodrigues

Mensagens : 1
Data de inscrição : 24/03/2013

Ver perfil do usuário

Voltar ao Topo Ir em baixo

Re: Analise de experimentos em parcelas subdivididas

Mensagem  Conteúdo patrocinado


Conteúdo patrocinado


Voltar ao Topo Ir em baixo

Ver o tópico anterior Ver o tópico seguinte Voltar ao Topo

- Tópicos similares

 
Permissão deste fórum:
Você não pode responder aos tópicos neste fórum