2013-04-26 12 views
11

Podczas korzystania z aplikacji Stargazer do tworzenia tabeli LaTeX na regresie logistycznym, standardowym zachowaniem jest wyprowadzanie wartości logitowych każdego modelu. Czy można zamiast tego uzyskać exp (logit)? Czy mogę zamiast tego uzyskać wskaźniki Odds?Ilości szans zamiast logit w stargazer() Wynik LaTeXa

W dokumentacji Stargazer następujący wzmianka o "Coef" -argument, ale nie rozumiem, czy to może włączyć exp (logits).

Coef: lista wektorów liczbowych, które zastąpią domyślne Coef wartości fi na po- dla każdego modelu. Nazwy elementów będą używane do dopasowywania współczynników do poszczególnych zmiennych towarzyszących i dlatego powinny być zgodne z nazwami współzmiennymi . Wektor NULL wskazuje, że dla danego modelu należy użyć domyślnego zestawu współczynników . Natomiast wektor NA oznacza, że ​​wszystkie współczynniki modelu powinny pozostać puste.

Odpowiedz

8

Zgodnie z komentarzem symbiotycznym w 2014 r. Nowsze wersje "gwiazdora" mają zastosowanie "". "" Dla "" coef "" se "" "t" "" p " "i" ci ", umożliwiając bezpośrednie przekształcenie tych statystyk.

apply.coef a function that will be applied to the coefficients. 
apply.se a function that will be applied to the standard errors. 
apply.t a function that will be applied to the test statistics. 
apply.p a function that will be applied to the p-values. 
apply.ci a function that will be applied to the lower and upper bounds of the confidence intervals. 

Znaczenie można bezpośrednio wykorzystać ...

stargazer(model, 
      apply.coef = exp, 
      apply.se = exp) 

EDIT: Zauważyłem jednak, że po prostu exponentiating cis robi nie dać tego, czego można się spodziewać.

EDYCJA: Możesz uzyskać prawidłowe CI za pomocą metody opisanej pod numerem here.

+2

pamiętać, że standardowy błąd ilorazu szans nie jest wykładnikiem s.e. z coeff. Raczej jest to se (OR) = exp (coeff) * se (coeff). Do zobaczenia, np. http://www.stata.com/statalist/archive/2005-09/msg00829.html – dzeltzer

7

stargazer pozwala na zastąpienie wielu rzeczy, zależnych zmiennych etykiet, etykiet zmiennych dodatkowych i tak dalej. Aby zastąpić te, które musisz dostarczyć wektor etykiet zmiennych, robi się to, aby domyślnie mieć publicznie dostępne nazwy wierszy zamiast nazw zmiennych od R.

Aby uzyskać współczynniki szans, musisz dostarczyć wektor ilorazów szans do stargazer. Jak zdobyć ten wektor? Bardzo łatwo, faktycznie. Załóżmy, że dany model jest nazywany model, wówczas kod jest:

coef.vector <- exp(model$coef) 
stargazer(model,coef=list(coef.vector)) 

Jeśli masz kilka modeli w tabeli, a następnie lista powinna zostać rozszerzona, na przykład coef=list(coef.vector1,coef.vector2,...), gdzie wszystkie wektory na liście będą pochodzić z podobnej potęgi jak powyżej.

+1

Świetnie, to działa dobrze dla współczynników, ale jak uzyskać prawidłowe błędy standardowe i Con przedziały ufności. – MYaseen208

+0

Aby uzyskać prawidłowe błędy standardowe, można użyć argumentu apply.se apply.se, aby uzyskać wszystkie potknięcia standardowe, które należy potęgować. Podobnie można użyć apply.ci dla przedziałów ufności – symbiotic

3

Problem polega na tym, że chcesz wyświetlić wskaźnik prawdopodobieństwa (bez logu), ale zachowaj statystyki testu na podstawie podstawowego modelu liniowego. Domyślnie, gdy używasz jednej z metod "apply", takich jak apply.coef = exp, stargazer ponownie obliczy statystyki t i p. Nie chcemy tego. Również standardowe błędy są w logu, ale nie możemy ich po prostu potęgować.Moje preferowane podejście jest:

  1. exponentiate z coefs w Stargazer
  2. wyłączyć auto P i auto t
  3. raport (nietransformowane) t-statystyka w tabeli zamiast standardowych błędów

W kodzie jest to:

stargazer(model, apply.coef=exp, t.auto=F, p.auto=F, report = "vct*") 
Powiązane problemy