# MANOVA with R and likeability data
library(foreign)
datafile <- read.spss(file.choose(),to.data.frame=TRUE)
# names(datafile)
# [1] "ParticipantNumber" "Condition" "Awkward"
# [4] "Inattentive" "Involved" "Bored"
# [7] "Friendly" "Cold" "Trustworthy"
#[10] "Positive" "Incompetent" "Perceptive"
#[13] "Generalimpression" "Similar" "Ratingofother"
#[16] "PartnerEvalYou"
datafile$Trustworthy[datafile$Trustworthy == -99] <- NA
attach(datafile)

Y <- cbind(Awkward,Inattentive,Involved,Bored,Friendly,Cold,Trustworthy,Positive,
Incompetent,Perceptive,Generalimpression,Similar,Ratingofother,PartnerEvalYou)
manatee <- manova(Y ~ Condition)
summary(manatee)

library(MASS)
(manatee2 <- lda(Condition ~ Y))
preds <- predict(manatee2)

plot(preds$x,col=as.numeric(Condition)+1,pch=19,xlab="LD1 86%",
ylab="LD2 14%",cex.lab=1.3,cex.axis=1.3,
main="LDs by Group and People")
meanX <- tapply(preds$x[,1],as.numeric(Condition[-2]),mean)
meanY <- tapply(preds$x[,2],as.numeric(Condition[-2]),mean)
points(meanX,meanY,cex=3,pch=8,col=2:4)
legend("bottomleft",pch=8,col=2:4,legend=c("Control","Like","Dislike"),
cex=1.3,text.col=2:4)

plot(manatee2$scaling[,1],manatee2$scaling[,2],pch=19,xlab="LD1 86%",
ylab="LD2 14%",cex.lab=1.3,cex.axis=1.3,main="LDs by Variables")
text(-.6,-.55,"Friendly")
text(.27,-.55,"Awkward")