主要依據《r實戰》,用自己的資料做探索,反正這文章也沒人看哈哈哈。
涉及到的包:psych、gparotation
#載入資料
data=read.csv("d:/from r/sem in r/study1.csv",header = true,sep=",")
data
str(data)
#生成相關係數矩陣
cor
library(psych)
#判斷需要提取的公共因子數,生成碎石圖
fa.parallel(cor,n.obs=209,fa="both",n.iter = 100,
main = "scree plots with parrallel analysis")
#修改x軸顯示引數
axis(1,at=seq(1,30,1))
#主軸迭代法提取因子,並進行斜交旋轉
library(gparotation)
fa.promax
fa.promax
#獲得因子荷載矩陣
fsm
else
}fsm(fa.promax)
#繪製正交或斜交結果的圖形
factor.plot(fa.promax,labels = rownames(fa.promax$loadings))
fa.diagram(fa.promax,******=f)
最後我的結果是這樣的:
很奇怪?和理論預設的不太一樣。。
並且按照《r實戰》的流程並沒有像mplus的結果一樣報告很多東西....
結構方程模型和r語言還要再繼續學習.....
R 語言 因子分析
因子分析 options digits 2 covariances ability.cov cov 將協方差矩陣轉化為相關係數矩陣 correlations cov2cor covariances 第一步 判斷需提取的公共因子數 library psych covariances ability.c...
資料的探索性分析及R實現
一.讀入資料city read.csv 2015city.csv 二.描述性統計data.outline function x input x a vector output the outline of data x 三.列聯表library vcd table a 用table 函式生成簡單的頻...
資料的探索性分析
1.簡單統計量分析 極差,最大值,最小值 2.3 原則 如果資料服從正態分佈,在3 原則下,異常值被定義為一組測定值中與平均值的偏差超過3倍標準差的值。如果資料不服從正態分佈,也可以用標準差作為 來進行倍數描述 3.箱型圖分析 在箱型圖中,異常值被定義為小於ql 1.5iqr或大於qu 1.5iqr...