[R] Sata and R users GLM methods translation

Jean-Baptiste Combes
2010-01-22T15:32:53Z

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