Posted at

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

More than 3 years have passed since last update.


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))

151116163721_Rplot.png


【今後起きる問題】


  • 日本語が文字バケる

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


  • 複雑なモデルを簡単に出力できるような工夫が必要になりそう。

    →Excelで自動的にモデルを出力したりみたいな?