r - CFA 漂亮的路径图

Question

当我从 psych 包中执行 EFA 时，我可以生成非常好的路径图，其中包含相关性和加载等:

通过以下代码使用 fa.diagram 函数:

library(psych)
covariances <- ability.cov$cov
correlations <- cov2cor(covariances)
fa.promax <- fa(correlations, nfactors=2, rotate="promax", fm="pa")
fa.diagram(fa.promax, simple=FALSE)

我更倾向于使用 CFA 而不是 EFA，并且希望能够以相同的方式绘制具有相关性的模型。我使用 SEM 包。我尝试了 pathDiagram 函数，但它只生成我可以从 Graphviz 粘贴到 GVedit 的代码。我这样做了，但相关性以及因素之间的相关性不存在。此外，我不喜欢剪切和粘贴到另一个程序(尽管来自 sem 的 pathDiagram 帮助文件显示:“要直接获取图形输出，点程序必须在系统搜索路径上。"; 我不知道到底要在实现这一目标的路径中放置什么，以及它是否还会有相关性？)

我希望能够制作一个路径图，类似于 psych 的 fa.diagram 为全民教育所做的。我希望能够在 R 中执行此操作，而不是剪切和粘贴。我不在乎您使用什么包来执行此操作，但我认为大多数人会同意越简单越好，如果可以的话在 R 中完成，它会给图形设备的选择带来很大的自由度。

以下是 CFA 模型的一些示例代码:

dat3 <- read.table(url("http://dl.dropbox.com/u/61803503/Proj_2b.dat"), 
            header=T, strip.white = T, as.is=FALSE, 
            na.strings= c("999", "NA", " "))

NHSDA.cov <- cov(dat3)

ete.mod <- specifyModel() 
F1 -> item3, lam1, NA 
F1 -> item5, lam2, NA 
F1 -> item9, lam3, NA 
F1 -> item10, lam4, NA 
F2 -> item4, lam5, NA 
F2 -> item13, lam6, NA 
F2 -> item14, lam7, NA 
F2 -> item15, lam8, NA 
F3 -> item1, lam9, NA 
F3 -> item6, lam10, NA 
F3 -> item7, lam11, NA 
F3 -> item11, lam12, NA 
F3 -> item12, lam13, NA 
F4 -> item2, lam14, NA 
F4 -> item8, lam15, NA 
F4 -> item16, lam16, NA 
F4 -> item17, lam17, NA 
item1 <-> item1, e1, NA 
item2 <-> item2, e2, NA 
item3 <-> item3, e3, NA 
item4 <-> item4, e4, NA 
item5 <-> item5, e5, NA 
item6 <-> item6, e6, NA 
item7 <-> item7, e7, NA 
item8 <-> item8, e8, NA 
item9 <-> item9, e9, NA 
item10 <-> item10, e10, NA 
item11 <-> item11, e11, NA 
item12 <-> item12, e12, NA 
item13 <-> item13, e13, NA 
item14 <-> item14, e14, NA 
item15 <-> item15, e15, NA 
item16 <-> item16, e16, NA 
item17 <-> item17, e17, NA 
F1 <-> F1, NA, 1
F2 <-> F2, NA, 1
F3 <-> F3, NA, 1
F4 <-> F4, NA, 1
F2 <-> F1, rF1F2 #1
F3 <-> F1, rF1F3 #2
F4 <-> F1, rF1F4 #3
F3 <-> F2, rF2F3 #4
F4 <-> F2, rF2F4 #5
F4 <-> F3, rF3F4 #6

ete.sem <- sem(ete.mod, NHSDA.cov, nrow(dat3)) 
(ete.SUM <- summary(ete.sem))

pathDiagram(ete.sem) #the attempt thus far

在 Win 7 上运行 R 2.14.2

Answer 1

像这样尝试:

pathDiagram(
              model=ete.sem
            , file="PathDiagram"
            , max.rank=paste("item", 1:17, sep="")
            , ignore.double=FALSE
            , edge.labels="values"
            , size=c(8, 8)
            , node.font=c("Helvetica", 10)
            , edge.font=c("Helvetica", 10)
            , rank.direction="LR"
            , digits=3
            , standardize=TRUE
            , graphics.fmt="pdf"
            )

这将在您的工作目录中生成 PathDiagram.pdf。反正我不知道在 R 中显示图形。