需要澄清Coef。以及圣厄尔. 混合lm结果中的随机参数

2024-05-26 20:48:05 发布

您现在位置:Python中文网/ 问答频道 /正文

我试图理解pythonstatsmodel包提供的混合线性模型的结果。我想避免数据分析和解释中的陷阱。问题在数据加载/输出代码块之后。在

加载数据和拟合模型:

import statsmodels.api as sm
import statsmodels.formula.api as smf
data = sm.datasets.get_rdataset("dietox", "geepack").data
md = smf.mixedlm("Weight ~ Time", data, groups=data["Pig"])
mdf = md.fit()
print mdf.summary()

Mixed Linear Model Regression Results
========================================================
Model:            MixedLM Dependent Variable: Weight    
No. Observations: 861     Method:             REML      
No. Groups:       72      Scale:              11.3669   
Min. group size:  11      Likelihood:         -2404.7753
Max. group size:  12      Converged:          Yes       
Mean group size:  12.0                                  
--------------------------------------------------------
             Coef.  Std.Err.    z    P>|z| [0.025 0.975]
--------------------------------------------------------
Intercept    15.724    0.788  19.952 0.000 14.179 17.268
Time          6.943    0.033 207.939 0.000  6.877  7.008
Group Var    40.394    2.149                            
========================================================

问题1。(a) 什么是群体Var系数(params)?我以为这是组变量(cov_params)的方差,但默认输出与内置方法输出不匹配。在

问题1。(b) “组变量”参数(params)是什么意思?在

^{pr2}$

问题2。(a) 组变量的标准误差(bse)是什么意思?为什么在默认输出中没有报告组Var估计值?这不重要吗?在

问题2。(b) 与标准方差误差(bse)有何不同?在

print "-----Standard errors-----"
print mdf.bse
print
print "-----Standard errors of random effects-----"
print mdf.bse_re

-----Standard errors-----
Intercept    0.788053
Time         0.033387
Group Var    0.637336
dtype: float64

-----Standard errors of random effects-----
Group Var    2.148771
dtype: float64

问题3。为什么summary()中没有报告随机参数的t值和p值?在

print "-----t-values (or z-values?)-----"
print mdf.tvalues
print
print "-----p-values-----"
print mdf.pvalues

-----t-values (or z-values?)-----
Intercept     19.952366
Time         207.938608
Group Var      5.575760
dtype: float64

-----p-values-----
Intercept    1.429597e-88
Time         0.000000e+00
Group Var    2.464519e-08
dtype: float64

参考号:https://www.statsmodels.org/dev/mixed_linear.html


Tags: datatimevargroupstandardbsemdfvalues

热门问题