[R] graph in R with grouping letters from the turkey test with agricolae package
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 TREAT REP VAR1VAR2VAR3VAR4 t1 1 16,10,805 0,966 1,449 t1 2 14,10,705 0,846 1,269 t1 3 5,8 0,290,348 0,522 t1 4 4,9 0,245 0,294 0,441 t2 1 30,81,541,848 2,772 t2 2 29 1,451,742,61 t2 3 28,81,441,728 2,592 t2 4 16,60,830,996 1,494 t3 1 31,31,565 1,878 2,817 t3 2 29,61,481,776 2,664 t3 3 28,11,405 1,686 2,529 t3 4 19,50,975 1,171,755 t4 1 21,45 1,0725 1,287 1,9305 t4 2 23,595 1,17975 1,4157 2,12355 t4 3 25,9545 1,2977251,55727 2,335905 t4 4 28,549951,4274975 1,7129972,5694955 t5 1 27,31,365 1,638 2,457 t5 2 13,70,685 0,822 1,233 t5 3 8,1 0,405 0,486 0,729 t5 4 0 0 0 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.
[R] agricolae package, graph in R with grouping letters from the Tukey, LSD, Duncan test
I would like to produce the attached graph (graph1) with the R package agricolae, could someone give me an example with the attached data (dataDoubleFactor)? I expect an adapted graph (graphDoubleFactor) with the data (data) Best regards TREAT CUT VAR1 t1 1 16.10 t1 2 14.10 t1 3 5.80 t1 4 4.90 t2 1 30.80 t2 2 29.00 t2 3 28.80 t2 4 16.60 t3 1 31.30 t3 2 29.60 t3 3 28.10 t3 4 19.50 t4 1 21.45 t4 2 23.60 t4 3 25.95 t4 4 28.55 t5 1 27.30 t5 2 13.70 t5 3 8.10 t5 4 0.00 __ 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.
Re: [R] graph in R with grouping letters from the turkey test with agricolae package
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
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.
[R] Request to Unsubscribe from R Forum Replies
Hello, I would like to request the removal of my email address from the R forum reply/update list. I no longer wish to receive notifications or messages related to forum discussions. Thank you for your attention and understanding. Best regards, [[alternative HTML version deleted]] __ 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 https://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
