Generates a random sample of y scores and x scores from a bivariate normal distributions with specified population means, standard deviations, and correlation. This function is useful for generating hypothetical data for classroom demonstrations.
random.yx(n, my, mx, sdy, sdx, cor, dec)
sample size
population mean of y scores
population mean of x scores
population standard deviation of y scores
population standard deviation of x scores
population correlation between x and y
number of decimal points
Returns n pairs of y and x scores
random.yx(10, 50, 20, 4, 2, .5, 1)
#> y x
#> 1 46.6 19.1
#> 2 52.6 23.4
#> 3 49.2 20.5
#> 4 47.3 20.2
#> 5 46.4 19.6
#> 6 45.7 21.3
#> 7 46.7 21.6
#> 8 56.4 23.2
#> 9 48.1 19.9
#> 10 51.9 20.6
# Should return:
# y x
# 1 50.3 21.6
# 2 52.0 21.6
# 3 53.0 22.7
# 4 46.9 21.3
# 5 56.3 23.8
# 6 50.4 20.3
# 7 44.6 19.9
# 8 49.9 18.3
# 9 49.4 18.5
# 10 42.3 20.2