南京邮电大学多元统计分析实验五-因子分析

使用平台：R、Rstudio

题目描述：分析评价全国35个中心城市的综合发展水平，无论是对城市自身的发展，还是对周边地区的进步，都具有十分重要的意义。应用因子分析模型，选取反应城市综合发展水平的12个指标作为原始变量，对全国35个中心城市的综合发展水平作分析评价。

> rm(list=ls())
> library(haven)      #用于read_csv函数
> library(tidyverse)  # 用于数据处理
> library(psych)      # 用于因子分析> data<-read_sav("因子分析（City）.sav")
> data=data[,-1]
> # 进行因子分析
> fa_results <- fa(data, 12, rotate = "varimax")  # 对12个指标进行因子分析，使用varimax旋转
> # 查看因子分析结果
> print(fa_results)
> factor_scores <- factor.scores(data, fa_results$loadings)  # 提取各城市在各因子上的得分
> # 打印各城市因子得分
> print(factor_scores)
> barplot(t(factor_scores$scores), beside = TRUE, col = rainbow(12), 
+         xlab = "Cities", ylab = "Factor Scores", 
+         main = "Factor Scores for 35 Cities")