Sub-setting data from ggplot2

tidyverse
,
1

Considering this figure, please I have been trying to get the number of data points that is above the 1.5:1 line and also the number of data point that is below the 1:1 line, but I don't know how to go about that. I have attached the sample dataset (first 50 observation) and code below:
I will be very grateful if someone can help me out

Packages

library(tidyverse)

Wrangling

D <- DF[!(DF$Tipo_Parc == "Extension" & DF$Ano == 2019), ]
#D <- DF[!(DF$Ano == 2016), ]

D$Genero <- as.factor(D$Genero)
D$Region <- as.factor(D$Region)
D$RainReg <- as.factor(D$RainReg)
D$Pais <- as.factor(D$Pais)
D$Department <- as.factor(D$Department)
D$ID_Prod <- as.factor(D$ID_Prod)
D$SubDep <- as.factor(D$SubDep)

#Ln Response Ration
D$RR <- (D$ASA / D$Testigo)
limit <- mean(D$RR) + 5 * sd(D$RR)
D <- D[!(D$RR > limit), ] #remove outlier yield differences
D <- D[!(D$RR < (limit * -1)), ] #remove outlier yield differences

#Remove no precip data and no location info
D <- D[!is.na(D$PriPET), ]
D <- D[!is.na(D$SubDep), ]

#percentage change
D$Dif <- D$ASA - D$Testigo # Absolute yield dif
D$Per <- (D$Dif) / D$Testigo * 100

length(unique(D$ID_Prod))#1271 producers

#Rainwater productivity
D$RWP.ASA <- NA
D$RWP.ASA <-
ifelse((D$Temporada == "Primera"),
D$ASA / D$PriPrecip,
D$ASA / D$PosPrecip)
D$RWP.Testigo <-
ifelse((D$Temporada == "Primera"),
D$Testigo / D$PriPrecip,
D$Testigo / D$PosPrecip)

D <- D %>% filter(Cultivo=='Maiz') %>% filter(Temporada=='Primera')

Modelling

ggplot(D, aes(Testigo, ASA, color = RR, group = 2)) +
geom_point() +
geom_abline(intercept = 0,
slope = 1,
linetype = "dashed") +
geom_abline(intercept = 0,
slope = 1.5,
linetype = "dashed") +
xlab(expression(Conventional ~ yield ~ (kg ~ ha ^ -1))) +
ylab(expression(WSA ~ yield ~ (kg ~ ha ^ -1))) +
theme(panel.grid.major = element_blank(),
panel.grid.minor = element_blank()) +
annotate(
geom = "text",
x = 6500,
y = 6850,
label = "1:1 line",
color = "black",
size=4
) +
annotate(
geom = "text",
x = 5000,
y = 8000,
label = "1.5:1 line",
color = "black",
size=4
) +
xlim(0, 8000) +
ylim(0, 8000) +
scale_colour_stepsn(
colors = c("#DCE319FF", "#3CBB75FF", "#20A387FF", "yellow3"),name = 'Yield response ratio',
#breaks = c(2,4,6),
guide = guide_coloursteps(even.steps = FALSE,
show.limits = TRUE,direction = 'horizontal',)
) +
theme_bw() +
theme(
# legend.position = "none",
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
text = element_text(
family = 'serif',
face = 'bold',
colour = 'black',
size = 15
),legend.position = c(0.2,0.8),
legend.title.position = 'top'
)

Rplot
Rplot1108×457 16.3 KB

2

Here is the dataset

