#Packages to install: Rcmdr, simpleboot, aplpack, MASS(VR), HH
#Bootstrap Review
x=rnorm(1000)
mean(x)
sd(x)

library(simpleboot) 
b=one.boot(x, mean, R=100)
hist(b,do.rug=T)
 
#Robust Correlation The following assumes the work ethic and publications has been imported
library(aplpack) 
attach(workethic) #workethic needs to be whatever you named the data upon import
bagplot(WORK_ETH,NUM_PUBS)
library(MASS)
cor(workethic)
cov.rob(workethic, cor=T,method="mcd")$cor
 
#Regression
mymodel = lm(NUM_PUBS~WORK_ETH) #mymodel is an arbitrary name
summary(mymodel)
 
#This alternative approach will automatically spit out the summary
(summary(lm(NUM_PUBS~WORK_ETH)))

#Obtain predicted values and plot intervals
mypredvalues=predict(mymodel, interval = "prediction") 
mypredvalues
library(HH)
ci.plot(mymodel)