inst/examples/R/data/theory/statistical_tests/fluctuation.R

In sigbertklinke/shinyExample: Create and run simple (interactive) programs as Shiny apps

x <- (0:1500)/100
y <- dchisq(x, 5)
cl <- 1+trunc(100*qchisq(0.05, 5))
cu <- 1+ceiling(100*qchisq(0.95, 5))          

pdf("fluctuation.pdf", width=5, height=2)
par(mar=c(0,0,3,0))
plot(x, y, type="l", axes=F, ylim=c(-0.05, max(y)), main=expression(paste(1-alpha, " fluctuation interval")))
polygon(c(x[cl:cu],x[cu:cl]), c(y[cl:cu],rep(0, cu-cl+1)), col="grey")
usr <- par("usr")
lines(usr[1:2], c(0,0))
lines(c(3,3), c(0, max(y)))
text(3,0, labels=expression(theta), pos=1)
lines(c(x[cl], x[cl]), c(0, y[cl]))
lines(c(x[cu], x[cu]), c(0, y[cu]))
text(x[cl],0, labels=expression(theta-c[l]), pos=1)
text(x[cu],0, labels=expression(theta+c[u]), pos=1)
text(x[cl], y[cl], labels=expression(alpha/2), pos=2)
text(x[co], y[cl], labels=expression(alpha/2), pos=4)
text(3, y[cl], labels=expression(1-alpha), pos=4)
dev.off()