Differences

This shows you the differences between two versions of the page.

--- r:path_analysis [2022/11/15 22:17] – [Lavaan in R: explanation] hkimscil
+++ r:path_analysis [2023/05/29 23:45] – hkimscil
@@ Line 1: / Line 1: @@
 ====== Path Analysis ======
+{{:r:pasted:20230529-234519.png}}
 ====== Introduction ======
 {{youtube>UGIVPtFKwc0}}
@@ Line 114: / Line 115: @@
 # my own
 # pbt model
-specmod5 <- "
+specmod5 <- '
     # Directional relations (path)
-    intention ~ attitude + norms + control
+    intention ~ a*attitude + b*norms + c*control
-    behavior ~ intention + norms
+    behavior ~ d*intention
     # Covariances
     attitude ~~ norms + control
     norms ~~ control
-"
+    ad := a*d
+    bd := b*d
+    cd := c*d
+'
 fitmod5 <- sem(specmod5, data=df)
 summary(fitmod5, fit.measures=TRUE, rsquare=TRUE)
@@ Line 747: / Line 751: @@
 CODING
 <code>
-processdata <- read.csv("http://commres.net/wiki/_media/r/path_analysis_datan_binw.csv")
+processdata<-read.csv("http://commres.net/wiki/_media/r/path_analysis_datan_binw.csv",
+                       header=TRUE, sep=",", fileEncoding="UTF-8-BOM")
 str(processdata)
 library(lavaan)
@@ Line 910: / Line 915: @@
 see [[https://www.rensvandeschoot.com/tutorials/lme4/|lme4 tutorial]]
+===== Exercise =====
+Using mtcars in R
+<code>
+?mtcars
+mtcars
+str(mtcars)
+df <- mtcars
+</code>
+<code>
+# model specfication
+model <-'
+  mpg ~ hp + gear + cyl + disp + carb + am + wt
+  hp ~ cyl + disp + carb
+'
+# model fit
+fit <- cfa(model, data = mtcars)
+summary(fit, fit.measures = TRUE, standardized=T, rsquare=T)
+semPaths(fit, 'std', layout = 'circle')
+</code>