#Here i have plotted bill length distribution as per their species and island library(ggplot2) head(penguins) ggplot(penguins,aes(species,bill_length_mm, fill= island))+ geom_boxplot() # in your case you will have following code ggplot(data,aes(variable,value, fill= tretament))+ geom_boxplot()

@@DevResearch Thank you. But I think you did not understand my question. For example, I have three treatments, each with 8 replications. For each treatment, I have 9 response variables (it can be more or less) like body weight, etc... But it is very difficult to write all the variable values for all the treatments and replications. Is there any short method or I will be imported the excel data into Rstudio. Does it recognize the software by itself for testing the normaality test? Once again thanks.