wd <- "C:/Users/JARVIS/Documents/R-wd/Multivariate Data Analysis/Assignment 3"
setwd(wd)
library(psych)
library(GPArotation)

data = read.table("twin.txt")
data = na.omit(data)
######### Task 2.1 #########
data_Task1 = data[2:1537,7:11]

# Convert data to numeric
indx <- sapply(data_Task1, is.factor)
data_Task1[indx] <- lapply(data_Task1[indx], function(x)
data_Task1 = as.numeric(as.character(x)))

cor(data_Task1)
factor2model <- principal(data_Task1,nfactors=2,rotate="oblimin")

loadings(factor2model)
factor2model$scores

######### Task 2.2 #########
data_Task2 = data[2:1537,4:6]

# Convert data to numeric
indx <- sapply(data_Task2, is.factor)
data_Task2[indx] <- lapply(data_Task2[indx], function(x)
data_Task2 = as.numeric(as.character(x)))

factor1model <- principal(data_Task2,nfactors=1,rotate="oblimin")

loadings(factor1model)
factor1model$scores

######### Task 2.3 #########
# Creating the 'new' data set
data_Task3 = cbind(factor1model$scores, factor2model$scores)