inst/exercise/course/GammaExerciseCode.R

In ocrug/AdvancedRegression: Advaced Regression Models With SAS And R

QIscore.data<- read.csv(file="./GammaExerciseData.csv",
header=TRUE, sep=',')

#plotting histogram with fitted normal density
library(rcompanion)
plotNormalHistogram(QIscore.data$score)

#fitting gamma regression
summary(fitted.model<- glm(score~ desgn+wrkyrs+priorQI, 
data=QIscore.data, family=Gamma(link=log)))

#using fitted model for prediction
print(pred.score<- predict(fitted.model, type='response', 
data.frame(desgn='nurse', wrkyrs=7, priorQI='yes')))