Вам может понадобиться erfinv
из pracma
, чтобы получить среднее значение mu
, например,
p <- 0.02
x <- 400
sgm <- 4
mu <- x - sgm*sqrt(2)*pracma::erfinv(2*p-1)
, такое, что
> mu
[1] 408.215
> pnorm(400,mu,sgm) # check the obtained value of mu
[1] 0.02
Более разумный подход - через qnorm
(спасибо @ user20650)
mu <- qnorm(p,x,sgm,FALSE)