Prezados,

Ao rodar a programação notei ainda que os g.l do bloco estão errados (tenho 3 blocos, logo g.l=2)

rm(list=ls())

dados <- read.table("G:/dados.txt", head=T)

head(dados)

__________________________

produtividade hibrido populacao bloco zona

1 101 30A37 60 1 1

2 131 30A37 60 1 1

3 138 30A37 60 1 1

4 128 30A37 60 1 1

5 101 30A37 60 1 1

6 89 30A37 60 1 1

____________________________

library(ExpDes)

fat3.rbd(dados$hibrido, dados$populacao, dados$zona, dados$bloco,

dados$produtividade, sigT=0.1, sigF=0.1)

____________________________

------------------------------------------------------------------------

Legend:

FACTOR 1: F1

FACTOR 2: F2

FACTOR 3: F3

------------------------------------------------------------------------

Analysis of Variance Table

------------------------------------------------------------------------

DF SS MS Fc Pr>Fc

Block 259.4167 7279.352 28.06046 0.1136 1

F1 3.0000 1149513.560 383171.18652 1551.6394 0

F2 4.0000 3798747.062 949686.76545 3845.7259 0

F3 2.0000 714265.142 357132.57111 1446.1968 0

F1*F2 12.0000 202406.001 16867.16674 68.303 0

F1*F3 6.0000 13306.876 2217.81275 8.981 0

F2*F3 8.0000 23012.734 2876.59176 11.6487 0

F1*F2*F3 24.0000 62022.721 2584.28002 10.465 0

Residuals 15563.0000 3843221.113 246.94603

Total 15622.0000 9813774.560 628.20219

------------------------------------------------------------------------

Erro em shapiro.test(anava$residuals) :

sample size must be between 3 and 5000

Em 14 de agosto de 2013 14:37, Tiago Souza Marçal <tiagosouzamarcal@hotmail.com> escreveu:

Natalia ente nos códigos da função e edite esta parte retirando o teste de Shapiro wilk.

Depois e só colar a função no R e utilizar a função que você tinha mostrado.

fat3.rbd(hibrido, populacao, zona, bloco, Prod, sigT = 0.1, sigF = 0.1) # A propósito o alpha que você esta utilizando é de 10%.

Att.

Tiago.

#################################################################

Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES

Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas

#################################################################

Date: Wed, 14 Aug 2013 14:03:08 -0300
From: nsmbarreto@gmail.com
To: R-br@listas.c3sl.ufpr.br
Subject: [R-br] ExpDes

Boa tarde prezados,
estou fazendo uma anova com 3 fatores e bloco

fat3.rbd(hibrido, populacao, zona, bloco, Prod, sigT = 0.1, sigF = 0.1)

No entanto, a analise nao 'e finalizada pois o teste de shapiro nao 'e aplicado para amostras maiores que 5000.

Por favor o que devo fazer para utilizar esse pacote?

Att

--

Natália da Silva Martins
Contato: (19) 8306-4743

_______________________________________________ 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.

--

Natália da Silva Martins
Bacharel em Estatística - Universidade Estadual de Maringá/ UEM
Mestra em Ciências: Area de concentração: Estatística e Experimentação Agronômica - ESALQ/ USP
Doutoranda em Ciências: Area de concentração: Estatística e Experimentação Agronômica - ESALQ/ USP
Contato: (19) 8306-4743