Hello,
I am learning R and I am fluent in Stata and I try to translate part of my
Stata code to R to check the reliability of the data under R. I have a
proportion variable as a dependent variable pQSfteHT . Independent variables
are dummies for two categorical variables called dQSvacrateHTQuali3 and
cluster_3. I am fitting a model with the Stata command below:
glm pQSfteHT dQSvacrateHTQuali3_2 dQSvacrateHTQuali3_3 dQSvacrateHTQuali3_4
dQSvacrateHTQuali3_5 cluster_32 cluster_33 cluster_34, link(probit)
family(binomial) robust
and the same (I expect) model with R with the command below:
nurse.model<-glm(pQSfteHT~dQSvacrateHTQuali3_2 + dQSvacrateHTQuali3_3 +
dQSvacrateHTQuali3_4 + dQSvacrateHTQuali3_5 + cluster_32 + cluster_33 +
cluster_34 ,family=binomial(link = "logit"))
I found some differences in the parameters, could it come from the "robust"
option in the Stata command? It sounds strange that a variance option would
lead to changes in parameters estimation but I am not an econometrician.
Is anyone bilingual in R and Stata and could have a look at the syntaxes
above?
Thank you in advance
Thank you also to the people answering my previous enquiry.
Jean-Baptiste
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
