匯入SPSS資料檔並實作因素分析(2)

• 為了要找出各變項在各因素之組型負荷量（pattern loadings），因此我們以psych套件中的factor.pa()函數來計算。未轉軸為fpex_n，使用varimax法轉軸為fpex_v。（下方的指令同時輸入估計分數）

fpex_n<-factor.pa(na.omit(ex3_2[,2:17]),3,rotate="none",scores=T,residuals=T)
fpex_v<-factor.pa(na.omit(ex3_2[,2:17]),3,rotate="varimax",scores=T,residuals=T)

• 未轉軸的組型負荷量。

• 使用varimax轉軸的組型負荷量。

• 和之前結果相同也就是使用3個因素可以解釋67%的變異量。
• 最後我們將轉軸後的估計分數加以輸出到新資料檔中以利後續的使用。

write.table(fpex_v$scores,"fpex_vscore.csv",sep=",",col.names=NA)