想做結(jié)構(gòu)方程模型,卻不知用Mplus還是Amos?都不用!

大家好,我是一個研究積極心理學(xué)的教育學(xué)博士,這是我在土澳的日常

前幾天的一個晚上,需要做一個結(jié)構(gòu)方程模型(SEM)來分析一個中介變量。我打開SPSS,把數(shù)據(jù)存成.dat格式,然后準(zhǔn)備用Mplus計算模型擬合。我編好程序,點擊運行的時候突然發(fā)現(xiàn)發(fā)現(xiàn),Mplus不能用了。這種情況以前也出現(xiàn)過,所以我不慌不忙地把早就準(zhǔn)備好的安裝包拿出來重裝了一遍。可以用倒是可以用了,但系統(tǒng)提示試用版只能計算2到3個變量。我的模型有四五個變量。于是我又按照著安裝包里的破解說明,重裝了兩三遍,還是不行。由于每次重裝都要重啟,這樣折騰了一圈下來,轉(zhuǎn)眼就過去了一個小時。

我很郁悶,以前這個安裝包是可以用的,難道因為我的系統(tǒng)升級了?于是我上網(wǎng)去找新的破解版。找到幾個都有各種問題。最終,我在一個論壇上找到的付費的安裝包,要注冊那個論壇,還要付一百多塊錢!倒不是因為沒有一百塊錢,而是非常不愿意為這樣的盜版東西付錢。

但是已經(jīng)快半夜兩點了,我看著時間一點一點流逝,非常氣憤。心想:人家R那么強大的軟件包都是免費的,一個盜版的Mplus憑什么要交一百塊錢?!

誒?R是免費的?我怎么早沒想到呢?

其實今年月份申請結(jié)束以后,我報了一個大數(shù)據(jù)分析的網(wǎng)課,利用一個月的時間學(xué)R語言和python,雖然都沒有學(xué)精,但我至少電腦里裝了這些軟件,并且我知道它們都是免費的!而且有很多現(xiàn)成的分析數(shù)據(jù)的package,也都是免費的!只是我學(xué)完以后就再也沒用過,所以當(dāng)時心里還是有點沒自信。不過現(xiàn)在逼上梁山了,只能硬著頭皮上了。

于是我搜索了一下做結(jié)構(gòu)方程模型的package,找到了一個叫做lavaan的包。
接著我在youtube上找了教學(xué)軟件,盡管是英文的,不過邊看操作邊做還是能跟上。于是我又花兩個多小時,裝好了安裝包(因為R很多安裝包的版本兼容性不好,所以試了好幾個版本花了一些時間),學(xué)會了用lavaan編結(jié)構(gòu)方程模型,以及用程序直接畫圖(盡管不是很美觀)。

凌晨4點多的時候,我終于搞定了!找免費Mplus的過程非常郁悶,但是學(xué)習(xí)R的時候還是很開心的。如果早點用R的話,我應(yīng)該2點多就可以睡覺了……

當(dāng)時很激動拍了張照片

搞定之后我也有一些反思:

1. 提前學(xué)習(xí)一些技能,不知道未來什么時候就會用到。

如果我不是說在四月份學(xué)過R,而當(dāng)天晚上從頭學(xué),兩個小時應(yīng)該是搞不定這些東西的。可能甚至連R程序怎么裝都沒搞明白。恩,是的,當(dāng)時是教數(shù)據(jù)分析的老師帶著裝的,因為不僅要裝R,還要再裝一個R studio,還有各種基礎(chǔ)的數(shù)據(jù)分析的包等等。幸好這些我都提前搞定了。而且現(xiàn)學(xué)的話,那些基礎(chǔ)代碼不明白,youtube的視屏也是看不懂的,因為視頻里只是從如何使用lavaan開始的。在那么緊張的凌晨,是絕對不可能淡定的從頭學(xué)起的。總之很幸運地在幾個月前學(xué)習(xí)了R語言,雖然當(dāng)時沒有用的,但是在若干個月后的某個夜晚,凌晨2點,我竟然用到了!

2. 放下自己學(xué)不會的念頭,編程也好,任何事情都好,入門沒你想的那么難。

