Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- > nuclear <- c(7, 20, 16, 6, 58, 9, 20, 50, 23, 33, 8, 10, 15, 16, 104)
- > quantile(nuclear)
- 0% 25% 50% 75% 100%
- 6.0 9.5 16.0 28.0 104.0
- quart <- function(x) {
- x <- sort(x)
- n <- length(x)
- m <- (n+1)/2
- if (floor(m) != m) {
- l <- m-1/2; u <- m+1/2
- } else {
- l <- m-1; u <- m+1
- }
- c(Q1=median(x[1:l]), Q3=median(x[u:n]))
- }
- Q1 Q3
- 9 33
- y <- matrix(NA, 2, 10)
- rownames(y) <- c("Q1", "Q3")
- colnames(y) <- c(1:9, "Quart")
- for (n in 3:5) {
- j <- 1
- for (i in 1:9) {
- y[, i] <- quantile(1:n, probs=c(1/4, 3/4), type=i)
- }
- y[, 10] <- quart(1:n)
- cat("n", n, ":n")
- print(y, digits=2)
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement