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

Pedro Campelo pcampelo.felix em gmail.com
Terça Junho 28 10:44:26 BRT 2016


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.
>>
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160628/fa22715b/attachment.html>


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