semLibraryを使用したモデル分析ことはじめ

1.進捗共有

【モデル分析】

• semLibraryを使用したモデル分析ことはじめ。
• とりあえず簡単なsampleを実行してみる。

参考サイト:Rで共分散構造分析・構造方程式モデル

library(sem)

#データ読み込み
dat <- read.csv(`ltest.csv`, header=TRUE, fileEncoding = `CP932`)

#相関行列 <- とりあえずこっち使う
r <- cor(dat)
#共分散行列
cov <- cov(dat)

#モデル作成
model <- specifyModel(
text=`
keyword_category -> rank,   a1, NA
keyword_category_p1 -> rank,    a2, NA
keyword_category_p1 <-> keyword_category,   a3, NA
keyword_category <-> keyword_category,  seed1,  NA
keyword_category_p1 <-> keyword_category_p1,    seed2,  NA
rank <-> rank,  e1, 1
`)

#sem分析をかける
ans <- sem(model,r,nrow(dat))

#結果をサマる
sumally(ans)

#結果を図に出力
pathDiagram(ans, `dat`, ignore.double=FALSE, edge.labels=`values`, digits=3, node.font=c(``, 8))

【今後起きる問題】

• 日本語が文字バケる
  →外部出力して別のアプリケーションで読みこめば行けるらしい。or数名を英字表記で対応

• 複雑なモデルを簡単に出力できるような工夫が必要になりそう。
  →Excelで自動的にモデルを出力したりみたいな？