2010-08-18 12 views
6

Si sta utilizzando Stata per eseguire regressioni binomiali negative in una replica. Non sei sicuro di cosa sia sotto il cofano su come Stata fa questo, ma vuoi sapere se c'è una funzione/pacchetto R che fa la stessa cosa? La R mi darà un'idea migliore di come funziona, dal momento che posso vedere il codice.c'è una funzione R per xtnbreg di Stata?

risposta

4

Controllare la funzione glm.nb nel pacchetto MASS. Se sei interessato a quello che sta succedendo "sotto il cofano", puoi vedere il codice sorgente della funzione semplicemente inserendo il suo nome al prompt dei comandi.

+1

ecco un collegamento a un esempio glm.nb: http://www.ats.ucla.edu/stat/r/dae/nbreg.htm –

1

Se sei più a tuo agio nell'usare R, allora questa è probabilmente la strada da percorrere; tuttavia, se siete interessati a ciò che è "sotto il cofano" in Stata, si può sempre vedere cosa sta succedendo più o meno allo stesso modo come in R utilizzando

set trace on 

per vedere quale codice è in esecuzione (o da scaricare tr di SSC) o utilizzando

viewsource xtnbreg.ado 

per vedere il codice effettivo che è gestito da xtnbreg.

Se sei interessato a come Stata sta calcolando i risultati in xtnbreg, c'è una discussione dettagliata della funzione di verosimiglianza nella pagina di manuale [XT] 367-370 con riferimenti inclusi.

Problemi correlati