df1 = 5
df2 = 20
x_vals= seq(0, 5, length.out = 100)
y_vals = df(x_vals, df1, df2)
f_data= data.frame(F_stat = x_vals, Density = y_vals)
alpha= 0.05
F_critical = qf(1 - alpha, df1, df2)
ggplot(f_data, aes(x = F_stat, y = Density)) +
geom_line(color = "black", linewidth = 1) +
geom_ribbon(aes(ymax = ifelse(F_stat >= F_critical, Density, NA),
ymin = 0), fill = "red", alpha = 0.5) + annotate("text",
x = F_critical - 0.5, y = 0.1, label = "Accept H0",
color = "green", size = 5) + annotate("text",
x = F_critical + 0.3, y = 0.05,
label = "Reject H0", color = "red", size = 5) +
labs(title = "F-Distribution Hypothesis Test",
x = "F Statistic", y = "P(F)")