DF <- structure(
list(
ID_Prod = c(
"200202401",
"200202401",
"200202401",
"200202401",
"200202401",
"200202402",
"200202402",
"200202402",
"200202402",
"200202402",
"200202402",
"200202402",
"200202402",
"200202404",
"200202404",
"200202404",
"200202404",
"200202404",
"200202405",
"200202405",
"200202405",
"200202405",
"200202406",
"200202406",
"200202406",
"200202406",
"200202406",
"200202407",
"200202407",
"200202407",
"200202407",
"200202407",
"200202408",
"200202408",
"200202408",
"200202408",
"200202408",
"200202409",
"200202409",
"200202409",
"200202409",
"200202409",
"200202410",
"200202410",
"200202410",
"200202410",
"200202410",
"200202411",
"200202411",
"200202411"
),
Ano = c(
2017L,
2018L,
2018L,
2016L,
2017L,
2018L,
2020L,
2016L,
2018L,
2019L,
2020L,
2019L,
2017L,
2018L,
2017L,
2018L,
2017L,
2016L,
2018L,
2017L,
2018L,
2017L,
2017L,
2016L,
2018L,
2017L,
2018L,
2018L,
2016L,
2018L,
2017L,
2017L,
2017L,
2018L,
2017L,
2018L,
2016L,
2016L,
2018L,
2018L,
2017L,
2017L,
2018L,
2017L,
2017L,
2018L,
2016L,
2017L,
2016L,
2017L
),
Pais = c(
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala",
"Guatemala"
),
Cultivo = c(
"Frijol",
"Frijol",
"Maiz",
"Frijol",
"Maiz",
"Maiz",
"Maiz",
"Frijol",
"Frijol",
"Frijol",
"Frijol",
"Maiz",
"Maiz",
"Frijol",
"Maiz",
"Maiz",
"Frijol",
"Frijol",
"Frijol",
"Frijol",
"Maiz",
"Maiz",
"Frijol",
"Frijol",
"Maiz",
"Maiz",
"Frijol",
"Frijol",
"Frijol",
"Maiz",
"Maiz",
"Frijol",
"Frijol",
"Frijol",
"Maiz",
"Maiz",
"Frijol",
"Frijol",
"Frijol",
"Maiz",
"Maiz",
"Frijol",
"Maiz",
"Maiz",
"Frijol",
"Frijol",
"Frijol",
"Maiz",
"Frijol",
"Frijol"
),
Temporada = c(
"Postrera",
"Postrera",
"Primera",
"Primera",
"Primera",
"Primera",
"Primera",
"Primera",
"Postrera",
"Postrera",
"Postrera",
"Primera",
"Primera",
"Postrera",
"Primera",
"Primera",
"Postrera",
"Primera",
"Postrera",
"Postrera",
"Primera",
"Primera",
"Postrera",
"Primera",
"Primera",
"Primera",
"Postrera",
"Postrera",
"Primera",
"Primera",
"Primera",
"Postrera",
"Postrera",
"Postrera",
"Primera",
"Primera",
"Primera",
"Primera",
"Postrera",
"Primera",
"Primera",
"Postrera",
"Primera",
"Primera",
"Postrera",
"Postrera",
"Primera",
"Primera",
"Postrera",
"Postrera"
),
Sistema = c(
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB",
"GB"
),
ASA = c(
2328.8,
1511.97293568445,
3645.83333333333,
204.54,
6133.275813,
1964.77272727273,
6508.515625,
159.09,
1889.96616960556,
1960.22727272727,
3073.95833333333,
3484.84848484848,
5208.85377,
2551.45432896751,
4487.901235,
3030.30303030303,
1199.111111,
272.73,
2173.4610950464,
1399.42029,
1515.15151515152,
4982.36102,
1366.075949,
136.36,
1254.73484848485,
4845.92201,
2078.96278656612,
2303.39626920678,
159.09,
639.015151515152,
3367.842514,
1202.823529,
1032.727273,
1795.46786112529,
3114.574315,
2107.00757575758,
1022.73,
954.54,
1984.46447808584,
7149.43181818182,
5692.957651,
1666.076942,
1704.54545454545,
5632.666667,
1136,
2840.85539868836,
409.09,
5847.058824,
331.81,
2272
),
Testigo = c(
1451.555556,
1039.48139328306,
2580.49242424242,
54.54,
4730.687654,
591.856060606061,
6242.5,
68.18,
1322.9763187239,
1524.62121212121,
1921.22395833333,
4621.21212121212,
3907.710876,
1795.46786112529,
4492.893338,
3408.90151515152,
326.4367816,
113.64,
1511.97293568445,
1234.782609,
639.015151515152,
4804.278976,
868.7058824,
90.91,
875.94696969697,
2545.658263,
1795.46786112529,
1612.37738844475,
90.91,
1136.36363636364,
2292.940566,
1044.597701,
826.1818182,
1133.97970176334,
3668.899522,
1799.24242424242,
681.82,
681.82,
1133.97970176334,
4545.26515151515,
3648.905673,
917.82,
1041.66666666667,
4039.296188,
568,
2456.95602048723,
90.91,
4538.852847,
195.45,
1231.807229
),
JulPET = c(
1.2277770996094,-78.2643783569336,
-78.2643783569336,
-18.6209518432617,
1.2277770996094,-78.9643783569336,
44.9514373779297,
-19.8209518432617,
-78.9643783569336,-68.0502014160156,
44.9514373779297,
-68.0502014160156,
0.627777099609403,-78.2643783569336,
1.2277770996094,
-78.2643783569336,
1.2277770996094,-18.6209518432617,
-78.9643783569336,
0.627777099609403,
-78.9643783569336,
0.627777099609403,
1.2277770996094,
-18.6209518432617,
-78.2643783569336,
1.2277770996094,
-78.2643783569336,
-78.9643783569336,
-19.8209518432617,-78.9643783569336,
0.627777099609403,
0.627777099609403,
1.32777709960939,-78.3643783569336,
1.32777709960939,
-78.3643783569336,
-18.6209518432617,-19.8209518432617,
-78.9643783569336,
-78.9643783569336,
0.627777099609403,
0.627777099609403,
-78.2643783569336,
1.2277770996094,
1.2277770996094,-78.2643783569336,
-18.6209518432617,
0.627777099609403,
-19.8209518432617,
0.627777099609403
),
PosPET = c(
156.561800765991,
-45.8779052734375,-45.8779052734375,
5.43591842651377,
156.561800765991,
-49.5779052734375,
258.977874755859,
1.63591842651374,
-49.5779052734375,
153.501736831665,
258.977874755859,
153.501736831665,
152.861800765991,
-45.8779052734375,
156.561800765991,
-45.8779052734375,
156.561800765991,
5.43591842651377,-49.5779052734375,
152.861800765991,
-49.5779052734375,
152.861800765991,
156.561800765991,
5.43591842651377,
-45.8779052734375,
156.561800765991,-45.8779052734375,
-49.5779052734375,
1.63591842651374,
-49.5779052734375,
152.861800765991,
152.861800765991,
155.961800765991,
-47.4779052734374,
155.961800765991,
-47.4779052734374,
4.73591842651375,
1.63591842651374,-49.5779052734375,
-49.5779052734375,
152.861800765991,
152.861800765991,-45.8779052734375,
156.561800765991,
156.561800765991,
-45.8779052734375,
5.43591842651377,
152.861800765991,
1.63591842651374,
152.861800765991
),
PriPET = c(
327.861569213867,
4.3409301757813,
4.3409301757813,
76.2892517089844,
327.861569213867,
1.94093017578132,
233.491961669922,
73.8892517089844,
1.94093017578132,
142.833575439453,
233.491961669922,
142.833575439453,
326.161569213867,
4.3409301757813,
327.861569213867,
4.3409301757813,
327.861569213867,
76.2892517089844,
1.94093017578132,
326.161569213867,
1.94093017578132,
326.161569213867,
327.861569213867,
76.2892517089844,
4.3409301757813,
327.861569213867,
4.3409301757813,
1.94093017578132,
73.8892517089844,
1.94093017578132,
326.161569213867,
326.161569213867,
331.361569213867,
6.84093017578132,
331.361569213867,
6.84093017578132,
79.7892517089844,
73.8892517089844,
1.94093017578132,
1.94093017578132,
326.161569213867,
326.161569213867,
4.3409301757813,
327.861569213867,
327.861569213867,
4.3409301757813,
76.2892517089844,
326.161569213867,
73.8892517089844,
326.161569213867
),
CanPET = c(
70.3991729736328,-72.7154342651367,
-72.7154342651367,
33.1415420532227,
70.3991729736328,-74.6154342651367,
103.168469238281,
30.8415420532227,
-74.6154342651367,-11.0964630126953,
103.168469238281,
-11.0964630126953,
68.7991729736328,-72.7154342651367,
70.3991729736328,
-72.7154342651367,
70.3991729736328,
33.1415420532227,
-74.6154342651367,
68.7991729736328,
-74.6154342651367,
68.7991729736328,
70.3991729736328,
33.1415420532227,
-72.7154342651367,
70.3991729736328,
-72.7154342651367,
-74.6154342651367,
30.8415420532227,-74.6154342651367,
68.7991729736328,
68.7991729736328,
70.8991729736328,-72.6154342651367,
70.8991729736328,
-72.6154342651367,
33.4415420532227,
30.8415420532227,
-74.6154342651367,
-74.6154342651367,
68.7991729736328,
68.7991729736328,
-72.7154342651367,
70.3991729736328,
70.3991729736328,-72.7154342651367,
33.1415420532227,
68.7991729736328,
30.8415420532227,
68.7991729736328
),
JulPrecip = c(
140.127777099609,
68.0356216430664,
68.0356216430664,
123.079048156738,
140.127777099609,
68.0356216430664,
184.55143737793,
123.079048156738,
68.0356216430664,
85.4497985839844,
184.55143737793,
85.4497985839844,
140.127777099609,
68.0356216430664,
140.127777099609,
68.0356216430664,
140.127777099609,
123.079048156738,
68.0356216430664,
140.127777099609,
68.0356216430664,
140.127777099609,
140.127777099609,
123.079048156738,
68.0356216430664,
140.127777099609,
68.0356216430664,
68.0356216430664,
123.079048156738,
68.0356216430664,
140.127777099609,
140.127777099609,
140.127777099609,
68.0356216430664,
140.127777099609,
68.0356216430664,
123.079048156738,
123.079048156738,
68.0356216430664,
68.0356216430664,
140.127777099609,
140.127777099609,
68.0356216430664,
140.127777099609,
140.127777099609,
68.0356216430664,
123.079048156738,
140.127777099609,
123.079048156738,
140.127777099609
),
PriPrecip = c(
854.061569213867,
543.640930175781,
543.640930175781,
627.389251708984,
854.061569213867,
543.640930175781,
776.791961669922,
627.389251708984,
543.640930175781,
733.133575439453,
776.791961669922,
733.133575439453,
854.061569213867,
543.640930175781,
854.061569213867,
543.640930175781,
854.061569213867,
627.389251708984,
543.640930175781,
854.061569213867,
543.640930175781,
854.061569213867,
854.061569213867,
627.389251708984,
543.640930175781,
854.061569213867,
543.640930175781,
543.640930175781,
627.389251708984,
543.640930175781,
854.061569213867,
854.061569213867,
854.061569213867,
543.640930175781,
854.061569213867,
543.640930175781,
627.389251708984,
627.389251708984,
543.640930175781,
543.640930175781,
854.061569213867,
854.061569213867,
543.640930175781,
854.061569213867,
854.061569213867,
543.640930175781,
627.389251708984,
854.061569213867,
627.389251708984,
854.061569213867
),
PosPrecip = c(
610.761800765991,
415.522094726563,
415.522094726563,
485.335918426514,
610.761800765991,
415.522094726563,
721.477874755859,
485.335918426514,
415.522094726563,
663.801736831665,
721.477874755859,
663.801736831665,
610.761800765991,
415.522094726563,
610.761800765991,
415.522094726563,
610.761800765991,
485.335918426514,
415.522094726563,
610.761800765991,
415.522094726563,
610.761800765991,
610.761800765991,
485.335918426514,
415.522094726563,
610.761800765991,
415.522094726563,
415.522094726563,
485.335918426514,
415.522094726563,
610.761800765991,
610.761800765991,
610.761800765991,
415.522094726563,
610.761800765991,
415.522094726563,
485.335918426514,
485.335918426514,
415.522094726563,
415.522094726563,
610.761800765991,
610.761800765991,
415.522094726563,
610.761800765991,
610.761800765991,
415.522094726563,
485.335918426514,
610.761800765991,
485.335918426514,
610.761800765991
),
CanPrecip = c(
344.199172973633,
209.984565734863,
209.984565734863,
307.541542053223,
344.199172973633,
209.984565734863,
376.368469238281,
307.541542053223,
209.984565734863,
290.703536987305,
376.368469238281,
290.703536987305,
344.199172973633,
209.984565734863,
344.199172973633,
209.984565734863,
344.199172973633,
307.541542053223,
209.984565734863,
344.199172973633,
209.984565734863,
344.199172973633,
344.199172973633,
307.541542053223,
209.984565734863,
344.199172973633,
209.984565734863,
209.984565734863,
307.541542053223,
209.984565734863,
344.199172973633,
344.199172973633,
344.199172973633,
209.984565734863,
344.199172973633,
209.984565734863,
307.541542053223,
307.541542053223,
209.984565734863,
209.984565734863,
344.199172973633,
344.199172973633,
209.984565734863,
344.199172973633,
344.199172973633,
209.984565734863,
307.541542053223,
344.199172973633,
307.541542053223,
344.199172973633
),
Alt = c(
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_,
NA_integer_
),
Genero = c(
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Femenino",
"Femenino",
"Femenino",
"Femenino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino",
"Masculino"
),
Edad = c(
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA,
NA
),
Region = c(
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L,
6L
),
RainReg = c(
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L,
800L
),
Tipo_Suelo = c(
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico",
"Cambisol_Cromico"
),
Localidad = c(
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON",
"EL RINCON"
),
Municipio = c(
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA",
"SAN JOSE LA ARADA"
),
Department = c(
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA",
"CHIQUIMULA"
),
SOC = c(
4.01492537313432,
4.01492537313432,
4.01492537313432,
4.01492537313432,
4.01492537313432,
2.76119402985074,
2.76119402985074,
2.76119402985074,
2.76119402985074,
2.76119402985074,
2.76119402985074,
2.76119402985074,
2.76119402985074,
3.51492537313432,
3.51492537313432,
3.51492537313432,
3.51492537313432,
3.51492537313432,
2.05223880597014,
2.05223880597014,
2.05223880597014,
2.05223880597014,
5.01492537313432,
5.01492537313432,
5.01492537313432,
5.01492537313432,
5.01492537313432,
5.30597014925373,
5.30597014925373,
5.30597014925373,
5.30597014925373,
5.30597014925373,
5.35074626865671,
5.35074626865671,
5.35074626865671,
5.35074626865671,
5.35074626865671,
1.50746268656716,
1.50746268656716,
1.50746268656716,
1.50746268656716,
1.50746268656716,
2.20895522388059,
2.20895522388059,
2.20895522388059,
2.20895522388059,
2.20895522388059,
2.16417910447761,
2.16417910447761,
2.16417910447761
),
pH = c(
4.8,
4.8,
4.8,
4.8,
4.8,
4.8,
4.8,
4.8,
4.8,
4.8,
4.8,
4.8,
4.8,
5.2,
5.2,
5.2,
5.2,
5.2,
5.1,
5.1,
5.1,
5.1,
5.1,
5.1,
5.1,
5.1,
5.1,
5,
5,
5,
5,
5,
5.9,
5.9,
5.9,
5.9,
5.9,
6.4,
6.4,
6.4,
6.4,
6.4,
5.2,
5.2,
5.2,
5.2,
5.2,
5.5,
5.5,
5.5
),
Arcilla = c(
51.7,
51.7,
51.7,
51.7,
51.7,
20.2,
20.2,
20.2,
20.2,
20.2,
20.2,
20.2,
20.2,
47.5,
47.5,
47.5,
47.5,
47.5,
45.4,
45.4,
45.4,
45.4,
24.4,
24.4,
24.4,
24.4,
24.4,
24.4,
24.4,
24.4,
24.4,
24.4,
53.8,
53.8,
53.8,
53.8,
53.8,
51.7,
51.7,
51.7,
51.7,
51.7,
55.9,
55.9,
55.9,
55.9,
55.9,
34.9,
34.9,
34.9
),
P = c(
18,
18,
18,
18,
18,
95,
95,
95,
95,
95,
95,
95,
95,
8.18,
8.18,
8.18,
8.18,
8.18,
9.79,
9.79,
9.79,
9.79,
16,
16,
16,
16,
16,
20,
20,
20,
20,
20,
6.37,
6.37,
6.37,
6.37,
6.37,
18,
18,
18,
18,
18,
107,
107,
107,
107,
107,
45,
45,
45
),
Year = c(
2L,
3L,
3L,
1L,
2L,
3L,
5L,
1L,
3L,
4L,
5L,
4L,
2L,
3L,
2L,
3L,
2L,
1L,
3L,
2L,
3L,
2L,
2L,
1L,
3L,
2L,
3L,
3L,
1L,
3L,
2L,
2L,
2L,
3L,
2L,
3L,
1L,
1L,
3L,
3L,
2L,
2L,
3L,
2L,
2L,
3L,
1L,
2L,
1L,
2L
),
Tipo_Parc = c(
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Pareada",
"Pareada",
"Pareada",
"Pareada",
"Pareada",
"Pareada",
"Pareada",
"Pareada",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension",
"Extension"
),
SubDep = c(
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1",
"GT1"
)
),
row.names = c(
"1",
"2",
"3",
"4",
"5",
"6",
"7",
"8",
"9",
"10",
"11",
"12",
"13",
"14",
"15",
"16",
"17",
"18",
"19",
"20",
"21",
"22",
"23",
"24",
"25",
"26",
"27",
"28",
"29",
"30",
"31",
"32",
"33",
"34",
"35",
"36",
"37",
"38",
"39",
"40",
"41",
"42",
"43",
"44",
"45",
"46",
"47",
"48",
"49",
"50"
),
class = "data.frame"
)

3

For your DF sample:

DF[DF$ASA > 1.5 * DF$Testigo, ] |> nrow()
# 21 strictly above the 1.5 slope line

DF[DF$ASA < DF$Testigo, ] |> nrow()
# 5 strictly below the 1.0 slope line
1 Like
4

Thank you very much for the help!

5

This topic was automatically closed 7 days after the last reply. New replies are no longer allowed.

If you have a query related to it or one of the replies, start a new topic and refer back with a link.