[R-br] ScottKnott

[Andre Oliveira]

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 respo
1      A   1     6
2      A   2    14
3      A   3     4
4      B   1     9
5      B   2    10
6      B   3     8
7      C   1    73
8      C   2    22
9      C   3    10
10     D   1     1
11     D   2    10
12     D   3     4
> library(doBy)
> summaryBy(respo~trata,dad,FUN=c(mean,sd))
  trata respo.mean  respo.sd
1     A          8  5.291503
2     B          9  1.000000
3     C         35 33.451457
4     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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131202/920784d5/attachment.html>