Skip to contents

This function will generate n random points from a chisquare distribution with a user provided, .df, .ncp, 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_chisquare(
  .n = 50,
  .df = 1,
  .ncp = 1,
  .num_sims = 1,
  .return_tibble = TRUE
)

Arguments

.n

The number of randomly generated points you want.

.df

Degrees of freedom (non-negative but can be non-integer)

.ncp

Non-centrality parameter, must be non-negative.

.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::rchisq(), and its underlying p, d, and q functions. For more information please see stats::rchisq()

Author

Steven P. Sanderson II, MPH

Examples

tidy_chisquare()
#> # A tibble: 50 × 7
#>    sim_number     x       y      dx      dy     p       q
#>    <fct>      <int>   <dbl>   <dbl>   <dbl> <dbl>   <dbl>
#>  1 1              1  9.62   -2.08   0.00142 0.982  9.62  
#>  2 1              2  0.884  -1.74   0.00578 0.450  0.884 
#>  3 1              3  1.37   -1.40   0.0189  0.553  1.37  
#>  4 1              4  2.95   -1.07   0.0498  0.760  2.95  
#>  5 1              5  0.345  -0.731  0.106   0.284  0.345 
#>  6 1              6  1.76   -0.395  0.182   0.618  1.76  
#>  7 1              7  0.0717 -0.0584 0.257   0.130  0.0717
#>  8 1              8  0.612   0.278  0.303   0.376  0.612 
#>  9 1              9 12.3     0.614  0.306   0.994 12.3   
#> 10 1             10  3.39    0.950  0.274   0.797  3.39  
#> # ℹ 40 more rows