Hi @Nuria
this is the code i used
ESHW9 <- integrate(gxhn,
0,
500,
sigma = DetectteamB@estimate)
ESHW9
ESHW9$value / 500
ESHW10 <- integrate(gxhn,
0,
500,
sigma = DetectteamA@estimate)
ESHW10
ESHW10$value / 500