실무보기 테이블에 자유도 표시

Question

Sweave 및 R으로 문서를 구성 할 때 테이블에 stargazer 라이브러리를 사용합니다.

stargazer을 사용할 때 glm으로 생성 된 모델의 잔류 편차와 관련된 자유도를 표시하는 메커니즘이 있습니까?

최소 코드 :

library(stargazer) 
set.seed(1234) 
data <- data.frame(x=1:10) 
data$y <- data$x + rnorm(10, 0, 0.2) 
model <- glm(y~x, data=data, family=gaussian) 
summary(model) 
stargazer(model,title="A test", align=T,label="Tab:test",style="all2") 

결과 몽 상가 테이블, 관찰을 가지고 가능성, AIC, 잔류 일탈 및 널 일탈하지만 D.F.를 기록합니다 나는 d.f. 운동 할 수있다. 이것이 직접적으로 표시 될 수 있다고 생각했을 것입니다. 참조 :

https://sites.google.com/site/marekhlavac/stargazer

업데이트 # 1 :

당신의 응답을 당신에게 마렉 감사드립니다. 패키지 내에서 - (I이 다시 올 것이다 4.5되지 않음) http://cran.r-project.org/src/contrib/Archive/stargazer/

에서

1. 버전 4.0을 얻이 발생하는 다른 사람의 이익을 위해, 여기 당신이 주변에 일을 형성 할 수있는 과정이다 R 아래의 디렉토리 구조를 편집하려면 아래의 지침에 따라 "stargazer-internal.R"을 편집하십시오.
2. 라이브러리는 당신의 R 세션
3. 는
4. 는 몽 상가 패키지의 편집 된 버전을 설치에서는 lib 기존 몽 상가를 제거했는지 확인에로드되지 않았 음을 확인합니다.
5. R에서 라이브러리를 다시로드하고 평소와 같이 컴파일하십시오. 명령 줄에서

   detach("package:stargazer", unload=TRUE) 
   remove.packages("stargazer") 
   

   : 여기

는 명령입니다 마지막으로

R CMD INSTALL -l <path to library directory> stargazer 

,

library(stargazer) 
stargazer(model6,model7,model8, title="Logistic model summary",align=T,label="Tab:logmod1", font.size="footnotesize", style="all2") 

결과를 (몇 모델 손이 가정) :

% Table created by stargazer v.4.0 by Marek Hlavac, Harvard University. E-mail: hlavac at fas.harvard.edu 
% Date and time: Tue, Sep 24, 2013 - 17:17:17 
% Requires LaTeX packages: dcolumn 
\begin{table}[!htbp] \centering 
    \caption{Logistic model summary} 
    \label{Tab:logmod1} 
\footnotesize 
\begin{tabular}{@{\extracolsep{5pt}}lD{.}{.}{-3} D{.}{.}{-3} D{.}{.}{-3} } 
\\[-1.8ex]\hline 
\hline \\[-1.8ex] 
& \multicolumn{3}{c}{\textit{Dependent variable:}} \\ 
\cline{2-4} 
\\[-1.8ex] & \multicolumn{3}{c}{whalesighted} \\ 
\\[-1.8ex] & \multicolumn{1}{c}{\textit{logistic}} & \multicolumn{1}{c}{\textit{probit}} & \multicolumn{1}{c}{\textit{glm: binomial}} \\ 
& \multicolumn{1}{c}{\textit{}} & \multicolumn{1}{c}{\textit{}} & \multicolumn{1}{c}{\textit{link = cloglog}} \\ 
\\[-1.8ex] & \multicolumn{1}{c}{(1)} & \multicolumn{1}{c}{(2)} & \multicolumn{1}{c}{(3)}\\ 
\hline \\[-1.8ex] 
visibility & 0.392^{***} & 0.226^{***} & 0.216^{***} \\ 
    & (0.051) & (0.027) & (0.026) \\ 
    Constant & -1.251^{***} & -0.745^{***} & -1.149^{***} \\ 
    & (0.246) & (0.144) & (0.182) \\ 
