str(Rehab)
attach(Rehab)
qqnorm(differ)
qqline(differ)
shapiro.test(differ)
boxplot(Psych1,Psych2,names=c("Psychiatrist 1","Psychiatrist 2"),
col=c("pink","lightblue"))
t.test(Psych1,Psych2,paired=TRUE)
detach(Rehab)
Run the code above in your browser using DataLab