[R-br] Dúvidas - Scott-Knott

Paulo Dick paulopcdick em gmail.com
Terça Junho 28 10:50:15 BRT 2016


Isso que a Daniela falou.

A função aqui "funcionou" assim: sk <- SK(x = as.factor(A$M), y = A$R,
model = 'y ~ x', which = 'x')

Digo "funcionou" porque desconheço esse método, não sei se a saída faz
sentido.


*Paulo Dick*
Estatístico
Mestre em Epidemiologia em Saúde Pública
Tel.: (55 21) 99591-2716

Em 28 de junho de 2016 10:48, Daniela Petruzalek <
daniela.petruzalek em gmail.com> escreveu:

> O modelo está reclamando que não existe uma variável categorica no seu
> data.frame. pelo o que me parece esta seria a M, correto?
>
> neste caso você pode converter A$M usando a rotina factor ou as.factor:
>
> A$M <- as.factor(A$M)
>
> Att.,
>
> Dani
>
> Em ter, 28 de jun de 2016 às 10:44, Pedro Campelo via R-br <
> r-br em listas.c3sl.ufpr.br> escreveu:
>
>> Paulo, segue o código.
>>
>> Os dados estão no data.frame anterior
>>
>>
>> #* Abrir pacote estatístico*
>> *require(agricolae)*
>> *require(ScottKnott)*
>> *require("ggplot2")*
>> *require("grid")*
>> *require("plyr")*
>> *require(reshape)*
>> *#Ler Dados*
>> *setwd("C:/Users/Pedro Campelo/Desktop")*
>> *A <- read.csv("C:/Users/Pedro Campelo/Desktop/A.csv", sep=";")*
>>
>> *data.frame(A)*
>>
>> *#ANOVA One Way*
>> *av<-aov(R ~ M, data.frame(A))*
>> *summary(av)*
>>
>> *#Plotar boxplot*
>> *par(font.axis=2)*
>> *boxplot(R~M,data.frame(A), xlab=expression(bold("Material")),
>> ylab=expression(bold("Value(%)")), col="lightgray")*
>>
>> *#Teste de Duncan*
>> *comparison <- duncan.test(av, 'M')*
>> *duncan.test(av, 'M', alpha = 0.05, console = TRUE)*
>>
>> *#Teste de Scott Knott*
>> *sk <- SK(x=av, y=A$R,  model='y~x', which='x', dispersion='se')*
>>
>> *summary(sk)*
>>
>>
>> *----*
>>
>> *Pedro Henrique Campelo Felix*
>> *Bachelor of Food Engineering (UNI-BH)*
>> *Msc. Chemical Engineering (UFMG)*
>> *PhD Student of Food Science (UFLA)*
>>
>> *+55 (31) 97560-1545 <%2B55%20%2831%29%2097560-1545> (TIM)*
>> *"**Nós somos aquilo que fazemos repetidamente. Excelência, então, não é
>> um modo de agir, mas um hábito." Aristóteles*
>>
>> Em 28 de junho de 2016 10:40, Paulo Dick <paulopcdick em gmail.com>
>> escreveu:
>>
>>> Pedro,
>>> Forneça um código mais reproduzível para podermos auxiliar, incluindo
>>> libraries
>>>
>>>
>>> *Paulo Dick*
>>> Estatístico
>>> Mestre em Epidemiologia em Saúde Pública
>>> Tel.: (55 21) 99591-2716
>>>
>>> 2016-06-28 10:17 GMT-03:00 Pedro Campelo via R-br <
>>> r-br em listas.c3sl.ufpr.br>:
>>>
>>>> Pessoal,
>>>>
>>>> Preciso de uma ajuda na execução desse script abaixo:
>>>>
>>>>
>>>> *#Data.Frame*
>>>>
>>>> *> data.frame(A)*
>>>> *   M     R*
>>>> *1  1 0.150*
>>>> *2  1 0.117*
>>>> *3  1 0.125*
>>>> *4  2 0.312*
>>>> *5  2 0.288*
>>>> *6  2 0.262*
>>>> *7  3 0.318*
>>>> *8  3 0.315*
>>>> *9  3 0.299*
>>>> *10 4 0.395*
>>>> *11 4 0.312*
>>>> *12 4 0.365*
>>>>
>>>> *#ANOVA One Way*
>>>> *av<-aov(R ~ M, data.frame(A))*
>>>> *summary(av)*
>>>>
>>>> *#Teste de Scott Knott*
>>>> *sk <- SK(x=av, y=A$R,  model='y~x', which='x', dispersion='se')*
>>>>
>>>> *summary(sk)*
>>>>
>>>> Quando vou executar a rotina, aparece esta mensagem:
>>>>
>>>> *> sk <- SK(x=av, y=A$R,  model='y~x', which='x', dispersion='se')*
>>>> *Error in SK.aov(x = av, y = A$R, model = "y~x", which = "x",
>>>> dispersion = "se") : *
>>>> *  No factors in the fitted model!*
>>>> *In addition: Warning message:*
>>>> *In replications(paste("~", xx), data = mf) : non-factors ignored: M*
>>>>
>>>>
>>>> Como poderia resolver isto?
>>>>
>>>> Obrigado
>>>>
>>>> *----*
>>>>
>>>> *Pedro Henrique Campelo Felix*
>>>> *Bachelor of Food Engineering (UNI-BH)*
>>>> *Msc. Chemical Engineering (UFMG)*
>>>> *PhD Student of Food Science (UFLA)*
>>>>
>>>> *+55 (31) 97560-1545 <%2B55%20%2831%29%2097560-1545> (TIM)*
>>>> *"**Nós somos aquilo que fazemos repetidamente. Excelência, então, não
>>>> é um modo de agir, mas um hábito." Aristóteles*
>>>>
>>>> _______________________________________________
>>>> 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.
>>>>
>>>
>>>
>> _______________________________________________
>> 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/20160628/0d787947/attachment.html>


Mais detalhes sobre a lista de discussão R-br