1

Stata と R の結果を比較しようとすると混乱します。Web ページhttp://www.ats.ucla.edu/stat/stata/webbooks/logistic/chapter2/default.htmにある例を使用しています。 最初にStataで次のコマンドを実行します

use http://www.ats.ucla.edu/stat/stata/webbooks/logistic/apilog, clear

次に、セクションに記載されている次のコマンドを使用します (2.2.2 主な効果と相互作用を伴う 2 x 2 レイアウト)

generate cred_ed = cred_hl*pared_hl
logit hiqual cred_hl pared_hl cred_ed

これら 2 つのコマンドは、Web ページに表示される結果を生成します。

そして、次のRコードを使用して同じ例を再現しました

Data<- read.csv("Book1.csv",header=T)
data.glm<-glm(hiqual~cred_hl + pared_hl + cred_hl*pared_hl,family=binomial,  data=Data)
summary(data.glm)

しかし、結果は一致しません！

R のデータ ファイルは、次のリンクからダウンロードできます。

https://spreadsheets.google.com/spreadsheet/ccc?key=0Ajt182RLsguldFlLQmd6Z1ZoczJCenJIdmREUkhxTFE&hl=en_US

注: 主効果のみのモデルの結果は一致しますが、交互作用を含めると一致しません。

前もって感謝します。

4

1 に答える 1

4

彼らは私に同じ結果をもたらします（uclaのデータを使用）。

library(foreign)
d1 <- read.dta('http://www.ats.ucla.edu/stat/stata/webbooks/logistic/apilog.dta')
m1 <- glm(hiqual~cred_hl + pared_hl + cred_hl*pared_hl,family=binomial,  data=d1)
于 2011-08-12T23:52:33.547 に答える