data(Overdose)
str(Overdose)
pairs(Overdose)
over.mlm <- lm(cbind(TCAD, AMI) ~ Gender + amount + BP + ECG_PR + ECG_QRS, data = Overdose)
coef(over.mlm)
# check for outliers
cqplot(over.mlm)
# HE plot shows that relations of responses to predictors are essentially one-dimensional
heplot(over.mlm)
# canonical correlation analysis
if(require(candisc)) {
cancor(cbind(TCAD, AMI) ~ as.numeric(Gender) + amount + BP + ECG_PR + ECG_QRS, data = Overdose)
}
Run the code above in your browser using DataLab