Yes, the data and the R code used are attached.
I would like to produce the attached graph (graph1) with the R package
agricolae, could someone give me an example with the attached data
(vermiwash and Rcode_vermiwash)?
Fig. 7, https://doi.org/10.1007/s42729-023-01295-3
I expect an adapted graph (graphDoubleFactor) with the data (vermiwash and
Rcode_vermiwash)
Best regards
Em ter., 12 de set. de 2023 às 18:54, Rui Barradas <ruipbarra...@sapo.pt>
escreveu:
> Às 16:24 de 12/09/2023, Loop Vinyl escreveu:
> > I would like to produce the attached graph (graph1) with the R package
> > agricolae, could someone give me an example with the attached data
> (data)?
> >
> > I expect an adapted graph (graph2) with the data (data)
> >
> > Best regards
> >
> >
> > ______________________________________________
> > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> > and provide commented, minimal, self-contained, reproducible code.
> Hello,
>
> There are no attached graphs, only data.
> Can you post the code have you tried?
>
> Hope this helps,
>
> Rui Barradas
>
>
library(agricolae)
dados <- read.table("vermiwash.txt", header = TRUE)
head(dados)
View(dados)
DOSE <- factor(dados$TREAT)
levels(DOSE)
CUT <- factor(dados$CUT)
levels(CUT)
#com a variável 1
model_var1<-aov(VAR1~DOSE*CUT, data=dados)
summary(model_var1)
cv.model(model_var1)
#
#var1
tukey_var1 <- HSD.test(model_var1, c("DOSE", "CUT"), main="VAR1 ~ DOSE*CUT",
console=TRUE)
head(tukey_var1)
#
x11()
bar.group(tukey_var1$groups,main="Teste de Tukey", ylim=c(0,50),
xlab="Tratamentos (Dosagem)")
#OR
#var1
compara.tukey.agri <- HSD.test(model_var1, "DOSE" ,group=T)
compara.tukey.agri
x11()
bar.group(compara.tukey.agri$groups,main="Teste de Tukey", ylim=c(0,50),
xlab="Tratamentos (DOSE)")
TREAT CUT REP VAR1
Z 1 1 16.1
Z 1 2 14.1
Z 1 3 5.8
Z 1 4 4.9
Z 2 1 16.1
Z 2 2 14.1
Z 2 3 5.8
Z 2 4 4.9
Z 3 1 16.1
Z 3 2 14.1
Z 3 3 5.8
Z 3 4 4.9
Z 4 1 16.1
Z 4 2 14.1
Z 4 3 5.8
Z 4 4 4.9
5 1 1 30.8
5 1 2 29
5 1 3 28.8
5 1 4 16.6
5 2 1 30.8
5 2 2 29
5 2 3 28.8
5 2 4 16.6
5 3 1 30.8
5 3 2 29
5 3 3 28.8
5 3 4 16.6
5 4 1 30.8
5 4 2 29
5 4 3 28.8
5 4 4 16.6
10 1 1 31.3
10 1 2 29.6
10 1 3 28.1
10 1 4 19.5
10 2 1 31.3
10 2 2 29.6
10 2 3 28.1
10 2 4 19.5
10 3 1 31.3
10 3 2 29.6
10 3 3 28.1
10 3 4 19.5
10 4 1 31.3
10 4 2 29.6
10 4 3 28.1
10 4 4 19.5
15 1 1 27.3
15 1 2 13.7
15 1 3 8.1
15 1 4 0
15 2 1 27.3
15 2 2 13.7
15 2 3 8.1
15 2 4 0
15 3 1 27.3
15 3 2 13.7
15 3 3 8.1
15 3 4 0
15 4 1 27.3
15 4 2 13.7
15 4 3 8.1
15 4 4 0
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
