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@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@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@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@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@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.