У меня есть фрейм данных:
priors <- data.frame(dist = c('lnorm', 'beta', 'gamma'),
a = c(0.5, 1, 10),
b = c(0.4, 25, 4),
n = c(100, 100, 100)
)
, и я хотел бы взять n образцов из распределения с параметрами a и b.
Я написал эту функцию:
pr.samp <- function(n,dist,a,b) {eval (parse (
text =
paste("r",dist,"(",n,",",a,",",b,")",sep = "")
))}
Я хотел бы знать:
Спасибо за advance!