# NOT RUN {
Lt <- Liland.test(12,1,1949,162)
print(Lt)
summary(Lt)
# Critical value
Liland.crit(1, 1949, 162)
# Power
plot(Liland.pow(1,1949,161, alpha = 0.05), type = 'l', xlab = '#(x<2)', ylab = 'power')
# }
Run the code above in your browser using DataLab