StackOverflow2013

Note that there are some explanatory texts on larger screens.

plurals

POStata and R results not matched in Logistic Regression with two categorical predictors and their interaction
text
Body
copied!<p>I am getting confused when i am trying to compare the results of Stata and R. I am using example given on the webpage <a href="http://www.ats.ucla.edu/stat/stata/webbooks/logistic/chapter2/default.htm" rel="nofollow">http://www.ats.ucla.edu/stat/stata/webbooks/logistic/chapter2/default.htm</a> First run the following command in Stata </p> <pre><code>use http://www.ats.ucla.edu/stat/stata/webbooks/logistic/apilog, clear </code></pre> <p>and then use following commands given in the section (2.2.2 A 2 by 2 Layout with Main Effects and Interaction)</p> <pre><code>generate cred_ed = cred_hl*pared_hl logit hiqual cred_hl pared_hl cred_ed </code></pre> <p>These two command will produce the results given on the webpage.</p> <p>And then i have used following R code to reproduce same example</p> <pre><code>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) </code></pre> <p>But results are not matched!</p> <p>Data file for R can be download from following link</p> <p><a href="https://spreadsheets.google.com/spreadsheet/ccc?key=0Ajt182RLsguldFlLQmd6Z1ZoczJCenJIdmREUkhxTFE&hl=en_US" rel="nofollow">https://spreadsheets.google.com/spreadsheet/ccc?key=0Ajt182RLsguldFlLQmd6Z1ZoczJCenJIdmREUkhxTFE&hl=en_US</a></p> <p>Note: Results for model with only main effects are matched but when we include interaction, it is not matched. </p> <p>Thanks in Advance.</p>

Querying!

Guidance

An individual column

Larger individual text columns get their own page to allow for proper reading.

SQuiL has stopped working due to an internal error.

If you are curious you may find further information in the browser console, which is accessible through the devtools (F12).

Reload