R/proc_glm.R

In procs: Recreates Some 'SAS®' Procedures in 'R'

# General Linear Model ----------------------------------------------------



#
# proc_glm <- function() {
#
#
#
#
#
# }



# @import sasLM
# get_lm_stats <- function() {
#
#
#
#
# }

# library(sasLM)
#
# BEdata = af(BEdata, c("SEQ", "SUBJ", "PRD", "TRT")) # Columns as factor
# formula1 = log(CMAX) ~ SEQ/SUBJ + PRD + TRT # Model
# GLM(formula1, BEdata) # ANOVA tables of Type I, II, III SS
# EMS(formula1, BEdata) # EMS table
# T3test(formula1, BEdata, Error="SEQ:SUBJ") # Hypothesis test
# ci0 = CIest(formula1, BEdata, "TRT", c(-1, 1), 0.90) # 90$ CI
# exp(ci0[, c("Estimate", "Lower CL", "Upper CL")]) # 90% CI of GMR