R/共分散構造分析/実験

原因２つが無相関

データ <- readMoments(diag=TRUE,names=c("原因1", "原因2", "結果"))
1
0.0  1
0.5  0.8 1
model <- specifyModel()
原因1 -> 結果, path1, NA
原因2 -> 結果, path2, NA
原因1 <-> 原因1, NA, 1
原因2 <-> 原因2, NA, 1
結果 <-> 結果, er, NA
ans <- sem(model, データ, N=100)
summary(ans)
pathDiagram(ans, ignore.double=FALSE, edge.labels="values", digits=3)

a graph image

原因２つに相関あり

データ <- readMoments(diag=TRUE,names=c("原因1", "原因2", "結果"))
1
0.3 1
0.5 0.8 1
model <- specifyModel()
原因1 -> 結果, path1, NA
原因2 -> 結果, path2, NA
原因1 <-> 原因1, NA, 1
原因2 <-> 原因2, NA, 1
結果 <-> 結果, er, NA
ans <- sem(model, データ, N=100)
summary(ans)
pathDiagram(ans, ignore.double=FALSE, edge.labels="values", digits=3)

a graph image

model <- specifyModel()
原因1 -> 結果, path1, NA
原因2 -> 結果, path2, NA
原因1 <-> 原因2, cor1, 1
結果 <-> 結果, er, NA

a graph image

model <- specify.model()
原因1 -> 結果, path1, NA
原因2 -> 結果, path2, NA
原因1 -> 原因2, path3, NA
原因1 <-> 原因1, NA, 1
原因2 <-> 原因2, NA, 1
結果 <-> 結果, er, NA

a graph image

model <- specify.model()
原因1 -> 結果, path1, NA
原因2 -> 結果, path2, NA
原因1 -> 原因2, path3, NA
原因1 <-> 原因1, NA, 1
結果 <-> 結果, er, NA

a graph image