add a function to sample from posterior

[annecori]

sample_posterior_R <- function(R, n_simul = 1000)
{
mu <- R$R$`Mean(R)`
sigma <- R$R$`Std(R)`
Rshapescale <- epitrix::gamma_mucv2shapescale(mu = mu, cv = sigma / mu)
rgamma(n_simul, shape = Rshapescale$shape, scale = Rshapescale$scale)
}

[zkamvar]

Given that epitrix::gamma_mucv2shapescale() looks like this:

epitrix::gamma_mucv2shapescale
#> function (mu, cv) 
#> {
#>     shape <- 1/(cv^2)
#>     scale <- mu * cv^2
#>     return(list(shape = shape, scale = scale))
#> }
#> <bytecode: 0x5629fdc76470>
#> <environment: namespace:epitrix>

I would suggest that we use this formulation instead of importing epitrix.