盡管中學(xué)時代學(xué)過一點Basic語言編程,但從本科開始就一直學(xué)心理,很久很久沒有碰過編程了。而且編程語言進化速度那么快,對程序員來說,很多時候都要重新學(xué)語言,別說我了。不過,這也說明了從頭學(xué)一個編程語言不是不可能的。
最開始聽到大數(shù)據(jù)分析什么的,我也覺得很恐懼,學(xué)心理的都懂得,統(tǒng)計是很多人不能說的傷痛。SPSS都用不熟,怎么可能還自己編程序呢?
但其實編程也不是那么難的。現(xiàn)在互聯(lián)網(wǎng)發(fā)達(dá)了,很多東西都不需要自己學(xué),都有人教,特別是入門階段,報一個輔導(dǎo)班,找一個老師帶著,很快就可以入門。剩下的內(nèi)容就可以用多少學(xué)多少,就像我學(xué)了基礎(chǔ)的R語言,然后老師也沒有教結(jié)構(gòu)方程,可是我就可以自己學(xué)了。
TED有個演講介紹20小時入門任何技能,我覺得這個還時間還是有一定準(zhǔn)確性的。

3. 不要花時間在沒有積累的事情上,而是用來學(xué)習(xí)和投資,讓努力可以積累。

同樣是2個小時,如果用來找免費資源,除了積累了一些搜索免費資源的技能,什么都學(xué)不到。而如果用這兩個小時學(xué)習(xí)如何使用R,如何做結(jié)構(gòu)方程模型,這些技能和努力在未來是可以發(fā)揮價值的。
以前互聯(lián)網(wǎng)還是免費的時代,很多人習(xí)慣了去尋找免費資源。但現(xiàn)在互聯(lián)網(wǎng)免費資源越來越少,如果你在這個時代還在找免費的資源,那浪費的時間的價值和所收獲的免費資源的價值絕對不成正比。
有句話叫做免費的才是最貴的,確實是這樣。使用免費的東西總是要付出一些代價的,可能是你為某個網(wǎng)站貢獻了流量,或者點擊了廣告。
總之,把你的時間和精力花在那些有積累效應(yīng)的事情上,而不是浪費在尋找免費資源上。

最后附上一段代碼:

***加載lavaan和畫圖的庫***

install.packages("lavaan", dependencies=TRUE)
install.packages("stringr")
install.packages("DiagrammeR")
install.packages("dplyr")
install.packages("semPlot")
install.packages("nloptr")

library("stringr")
library('lavaan')
library("DiagrammeR")
library("dplyr")
library('nloptr')
library("semPlot")

***讀入數(shù)據(jù)集***
dat1 <- read.csv("C:/Desktop/DT.csv")

***設(shè)定結(jié)構(gòu)方程模型***
Model1 <- '
grit =~ GritInterest + GritPersistence
GritInterest ~~ GritPersistence
grit ~ External + Introjected + Identified + Intrinsic + GM
External ~ GM
Introjected ~ GM
Identified ~ GM
Intrinsic ~ GM
External ~~ Introjected + Identified + Intrinsic
Introjected ~~ Identified + Intrinsic
Identified ~~ Intrinsic
'

Model2 <- '
grit =~ GritInterest + GritPersistence
grit ~ External + Introjected + Identified + Intrinsic + GM
External ~ GM
Introjected ~ GM
Identified ~ GM
Intrinsic ~ GM
'

***計算擬合***
fit1 <- sem(model=Model1,data=dat1)
fit2 <- sem(model=Model2,data=dat1)

***輸出擬合結(jié)果***
summary(fit1,standardized=T,fit.measures=T,rsq=T,modindices = TRUE)
summary(fit2,standardized=T,fit.measures=T,rsq=T,modindices = TRUE)

***畫圖***
semPaths(fit1, intercept = FALSE, whatLabel = "est",
         residuals = FALSE, exoCov = FALSE)
這個圖是另一個模型,不是上邊的代碼
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務(wù)。

推薦閱讀更多精彩內(nèi)容

  • Android 自定義View的各種姿勢1 Activity的顯示之ViewRootImpl詳解 Activity...
    passiontim閱讀 173,581評論 25 708
  • 宿柒閱讀 124評論 0 0
  • 電梯演講 電梯演講是指短時間內(nèi)表述結(jié)果內(nèi)容。 中文名 電梯演講 外文名 The elevator speech 機...
    李聘2023閱讀 5,409評論 0 0
  • 正如標(biāo)題所示,我不懂婚姻是什么,我今年24歲,和男友在一起一年半了,我們兩個發(fā)展很快,在一起一年就訂婚了,今年2-...
    無題2014閱讀 577評論 0 0
  • 那年夏天,我回到家里,他沒在家,我問娘他去哪了,娘說他去工地干活了,我又問他什么時候回來,娘沉默了一會兒,說:“快...
    牧野橋夫閱讀 258評論 2 1