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