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

[Daniela Petruzalek]

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 (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/12204aa5/attachment.html>