buildDist <- function() { pcomm <<- sample(1:5,500,TRUE,c(.1,.1,.2,.3,.3)) pprob <<- sample(1:5,500,TRUE,c(.1,.1,.1,.3,.4)) pntgrty <<- sample(1:5,500,TRUE,c(.1,.1,.1,.3,.4)) port <- .2*pcomm + .3*pprob + .5*pntgrty } op <- par(mfrow = c(2,2)) hist(pprob,main="Problem Solving",xlab="",col="Gray") hist(fitted(fit),main="Predicted Performance",xlab="",col="Gray") par(op)