Cesar Rabak e demais colegas, boa tarde.
Agradeço o retorno e peço desculpas pela demora em responder esta postagem.
Respondendo os questionamentos:
- Por que então que você constrói (especifica) o seu termo Error como PARC:REP?
R: o modelo da anova que eu utilizei [mod<-aov(var~PARC*SubP+Error(PARC:REP))]
retirei de outras postagens feitas no fórum e também o vi em outras consultas,
portanto, não sei o porquê do "Error como PARC:REP" (e gostaria de entender),
mas entendo que seria o erro referente às parcelas ('Erro a' do resultado da
anova). Em todas as minhas consultas o modelo para parcela subdividida é sempre
este, mas comigo não funcionou (ou melhor, funcionou, mas apareceu uma mensagem
de erro, como mostrado anteriormente). Entendo que se eu não inserir o termo do
'Error', a análise será gerada considerando um esquema fatorial.
- Ou fazendo uma pergunta mais próxima da teoria da análise que você quer fazer
(uma ANOVA com medidas repetidas), qual é o erro "dentro" que você quer levar
em conta na sua ANOVA?
R: como dito anteriormente, o erro referente à Parcela (Erro a), visto que no
resultado da anova devo ter dois Erros (a e b).
- Outra perguntinha: você entende que ao especificar o termo erro você reduziu
o número de pontos para ajuste para dez e esse é o número de fatores PARC vezes
o número de parcelas SubP?
R: não consigo compreender isto.
Se você, ou alguém aqui do fórum, puder me ajudar, agradeço muito, estou
ansiosa para resolver este problema. Se eu estiver errada em algum ponto, por
favor, corrijam-me. Se o problema está no modelo, qual seria o que eu poderia
utilizar para parcela subdividida (PS) (tenho um experimento em DIC, com
parcelas subdividas (2 x 5) e 4 repetições)?
Em suma, meu objetivo é gerar o resíduo (que não sei se utilizei o comando
correto, visto que é uma PS) para poder realizar o teste de Shapiro-Wilk (SW),
pois o quadro da anova que utilizo é aquele do pacote 'ExpDes.pt' (e,
infelizmente, para o caso de PS este pacote não apresenta o resultado do teste
de SW). Assim, se houver uma solução ainda mais simples para que possa atingir
este objetivo, peço, por gentileza, que me apresentem.
Muito obrigada.
Aguardo retorno.
Att.,
Angélica.
________________________________
De: Cesar Rabak <[email protected]>
Enviado: sexta-feira, 9 de setembro de 2016 16:14
Para: Angélica Ricarte
Assunto: Re: [R-br] Parcela Subdividida
Olá Angélica,
CMR (que está descrito no FAQ deste grupo) é "código mínimo reproduzível".
A ideia é permitir que a gente refaça seus passos e acelere a ajuda.
Seu código acima não é exatamente reproduzível porque esta linha:
> PS<-read.xlsx(file.choose(),sheetName="ParcSub_DIC",h=T)
diz respeito a um arquivo que está na sua máquina. Contudo, como você imprimiu
o total do dataframe importado co o comando seguinte, está quase!
OK, então você tem uma resposta "var" para o fator PARC e aninhado em PARC você
tem quatro subparcelas, com quatro repetições por medida.
Por que então que você constrói (especifica) o seu termo Error como PARC:REP?
Ou fazendo uma pergunta mais próxima da teoria da análise que você quer fazer
(uma ANOVA com medidas repetidas), qual é o erro "dentro" que você quer levar
em conta na sua ANOVA?
Outra perguntinha: você entende que ao especificar o termo erro você reduziu o
número de pontos para ajuste para dez e esse é o número de fatores PARC vezes o
número de parcelas SubP?
HTH
--
Cesar Rabak
2016-09-08 14:27 GMT-03:00 Angélica Ricarte
<[email protected]<mailto:[email protected]>>:
Boa tarde, Cezar Rabak.
Entendi melhor, porém ainda não por completo.
Mas, o meu problema agora é que não consigo mais obter o mesmo resultado. Estou
importando exatamente os mesmos dados e utilizando exatamente os mesmos
comandos, porém não consigo mais o mesmo resultado que enviei no fórum dia
04/09/16. Quando realizo a ANOVA o resultado gerado está sendo diferente (e
está errado) daquele que apresentei anteriormente (este está correto). Gostaria
de saber como posso obter o "resíduo" para inserir no comando para o teste de
Shapiro Willk: "shapiro.test(resíduo)".
Seguem abaixo os comandos utilizados (acho que é o 'CMR' a que você se
referiu) e o dados, alguém poderia me ajudar a encontrar o meu erro e mostrar
uma solução? Já procurei muito, mas não encontro.
_____________________________________________________________
> PS<-read.xlsx(file.choose(),sheetName="ParcSub_DIC",h=T);PS
PARC SubP REP var
1 C 30 1 12
2 C 30 2 13
3 C 30 3 14
4 C 30 4 12
5 C 45 1 18
6 C 45 2 17
7 C 45 3 16
8 C 45 4 17
9 C 60 1 22
10 C 60 2 21
11 C 60 3 24
12 C 60 4 23
13 C 75 1 24
14 C 75 2 23
15 C 75 3 22
16 C 75 4 24
17 C 90 1 21
18 C 90 2 20
19 C 90 3 21
20 C 90 4 19
21 S 30 1 22
22 S 30 2 21
23 S 30 3 24
24 S 30 4 23
25 S 45 1 21
26 S 45 2 22
27 S 45 3 23
28 S 45 4 21
29 S 60 1 20
30 S 60 2 18
31 S 60 3 19
32 S 60 4 20
33 S 75 1 17
34 S 75 2 15
35 S 75 3 16
36 S 75 4 15
37 S 90 1 12
38 S 90 2 13
39 S 90 3 12
40 S 90 4 11
> attach(PS)
> SubP<-as.factor(SubP);is.factor(SubP)
[1] TRUE
> mod<-aov(var~PARC*SubP+Error(PARC:REP));summary(mod)
Error: PARC:REP
Df Sum Sq Mean Sq F value Pr(>F)
PARC 1 6.453 6.453 80.67 0.0706 .
Residuals 1 0.080 0.080
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
PARC 1 1.7 1.67 1.535 0.226
SubP 4 108.9 27.23 25.076 6.64e-09 ***
PARC:SubP 4 496.9 124.22 114.418 < 2e-16 ***
Residuals 28 30.4 1.09
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> resíduo<-residuals(mod);resíduo
NULL
___________________________________________________________________
Aguardo retorno.
Obrigada.
Att.,
Angélica.
________________________________
From: [email protected]<mailto:[email protected]>
Date: Tue, 6 Sep 2016 21:49:25 -0300
Subject: Re: [R-br] Parcela Subdividida
To: [email protected]<mailto:[email protected]>;
[email protected]<mailto:[email protected]>
Angélica,
Como você deu nomes "crípticos" para seu dataframe, e não postou um CMR, dá só
para fazer uma 'interpretação' das mensagens.
Aparentemente você quer a resposta PS[,4] para a interação entre as variáveis
(grupos[?]) PS[,1] e PS[,2] e considera a combinação dos fatores PS[,1] e
PS[,3] como medidas repetidas e tem 28 casos ao todo.
Além do mais você considera o fator PS[,1] como resposta e Error() ou seja
parte da especificação das medidas repetidas.
A mensagem sobre o modelo Error() ser singular ocorre quando você tenta ajustar
mais parâmetros que seus dados permitem...
Se você continuar com dificuldades, você precisa explicar melhor o que deseja
fazer ou consultar uma obra sobre ANOVA com medidas repetidas.
HTH
--
Cesar Rabak
2016-09-05 21:05 GMT-03:00 Angélica Ricarte via R-br
<[email protected]<mailto:[email protected]>>:
Prezado Cesar,
Não consegui compreender a mensagem de erro e nem o porquê que minha matriz é
nula. Você pode ajudar?
Agradeço e aguardo retorno.
Obrigada.
Att.,
Angélica
________________________________
Date: Mon, 5 Sep 2016 12:52:41 -0700
From:
[email protected]<mailto:ml-node%[email protected]>
To: [email protected]<mailto:[email protected]>
Subject: Re: [R-br] Parcela Subdividida
Angélica,
Você entendeu a mensagem de erro?
Não é possível analisar uma matriz nula....
2016-09-04 1:22 GMT-03:00 Angélica Ricarte via R-br <[hidden email]>:
Prezados,
Estou tentando realizar uma análise de variância para um experimentos em
parcelas subdivididas em DIC, mas, apesar de ser gerado um resultado, aparece
uma mensagem de erro. Quando tento realizar o cálculo dos resíduos também dá
erro e, consequentemente, não consigo fazer o teste de normalidade. Alguém pode
me ajudar? Vejam como realizei as análises e os seus resultados:
#ANOVA:
> mod<-aov(PS[,4]~PS[,1]*PS[,2]+Error(PS[,1]:PS[,3]));summary(mod)
Warning message:
In aov(PS[, 4] ~ PS[, 1] * PS[, 2] + Error(PS[, 1]:PS[, 3])) :
modelo Error() é singular
Error: PS[, 1]:PS[, 3]
Df Sum Sq Mean Sq F value Pr(>F)
PS[, 1] 1 8.1 8.100 11.3 0.0152 *
Residuals 6 4.3 0.717
---
Signif. codes:
0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
PS[, 2] 4 108.9 27.23 24.94 3.02e-08 ***
PS[, 1]:PS[, 2] 4 496.9 124.23 113.79 3.09e-15 ***
Residuals 24 26.2 1.09
---
Signif. codes:
0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
#Resíduos:
> residuals(mod)
NULL
#Teste de Shapiro-Wilk:
> res<-residuals(mod)
> shapiro.test(res)
Error: is.numeric(x) is not TRUE
Obrigada.
Att.,
Angélica.
_______________________________________________
R-br mailing list
[hidden email]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
_______________________________________________
R-br mailing list
[hidden email]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne?a c?digo
m?nimo reproduz?vel.
________________________________
If you reply to this email, your message will be added to the discussion below:
http://r-br.2285057.n4.nabble.com/R-br-Parcela-Subdividida-tp4666613p4666617.html
To unsubscribe from R-br, click
here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=YW5nZWxpY2EucmljYXJ0ZUBob3RtYWlsLmNvbXwzMzU3OTgyfDEwODg4MzAyMzk=>.
NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble:email.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble:email.naml-instant_emails%21nabble:email.naml-send_instant_email%21nabble:email.naml>
_______________________________________________
R-br mailing list
[email protected]<mailto:[email protected]>
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo
m�nimo reproduz�vel.
