inst/examples/R/stat/Quantil/Lampen2.R

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

gl <- read.csv("gluehlampen.csv")
F_hat = c(0, 1, 7, 16, 19, 19)
plot(c(gl$k_u, gl$k_o[4], 2500), F_hat, type = "l", axes = FALSE, ylab = "",
     xlab = "", xpd = TRUE, lwd = 3)
axis(1, at =c(1:2500), labels = rep("",2500), xpd = TRUE, tck = 0, pos = 0)
axis(2, at = c(0:21), labels = rep("",22), pos = 0, xpd = TRUE, tck = 0)
points(0, 21, pch =17, xpd = TRUE)
text(x = 666.67, y = -1, labels=expression(x["0,5"]), xpd = TRUE)
text(x = 2500, y = -1, labels="x", xpd = TRUE)
text(x = -80, y = 22, labels= "F(x)", xpd = TRUE)
text(x = -80, y = 10, labels= "0,5", xpd = TRUE)
text(x = -80, y = 19, labels= "1", xpd = TRUE)
lines(x = c(666.67,666.67), y = c(0, 10), lty=2)
lines(x = c(0,666.67), y = c(10, 10), lty=2)
lines(x = c(0,2500), y = c(19, 19), lty=2, xpd = TRUE)