# family of distributions parameterized by TWO df1 values f_values <- rf(n = 10000, df1 = 15, df2 = 50) hist(f_values, breaks = 30)
# normalize distribution for QQ f_values <- (f_values - mean(f_values))/sd(f_values) qqnorm(f_values); qqline(f_values, col = 2)