# Error type I when we use a random sample of size 50 from a Inverse Normal distribution,
# a critical value c = 0.5 and R = 100 to test H_0: theta = 3 vs H_1: theta != 3
library(gamlss.dist)
invnormal_errorI(0.5,n=50,theta0=3,R=100)
Run the code above in your browser using DataLab