Good morning
Please, is it a package named 'biology' in R?
I did the anova test and would like to use a bargraph to appreciate the difference between factors. the post-hoc Tukey"s test is associated . something like this:

So, I applied the instructions (script as indicated in the handbook) with my data as you can see with the script bellow but it does't work:
library(biology)
Mbargraph(dataframe1$NestGroupSize, dataframe1$VegetationType2, symbols = c("A", "AB", +"B", "AB"), ylab = "Mean number of nests", xlab = "Vetgetation type")

but when I run the script this is the results in red:

library(biology)
Error in library(biology) : aucun package nommé ‘biology’ n'est trouvé
Mbargraph(dataframe1$NestGroupSize, dataframe1$VegetationType2, symbols = c("A", "AB", "B", "AB"), ylab = "Mean number of nests", xlab = "Vetgetation type")
Error in Mbargraph(dataframe1$NestGroupSize, dataframe1$VegetationType2, :
impossible de trouver la fonction "Mbargraph"
if ther is anothe roption please let me know

It seems like biology package is not available at CRAN or any mainstream package source, look in your documentation of you can find a reference to the source or the author.

I took the data available for the book and made a graph with the following code. Most of the code is just wrangling the data to facilitate the layout in your example.

[quote="FJCC, post:4, topic:62035"]
Thanks you FJCC for the suggestion. However, whenI apply it to my own that it does't work. I instead have the following result:

> FIT <- lm(NestGroupSize ~ VegetationType2, data = dataframe1)
> library(broom)
> TidyFit <- tidy(FIT)
> TidyFit
# A tibble: 4 x 5
term estimate std.error statistic p.value
<chr> <dbl> <dbl> <dbl> <dbl>
1 (Intercept) 3.81 0.267 14.3 1.35e-30
2 VegetationType2Open forest 1.14 0.447 2.55 1.16e- 2
3 VegetationType2Wooding sav~ 0.332 0.461 0.720 4.73e- 1
4 VegetationType2Secondary f~ -0.815 2.42 -0.337 7.37e- 1
> TidyFit$term <- c("Mature forest", " Wooding savannah", " Open forest", "Secondary forest", "Grass savannah", "Swampy Zone"," Semi-Swampy zone")
Error in `$<-.data.frame`(`*tmp*`, term, value = c("Mature forest", " Wooding savannah", :
replacement has 7 rows, data has 4
> TidyFit$term <- factor(TidyFit$term, levels = c("Mature forest", "Wooding savannah", " Open forest", "Secondary forest", "Grass savannah", "Swampy Zone", "Semi-Swampy zone"), ordered = TRUE)
> Offset <- TidyFit[1, "estimate", drop = TRUE]
> TidyFit[4:7, "estimate"] <- TidyFit[4:7, "estimate"] + Offset
> TidyFit$Symbols <- c( "A", "B", "C","AB","BA","CB", "AC")
> TidyFit
# A tibble: 7 x 6
term estimate std.error statistic p.value Symbols
* <ord> <dbl> <dbl> <dbl> <dbl> <chr>
1 NA 3.81 0.267 14.3 1.35e-30 A
2 NA 1.14 0.447 2.55 1.16e- 2 B
3 NA 0.332 0.461 0.720 4.73e- 1 C
4 NA 3.00 2.42 -0.337 7.37e- 1 AB
5 NA NA NA NA NA BA
6 NA NA NA NA NA CB
7 NA NA NA NA NA AC

I would like to remind that my goal cpnsist to compare factors means from anova with a bargraph associated to leters. Factors in the graph previously attached with same letters indicate that there is not signifacnt difference among them.
Thanks
Thanks

The package is available from the authors homepage: http://users.monash.edu.au/~murray/BDAR/
But be warned! This package is 10+ Years old and was written for a very out of date version of R.
So do not expect it to work flawlessly with a modern setup!

Thanks AssignUser for the biology package. I got it but as you said it is not compatible with my R version (3.6).
Concerning the reprex I don't understand how it works. however this is my worry.
I have survey a forest to count ape nests and these nests fell in seven vegetation types (that I considered as factors for this variable). ANOVA help me to determine significant differences between the vegetation types and Tukey test have helped me to compare the respective difference of each factors to another. I now want to show this differences in a graphic like the one from Murray book above . Unfortunately, the procedure that Murray used did not help me. Murray suggested in his book to use library(biology) package and the script Mbargraph(medley$DIVERSITY, medley$ZINC, symbols = c("A", "AB","B", "AB"), ylab = "Mean diatom diversity",xlab = "Zinc concentration")
is it any option that can help me to meet my goal? looking forward to hear from you.
Thanks

I have no idea what you are trying to do with this. You have a 4 row dataframe as a result of fitting your model, and you then try adding a 7row column to it ? What and why...

Hello Dear nirgrahamuk
I have survey a forest to count ape nests and these nests fell in seven vegetation types (that I considered as factors for this variable). ANOVA help me to determine significant differences between the vegetation types and Tukey test have helped me to compare the respective difference of each factors to another. I am trying to show this differences in a graphic as Logan did with the example about Zinc concentration with diatom diversity

Don't they show that you only have 3 values of Vegetation Type 2 that would have estimates different from the base case (Intercept) - so at most you can plot charts of 4 different bars, and not 7 .
(well I suppose you could plot 7 but there would be 4 identical bars only differing by the VegetationType2 label, and not by prediction.