\hline \\[-1.8ex] 
Observations & \multicolumn{1}{c}{232} & \multicolumn{1}{c}{232} & \multicolumn{1}{c}{232} \\ 
Log Likelihood & \multicolumn{1}{c}{-110.485} & \multicolumn{1}{c}{-110.888} & \multicolumn{1}{c}{-112.694} \\ 
Akaike Inf. Crit. & \multicolumn{1}{c}{224.970} & \multicolumn{1}{c}{225.775} & \multicolumn{1}{c}{229.388} \\ 
Residual Deviance (df = 230) & \multicolumn{1}{c}{220.970} & \multicolumn{1}{c}{221.775} & \multicolumn{1}{c}{225.388} \\ 
Null Deviance (df = 231) & \multicolumn{1}{c}{310.759} & \multicolumn{1}{c}{310.759} & \multicolumn{1}{c}{310.759} \\ 
\hline 
\hline \\[-1.8ex] 
\textit{Note:} & \multicolumn{3}{r}{$^{*}$p$<$0.1; $^{**}$p$<$0.05; $^{***}$p$<$0.01} \\ 
\normalsize 
\end{tabular} 
\end{table} 

오류로 돌아 가기 4.5 코드를 기반으로 해결 방법을 구현할 때 얻을 수 있습니다. 나는 맥 바이너리 (4.5.1 버전) (http://cran.r-project.org/web/packages/stargazer/index.html)에서 설치할 때 실제로는 같은 오류가 발생하고, 단순히 실습을 사용하려고한다. 아래를 보라.

> install.packages("stargazer") 
trying URL 'http://cran.ms.unimelb.edu.au/bin/macosx/contrib/3.0/stargazer_4.5.1.tgz' 
Content type 'application/x-tar' length 332917 bytes (325 Kb) 
opened URL 
================================================== 
downloaded 325 Kb 

> stargazer(model6,model7,model8,  
+   title="Logistic model summary", 
+   align=T, 
+   label="Tab:logmod1", 
+   font.size="footnotesize", 
+   style="all2") 
Error in `rownames<-`(`*tmp*`, value = "visibility") : 
    length of 'dimnames' [1] not equal to array extent 

마렉, 참조 용으로 추적 결과를 이메일로 보내 드리겠습니다. 건배.

Answer 1

여기 시인. 패키지의 기본값은 잔여 및 null 편차에 대한 자유도를 출력하지 않는 것 같습니다. 다음 릴리스에서 기본값을 변경하는 것을 고려할 것입니다.

빠른 해결책으로, (CRAN의) 소스 패키지를 사용하고 stargazer-internal에서 .adjust.settings.style 함수를 수정하는 것이 좋습니다.R는 다음을 포함합니다 :

if (style == "all") { 
     .format.table.parts <<- c("=!","dependent variable label","dependent variables","models","columns","numbers","-","coefficients","-","omit","-","additional","N","R-squared","adjusted R-squared","max R-squared","log likelihood","sigma2","theta(se)*(p)", "SER(df)","F statistic(df)*(p)","chi2(df)*(p)","Wald(df)*(p)","LR(df)*(p)","logrank(df)*(p)","AIC","BIC","UBRE","rho(se)*(p)","Mills(se)*(p)","residual deviance(df)*","null deviance(df)*","=!","notes") 
     .format.coefficient.table.parts <<- c("variable name","coefficient*","standard error","t-stat","p-value") 
    } 

    else if (style == "all2") { 
     .format.table.parts <<- c("=!","dependent variable label","dependent variables","models","columns","numbers","-","coefficients","-","omit","-","additional","N","R-squared","adjusted R-squared","max R-squared","log likelihood","sigma2","theta(se)*(p)", "SER(df)","F statistic(df)*(p)","chi2(df)*(p)","Wald(df)*(p)","LR(df)*(p)","logrank(df)*(p)","AIC","BIC","UBRE","rho(se)*(p)","Mills(se)*(p)","residual deviance(df)*","null deviance(df)*","=!","notes") 
     .format.coefficient.table.parts <<- c("variable name","coefficient*","standard error") 
    } 

참고 여기에 유일한 변화는 내가 "(DF) *"를 "잔류 일탈 '과'널 일탈"에 추가 한 것을.