[R-br] ScottKnott

[Tiago Souza Marçal]

Caso você queira conferir tente usar a função do pacote EXPDES.pt::dic(trata,respo,quali=c(TRUE),mcomp="sk").
 
Tente executar seu comando da seguinte forma:
 
 sk1<-with(dad,SK(y=respo, model="respo~trata",which="trata"))
 
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: Mon, 2 Dec 2013 17:39:25 -0800
From: andreolsouza em yahoo.com.br
To: r-br em listas.c3sl.ufpr.br
Subject: [R-br] ScottKnott

Pessoal depois de pesquisar na grupo e no manual da função ainda não resolvi o que quero. Estou com esta estrutura. 


> setwd("/media/Meus arquivos/Pendrive/KINGSTON/sweave/apostilaR2")> require(gdata)> dad<-read.xls("dic.xls")> str(dad)'data.frame':	12 obs. of  3 variables: $ trata: Factor w/ 4 levels "A","B","C","D": 1 1 1 2 2 2 3 3 3 4 ... $ rep  : int  1 2 3 1 2 3 1 2 3 1 ... $ respo: int  6 14 4 9 10 8 73 22 10 1 ...> attach(dad)> dad   trata rep respo1      A   1     62      A   2    143      A   3     44      B   1     95      B   2    106      B   3     87      C   1    738      C   2    229      C   3    1010     D   1     111     D   2    1012     D   3     4> library(doBy)> summaryBy(respo~trata,dad,FUN=c(mean,sd))  trata respo.mean  respo.sd1     A          8  5.2915032     B          9  1.0000003     C         35 33.4514574     D          5  4.582576> require(ScottKnott)> sk1<-SK(dad,y=respo,model="respo~trata",which="trata")> summary(sk1) Levels Means SK(5%)      C    35      a      B     9      a      A     8      a      D     5      a> plot(sk1, title = "Treatamentos", col = rainbow(3))
média iguais ? 
 Minha pergunta é o que fiz de errado? Declarei o modelo errado?  o correto seria sk1<-SK(x=dad,y=respo,model="respo~trata",which="trata")? Não mudou nada na resposta summary(sk1) Levels Means SK(5%)      C    35      a      B     9      a      A     8      a      D     5      a
obrigado


_______________________________________________
R-br mailing list
R-br em 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. 		 	   		  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131203/1a5c977a/attachment.html>