stata - Fractional logit models in SAS -


institutionally constrained using sas (yes, know). have basic specification run in stata/r no problem: fractional logit model (papke wooldridge 1996). it's glm binomial distribution assumption , logit link function. data context stationary time series in unit interval—percentage data.

in stata run

glm y x, family(binomial) link(logit) 

in r

amodel <- glm(y ~ x, family=binomial(link=logit), data = adataframe) 

attempting in sas using proc glimmix:

proc glimmix data =adatatable method = rspl; class someclassifier anotherclassifier; model y = x / dist = binomial link = logit solution;  random _residual_;  run; 

i'm dealing panel dataset, doesn't matter in r or stata syntax appears needed information proc glimmix, hence inclusion of 'class' line. able fit models close original stata/r differ in non trivial ways when @ individual parameters or predicted values (correlation between different predicted values .97). can advise on proper way fractional logit in sas? think inclusion of "random" line have above 1 source of trouble, seems add random effects model via matrix * vector operation.

solution simple turns out. need use: method = quad use quasi maximum likelihood estimation, same used in stata , r.


Comments

Popular posts from this blog

neo4j - finding mutual friends in a cypher statement starting with three or more persons -

php - How to remove letter in front of the word laravel -

minify - Minimizing css files -