Re: [R-br] Parcela Subdividida

2016-09-08 Por tôpico Angélica Ricarte via R-br 
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
10C   60   2  21
11C   60   3  24
12C   60   4  23
13C   75   1  24
14C   75   2  23
15C   75   3  22
16C   75   4  24
17C   90   1  21
18C   90   2  20
19C   90   3  21
20C   90   4  19
21S   30   1  22
22S   30   2  21
23S   30   3  24
24S   30   4  23
25S   45   1  21
26S   45   2  22
27S   45   3  23
28S   45   4  21
29S   60   1  20
30S   60   2  18
31S   60   3  19
32S   60   4  20
33S   75   1  17
34S   75   2  15
35S   75   3  16
36S   75   4  15
37S   90   1  12
38S   90   2  13
39S   90   3  12
40S   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   11.71.67   1.5350.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.41.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: cesar.ra...@gmail.com
Date: Tue, 6 Sep 2016 21:49:25 -0300
Subject: Re: [R-br] Parcela Subdividida
To: angelica.rica...@hotmail.com; r-br@listas.c3sl.ufpr.br

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 :



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: ml-node+s2285057n417...@n4.nabble.com
To: angelica.rica...@hotmail.com
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]18.1   8.10011.3 0.0152 *
Residuals  64.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

Re: [R-br] Parâmetros do plot

2016-09-08 Por tôpico Andre Oliveira via R-br 
Olhei sim Cesar,  mas não obtive exito para tirar os vetores e deixar só os 
tratamentos.  André Oliveira Souza. Graduação em Matemática, mestrado em 
estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do 
Espirito Santo.  IFES 

Em Quarta-feira, 7 de Setembro de 2016 15:27, Cesar Rabak 
 escreveu:
 

 Você viu a documentação do plot.bpca2d no help do pacote?


2016-09-07 8:34 GMT-03:00 Andre Oliveira via R-br :

Bom dia,pesquisei e conversei com algumas pessoas e não resolvi o meu problema. 
Como tirar estes vetores desta imagem feita com o plot?require(bpca)
fit=bpcaplot(bpca(iris[2:4],center =2),var.factor = 0.2)
obrigado
 André Oliveira Souza. Graduação em Matemática, mestrado em estatística 
aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. 
 IFES
__ _
R-br mailing list
R-br@listas.c3sl.ufpr.br
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
R-br@listas.c3sl.ufpr.br
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.