日々の記録簿

日々のできごとの記録

日記

Multivariate Regression Analysis に関するメモ

投稿日:

Multivariate Regression Analysis (多変量回帰分析)に関するメモです。

多変量回帰にもいろいろあるようですが、ここでは「独立変数は同じだが従属変数が複数あるケース」を想定しています。

統計ソフトによる分析例

Mplus

Mplus Data Analysis Examples Multivariate Regression Analysis

Mplusによる分析例と検定の方法が掲載されている。

多変量回帰分析のイメージをつかむためには有効。ただし、分析方法しか載っていないので、背後にある理屈については別に勉強する必要がある。

Stata

Stata Data Analysis Examples Multivariate Regression Analysis

Stata版。先ほどと同じデータを使っている。Mplusよりこちらのほうが手順がやや丁寧。ただし、やはり細かい理屈は説明されていない。

R

Multivariate Linear Models in R(PDF)

Rによる解説。ある程度、数式も使って説明されている。

Multivariate multiple regression

Rによるコード。ほとんど説明はない。

Multivariate multiple regression in R

manovaとcar package内にあるManovaを使用した際に、なぜ出力がちがうのかを説明している。上記のサイトでサラッと流されている点を丁寧に解説してくれている。

Rによる分析

RでStataやMplusによる分析結果を再現してみる。

# データの準備
mvreg <- read.table("http://www.ats.ucla.edu/stat/mplus/dae/mvregex.dat", sep=",")
names(mvreg) <- c("locus", "self", "motiv", "read", "write", "science", "prog", "prog1", "prog2", "prog3")
mvreg$prog <- factor(mvreg$prog1 + 2 * mvreg$prog2 + 3 * mvreg$prog3)

# car packageを使う
library(car)

# 多変量回帰
mod.mvreg <- lm(cbind(locus, self, motiv) ~ read + write + science + prog, data=mvreg)

# 係数の出力
mod.mvreg
summary(mod.mvreg)

# 検定1
manova.mvreg <- Anova(mod.mvreg)
summary(manova.mvreg)

# 検定2。1と同じになる
summary(Manova(mod.mvreg))

linearHypothesis(mod.mvreg, "read=0", verbose=TRUE)

文献

細かい理屈がわからないので、いろいろ読んでみる。感想は日をあらためて書く予定。

スポンサーリンク

スポンサーリンク

-日記

執筆者:


comment

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

関連記事

vim-r-pluginの導入方法

vimからRを使う、vim-r-pluginの紹介。 まあ、vimにこだわりがないなら、RStudio使ったほうが幸せになれると思います。 環境は、Arch Linuxです。 インストール vim-r …

FactoMineRからベンゼクリの修正割合(Benzecri’s modified rates)を計算する

最近、多重対応分析(Multiple Correspondence Analysis)を勉強していた。RではFactoMineRのMCA関数がいい感じである。 ただ、ベンゼクリの修正割合を計算できると …

Thinkpad X250を購入した!

物欲に負けて、Thinkpad X250を購入してしまいました・・・。 さっそく、Linuxをインストールしたのですが、いくつか設定が必要なことがあったので、そのメモです。 トラックポイント/タッチパ …

Nvim-Rの導入方法(その2)

vimからRを使う場合、Nvim-Rが便利です。 以前にも紹介しましたが、設定をいくつか間違えていたので修正版です。 インストール まず、neovimをインストールする。このとき、deopleteを使 …

ThinkPad X250 のトラックポイントが使いにくい!

X250を1ヶ月ほど使っていますが、最近いまいち旧機種に比べて使いにくいと思うようになりました。 で、気づいたのが、X250はトラックポイントの反応が鈍いということ。旧機種に比べて、操作に力が要ります …