Skip to contents

This function will generate n random points from a paralogistic distribution with a user provided, .shape, .rate, .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_paralogistic(
  .n = 50,
  .shape = 1,
  .rate = 1,
  .scale = 1/.rate,
  .num_sims = 1
)

Arguments

.n

The number of randomly generated points you want.

.shape

Must be strictly positive.

.rate

An alternative way to specify the .scale

.scale

Must be strictly positive.

.num_sims

The number of randomly generated simulations you want.

Value

A tibble of randomly generated data.

Details

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

Author

Steven P. Sanderson II, MPH

Examples

tidy_paralogistic()
#> # A tibble: 50 × 7
#>    sim_number     x       y    dx       dy      p       q
#>    <fct>      <int>   <dbl> <dbl>    <dbl>  <dbl>   <dbl>
#>  1 1              1  0.434  -2.40 1.59e- 3 0.302   0.434 
#>  2 1              2  0.251   1.45 2.08e- 1 0.200   0.251 
#>  3 1              3 10.2     5.29 3.47e- 2 0.911  10.2   
#>  4 1              4  4.58    9.14 5.44e- 3 0.821   4.58  
#>  5 1              5  0.326  13.0  7.02e- 4 0.246   0.326 
#>  6 1              6  3.09   16.8  6.87e- 4 0.756   3.09  
#>  7 1              7  0.242  20.7  4.79e- 3 0.195   0.242 
#>  8 1              8  0.634  24.5  3.56e-10 0.388   0.634 
#>  9 1              9  0.0413 28.4  0        0.0397  0.0413
#> 10 1             10  0.369  32.2  1.47e-18 0.269   0.369 
#> # ℹ 40 more rows