Multiple boxplots

[Bayan]

How can I do multiple box plot that shows for each biofilter two RT's 22 and 70 values of teh TOC_RE.
R Biofilter TOC_RE
1 22 D 47.67797
2 70 D 57.04705
3 22 D 26.04057
4 70 D 24.91069
5 22 D 30.22899
6 70 D 27.60135
7 22 D -59.64255
8 70 D 32.07621
9 22 D 37.97165
10 70 D 49.66546
11 22 D 62.18069
12 22 CB 46.01593
13 70 CB 55.68262
14 22 CB 26.04057
15 70 CB 37.05162
16 22 CB 42.99045
17 70 CB 35.48087
18 22 CB 18.31099
19 70 CB 37.30960
20 22 CB 45.46277
21 70 CB 47.39043
22 22 CB 62.16136
23 22 W 14.07349
24 70 W 29.45998
25 22 W 16.85578
26 70 W 17.40317
27 22 W 32.02720
28 70 W 30.88080
29 22 W 21.16420
30 70 W 30.01503
31 22 W 46.56881
32 70 W 44.99192
33 22 W 62.77956
34 22 M 43.56143
35 70 M 53.66764
36 22 M 33.58706
37 70 M 35.95476
38 22 M 37.82223
39 70 M 35.27777
40 22 M 26.82170
41 70 M 34.51255
42 22 M 46.72077
43 70 M 50.49602
44 22 M 63.40935
45 22 LS 30.92150
46 70 LS 43.29109
47 22 LS 27.03161
48 70 LS 27.85037
49 22 LS 40.01183
50 70 LS 34.15952
51 22 LS 24.03189
52 70 LS 35.79513
53 22 LS 44.65384
54 70 LS 48.08845
55 22 LS 62.19986
56 22 B5 56.31297
57 70 B5 64.13582
58 22 B5 56.20525
59 70 B5 51.39499
60 22 B5 51.14557
61 70 B5 47.93902
62 22 B5 43.99362
63 70 B5 47.54421
64 22 B5 54.15737
65 70 B5 54.79393
66 22 B5 68.21056
67 22 B10 62.82134
68 70 B10 69.47876
69 22 B10 64.47490
70 70 B10 59.62507
71 22 B10 61.58454
72 70 B10 57.91203
73 22 B10 61.06062
74 70 B10 63.16541
75 22 B10 69.66219
76 70 B10 68.52537
77 22 B10 76.99970

[FJCC]

Is this the kind of plot you want?

library(ggplot2)
DF <- read.csv("~/R/Play/Dummy.csv",sep=" ")
ggplot(DF, aes(x = Biofilter, y = TOC_RE, fill = as.factor(R))) + 
  geom_boxplot() +  labs(fill = "R")

^{Created on 2022-07-20 by the reprex package (v2.0.1)}

[Bayan]

Thank you so much !!! yes it is

[Ven]

You can also flip the plot with coord_flip(), given your x-axis is qualitative, It may be easier to read on the y-axis.

@FJCC

ggplot(DF, aes(x = Biofilter, y = TOC_RE, fill = as.factor(R))) + 
  geom_boxplot() + 
  coord_flip() + 
  labs(fill = "R")