Skip to contents

This function will generate n random points from a gamma distribution with a user provided, .shape, .scale, and number of random simulations to be produced. The function returns a tibble with the simulation number column the x column which corresponds to the n randomly generated points, the d_, p_ and q_ data points as well.

The data is returned un-grouped.

The columns that are output are:

  • sim_number The current simulation number.

  • x The current value of n for the current simulation.

  • y The randomly generated data point.

  • dx The x value from the stats::density() function.

  • dy The y value from the stats::density() function.

  • p The values from the resulting p_ function of the distribution family.

  • q The values from the resulting q_ function of the distribution family.

Usage

tidy_gamma(
  .n = 50,
  .shape = 1,
  .scale = 0.3,
  .num_sims = 1,
  .return_tibble = TRUE
)

Arguments

.n

The number of randomly generated points you want.

.shape

This is strictly 0 to infinity.

.scale

The standard deviation of the randomly generated data. This is strictly from 0 to infinity.

.num_sims

The number of randomly generated simulations you want.

.return_tibble

A logical value indicating whether to return the result as a tibble. Default is TRUE.

Value

A tibble of randomly generated data.

Details

This function uses the underlying stats::rgamma(), and its underlying p, d, and q functions. For more information please see stats::rgamma()

Author

Steven P. Sanderson II, MPH

Examples

tidy_gamma()
#> # A tibble: 50 × 6
#>    sim_number     x       y      dx      dy       p
#>    <fct>      <int>   <dbl>   <dbl>   <dbl>   <dbl>
#>  1 1              1 0.238   -0.254  0.00518 0.238  
#>  2 1              2 0.175   -0.215  0.0199  0.175  
#>  3 1              3 0.0363  -0.177  0.0630  0.0363 
#>  4 1              4 0.187   -0.139  0.167   0.187  
#>  5 1              5 0.394   -0.101  0.369   0.394  
#>  6 1              6 0.0175  -0.0624 0.690   0.0175 
#>  7 1              7 0.526   -0.0241 1.09    0.526  
#>  8 1              8 0.229    0.0142 1.48    0.229  
#>  9 1              9 0.0430   0.0525 1.76    0.0430 
#> 10 1             10 0.00369  0.0907 1.87    0.00369
#> # ℹ 40 more rows