SEMに関する覚え書き

RでSEMをする際の覚え書き。だいたい理解できた。

Sys.setenv("http_proxy"="http://address:10080")
source("all.R")
factanal2(dataset[c(21,22,23,24,25,26,27,28)],factors=3,rotation="promax")

cor1 <- cor(dataset[c(85,86,87,94,95,96,97,98,99,100,101,102)])
model1 <- specify.model()
ans1 <- sem(model1, cor1, N=847)
summary(ans1)
std.coef(ans1)
path.diagram(ans1, ignore.double=FALSE, edge.labels="values", digits=3, standardize=TRUE)

Graphvizでノードをまとめる方法。

graph [label = "キャプションを書く", labelloc = t];
{ rank = same; "AT" "NEX" "PEX"}
{ rank = min; "Auniv" "Cuniv" "Huniv"}
{ rank = max; "exC" "exN1"}