load("C:/Users/Bruno/Livre AQR/Data SMP/dataAQR")
mhp.mod$direct <- 2-mhp.mod$determ
mhp.mod$direct.1 <- mhp.mod$determ.1
mhp.mod$direct.2 <- mhp.mod$determ.2
mhp.mod$direct.3 <- mhp.mod$determ.3
mhp.mod$direct.4 <- mhp.mod$determ.4
mhp.mod$direct.5 <- mhp.mod$determ.5
ns <- c("ns1","ns2","ns3","ns4","ns5")
mhp.mod$score.ns <- apply(mhp.mod[,ns],1,sum)
sum(is.na(mhp.mod$ns))


charac <- c("direct.1","direct.2","direct.3","direct.4","direct.5","coop.1","coop.2",
"coop.3","coop.4","coop.5","emostab.1","emostab.2","emostab.3","emostab.4",
"emostab.5","transc.1","transc.2","transc.3","transc.4","transc.5")

temp <- c("ns1","ns2","ns3","ns4","ns5","ha1","ha2","ha3","ha4","ha5","rd1","rd2","rd3","rd4","rd5")
par(mfrow=c(3,5))

for(i in temp) barplot(table(mhp.mod[,i],useNA="ifany"), main=i, space=0.6, names.arg=c("1","2","3","NA"), ylim=c(0,500), col=c("white","white","white","black"))


library(psy)
scree.plot(mhp.mod[,charac], simu=20)
scree.plot(mhp.mod[,temp], simu=20)

fa.res1 <- factanal(na.omit(mhp.mod[,temp]), factors=3)
fa.res1

library(psy)
cronbach(mhp.mod[,charac])












#######################################################################################
ns <- c("ns1","ns2","ns3","ns4","ns5")
ha <- c("ha1","ha2","ha3","ha4","ha5")
rd <- c("rd1","rd2","rd3","rd4","rd5")
par(mfrow=c(1,3))
mtmm(mhp.mod,list(ns,ha,rd), itemtot=TRUE, namesDim=c("Novelty Seeking","Harm Avoidance","Reward Dependence"))
par(mfrow=c(1,3))
mtmm(mhp.mod,list(ns,ha,rd), itemtot=TRUE, stripChart=TRUE, namesDim=c("Novelty Seeking","Harm Avoidance","Reward Dependence"))

direct <- c("direct.1","direct.2","direct.3","direct.4","direct.5")
coop <- c("coop.1","coop.2","coop.3","coop.4","coop.5")
emostab <- c("emostab.1","emostab.2","emostab.3","emostab.4","emostab.5")
transc <- c("transc.1","transc.2","transc.3","transc.4","transc.5")
par(mfrow=c(1,4))
mtmm(mhp.mod,list(direct,coop,emostab,transc), itemTot=TRUE, namesDim=c("Self-directedness","Cooperativeness","Affective Stability","Self-transcendence"))