R for Statistical Inference course Question 5(b)

# Set the range of p-values
p_values <- seq(0.001, 0.1, by = 0.001)

# Calculate the corresponding x values for a standard normal distribution
x_values <- qnorm(1 - p_values)

# Calculate the likelihood ratio for each p-value
lr_values <- exp(-0.5 * x_values^2)

# Plot the likelihood ratio
plot(p_values, lr_values, type = "b", log = "y", xlab = "p-value", ylab = "Likelihood Ratio (log scale)",
     main = "Likelihood Ratio vs. p-value")