mydata <- read.table("bigfive.csv", sep=",", header=TRUE) mydata[,1] mydata[,2] cor(mydata[,1], mydata[,2]) cor(mydata) mycorrelation <- cor(mydata) library(psych) oitems <- mydata[ , c(1, 2, 3, 4, 5, 6) ] oitems <- mydata[ , 1:6 ] alpha( oitems ) ototal <- oitems[,1] + oitems[,2] + oitems[,3] + oitems[,4] + oitems[,5] + oitems[,6] oitems2 <- data.frame(oitems, ototal)