R/fp1.R

In austinragotzy/mathstat: Introduction to Mathematical Statistics R Package

#' @title F-Distribution Table for Upper 0.05 Critical Points
#'
#' @description Produces the F-distribution table for the 0.95 quantile
#' of the F-distribution.
#'
#' @details This table more specifically produces the F-distribution table for
#' when alpha equals 0.05. The initial row indicates the first nine degrees
#' of freedom in the numerator and the inital column indicates degrees of freedom in
#' the denominator. See page 711 for more details.
#'
#' @return \code{xmat} is the F-distribution table for the 95th quantile.
#'
#' @references Hogg, R., McKean, J., Craig, A. (2018) Introduction to Mathematical
#'             Statistics, 8th Ed. Boston: Pearson.
#'
#' @seealso fp2() for the F-distribution table of the 0.95 quantile with the numerator
#' having greater than nine degrees of freedom., as well as fp3() and
#' fp4() for the F-distribution table of the 0.99 quantile.
#'
#' @examples
#' fp1()
#'
#' @export fp1

fp1 <- function() {
  
  numdf <- 1:9
  denomdf <- c(1:30, 40, 60, 120)
  
  ic <- length(denomdf)
  
  xmat <- matrix(rep(0, 9 * ic), ncol = 9)
  
  for (i in 1:ic) {
    for (j in 1:9) {
      xmat[i, j] <- round(100 * qf(0.95, j, denomdf[i]))/100
    }
  }
  
  xmat <- cbind(denomdf, xmat)
  xmat <- rbind(c(0, numdf), xmat)
  
  vec <- c(0, round(100 * qchisq(0.95, numdf)/numdf)/100)
  
  xmat <- rbind(xmat, vec)
  
  return(xmat)
}