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

[Pedro Campelo]

Daniela e Paulo,

Deu certo. 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 (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:50, Paulo Dick <paulopcdick em gmail.com> escreveu:

> 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/22f4dc43/attachment.html>