Używałem Stata do uruchamiania negatywnej regresji dwumianowej w replikacji. Nie wiesz, co kryje się pod tym, w jaki sposób Stata to robi, ale chciał wiedzieć, czy istnieje funkcja/pakiet R, który robi to samo? R da mi lepsze pojęcie o tym, jak to działa, ponieważ widzę kod.jest funkcja R dla xtnbreg Stata?
Zajrzyj do funkcji glm.nb w pakiecie MASS. Jeśli interesuje Cię to, co dzieje się "pod maską", możesz zobaczyć kod źródłowy tej funkcji, wpisując jej nazwę w wierszu polecenia.
Jeśli wygodniej jest używać R, prawdopodobnie jest to sposób na wykonanie; Jednakże, jeśli jesteś zainteresowany tym, co znajduje się „pod maską” w Stata, zawsze można zobaczyć, co dzieje się w taki sam sposób jak w R stosując
set trace on
aby zobaczyć, co kod jest uruchomiony (lub do pobrania tr z SSC) lub przy użyciu
viewsource xtnbreg.ado
aby zobaczyć rzeczywisty kod, który jest prowadzony przez xtnbreg.
Jeśli interesuje Cię sposób, w jaki Stata oblicza wyniki w xtnbreg, jest szczegółowa dyskusja na temat funkcji wiarygodności w podręczniku [XT] strona 367-370 wraz z odnośnikami.
Oto link do przykładu glm.nb: http://www.ats.ucla.edu/stat/r/dae/nbreg.htm –