#plotting the density curve(df(x, 2, 3), from=0, to=10) curve(df(x, 3, 2), from=0, to=10) #simulating from F n.sim=100000 data=rf(n.sim, 3, 2) hist(data) summary(data) boxplot(data) #95th quantile qf(0.95, 3, 2) #say with a data set with n=20, and p=3 qf(0.95, 3, 20-3-1) curve(df(x, 3, 20-3-1), from=0, to=10)