1

我想问如何计算inf。包中的固定效应 logit 模型的标准，例如 AIC 等 bife。

基本summmary输出不包括 AIC，但是在查看时：Goodness-of-fit for fixed effect logit model using 'bife' package

计算了 AIC 标准。我怎么在我的摘要输出和对数似然中都没有它。

dta = bife::psid
mod_logit <- bife(LFP ~ AGE + I(INCH / 1000) + KID1 + KID2 + KID3 | ID, data = dta, bias_corr = "ana")
summary(mod_logit)
4

1 回答 1

1

如果您检查bife代码，AIC 至少在 0.5 版本中是在早期版本中计算的。您可能正在使用不再包含 AIC 的当前版本 0.6。

如果您不介意使用旧版本，请尝试以下操作：

  1. 从您的库中删除当前版本。

  2. 从 CRAN 网站下载 0.5 版：https ://cran.r-project.org/src/contrib/Archive/bife/

  3. 安装到您的计算机：install.packages("D:\\bife_0.5.tar.gz", repos = NULL, type="source").假设它存储在D:驱动器上。

或者：

require(devtools)
install_version("bife", version = "0.5", repos = "http://cran.us.r-project.org")

如果安装成功，请在包含 AIC 的情况下运行以下命令：

library(bife)
dta = bife::psid
mod_logit <- bife(LFP ~ AGE + I(INCH / 1000) + KID1 + KID2 + KID3 | ID, data = dta, bias_corr = "ana")
summary(mod_logit)
#> ---------------------------------------------------------------
#> Fixed effects logit model
#> with analytical bias-correction
#> 
#> Estimated model:
#> LFP ~ AGE + I(INCH/1000) + KID1 + KID2 + KID3 | ID
#> 
#> Log-Likelihood= -3045.505 
#> n= 13149, number of events= 9516
#> Demeaning converged after 5 iteration(s)
#> Offset converged after 3 iteration(s)
#> 
#> Corrected structural parameter(s):
#> 
#>               Estimate Std. error t-value  Pr(> t)    
#> AGE           0.033945   0.012990   2.613  0.00898 ** 
#> I(INCH/1000) -0.007630   0.001993  -3.829  0.00013 ***
#> KID1         -1.052985   0.096186 -10.947  < 2e-16 ***
#> KID2         -0.509178   0.084510  -6.025 1.74e-09 ***
#> KID3         -0.010562   0.060413  -0.175  0.86121    
#> ---
#> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#> 
#> AIC=  9023.011 , BIC=  19994.7 
#> 
#> 
#> Average individual fixed effects= 0.0122
#> ---------------------------------------------------------------

由reprex 包于 2020-01-09 创建(v0.3.0)

于 2020-01-09T06:23:24.797 回答