Question

当我尝试使用.rnw包中由specifyModel()指定的模型编译sem文件时，我一直收到错误消息。 knitr不喜欢F1 -> X1, lam1, NA中的逗号。

以下是重现错误的示例：

\documentclass{article}

\begin{document}

<<specify, include=FALSE, tidy=FALSE>>=
# Simple CFA Model
  # from: http://www.statmethods.net/advstats/factor.html
    library(sem)
    model.mydata <- specifyModel() 
      F1 ->  X1, lam1, NA
      F1 ->  X2, lam2, NA 
      F1 ->  X3, lam3, NA 
      F2 ->  X4, lam4, NA 
      F2 ->  X5, lam5, NA 
      F2 ->  X6, lam6, NA 
      X1 <-> X1, e1,   NA 
      X2 <-> X2, e2,   NA 
      X3 <-> X3, e3,   NA 
      X4 <-> X4, e4,   NA 
      X5 <-> X5, e5,   NA 
      X6 <-> X6, e6,   NA 
      F1 <-> F1, NA,    1 
      F2 <-> F2, NA,    1 
      F1 <-> F2, F1F2, NA

    model.mydata
# end
@

\end{document}

这是错误：

#Quitting from lines 5-27 (test.rnw) 
#Error in parse(text = x, srcfile = src) : <text>:5:16: unexpected ','
#5:     model.mydata <- specifyModel() 
#6:       F1 ->  X1,

有什么想法吗？

Answer 1

如果你没有在specifyModel中指定文件，它将尝试从标准输入流（stdin）中读取 - 虽然这可能在knitr范围内，我认为从文件中更容易阅读模型规范。

例如

model <- readLines(con = textConnection('F1 ->  X1, lam1, NA
F1 ->  X2, lam2, NA 
F1 ->  X3, lam3, NA 
F2 ->  X4, lam4, NA 
F2 ->  X5, lam5, NA 
F2 ->  X6, lam6, NA 
X1 <-> X1, e1,   NA 
X2 <-> X2, e2,   NA 
X3 <-> X3, e3,   NA 
X4 <-> X4, e4,   NA 
X5 <-> X5, e5,   NA 
X6 <-> X6, e6,   NA 
F1 <-> F1, NA,    1 
F2 <-> F2, NA,    1 
F1 <-> F2, F1F2, NA'))
cat(model, file = 'model.spec', sep = '\n')
model.mydata <- specifyModel(file = 'model.spec')

或者，您只需创建model.spec文件并将其读入即可。

sem包的knitr和specifyModel（）

1 个答案: