# Spring 2006 R code for FACTOR ANALYSIS 06sR91.txt 4/11/06

T84=matrix(scan(file='D:/T8-4.dat'),byrow=T,ncol=5)

T84[1,]

n=dim(T84)[1]

p=dim(T84)[2]

#dimnames(T84)=list(NULL,c("AlliedChemical","DuPond","UnionCarbide","Exxon","Texco"))

#-------------------1 factor-----------

T84.fa=factanal(T84, factors=1)

#MLE and VARIMAX are default

print(T84.fa)

#T84.fa=factanal(T84,factors=1,method="mle");print(T84.fa);print(summary(T84.fa))

#plot(1:n,T84.fa$scores[,1]); plot(1:n,T84.fa$scores[,1])

factanal(T84,factors=1,scores="Bartlett")$scores

#-------------------2 factors--------

T84.fa=factanal(T84,factors=2,scores="regression")

print(T84.fa)

plot(loadings(T84.fa), type="n")

text(loadings(T84.fa)[,1],loadings(T84.fa)[,2],c("AlliedChemical","DuPond","UnionCarbide","Exxon","Texco"))

plot(1:n,T84.fa$scores[,1],xlim=c(-1,101), ylim=c(-3,3), pch=19,col="blue",ylab="scores for Factor 1")

plot(T84.fa$scores, main="Scatter plot of factor 1 and factor 2",xlab="factor 1",ylab="factor 2",pch=19,col="red")

T84.fa=factanal(T84,factors=2,method="mle",rotation="varimax",scores="regression")

print(T84.fa)

plot(loadings(T84.fa), type="n")

text(loadings(T84.fa)[,1],loadings(T84.fa)[,2],c("AlliedChemical","DuPond","UnionCarbide","Exxon","Texco"))