医学性别、年龄、受教育年限的校正方法与代码

条件漫步

已于 2023-03-19 14:29:08 修改

阅读量776

点赞数

分类专栏：统计知识文章标签： python 机器学习医学统计检验

于 2022-07-17 13:42:47 首次发布

本文链接：https://blog.csdn.net/chenhepg/article/details/125831044

版权

统计知识专栏收录该内容

9 篇文章 0 订阅

订阅专栏

@创建于：2022.07.17
@修改于：2022.07.17

文章目录

1、参数介绍

statsmodels.regression.linear_model.OLS(endog, exog=None, missing='none', hasconst=None, **kwargs)


Parameters

    endog: array_like
        A 1-d endogenous response variable. The dependent variable.
        
    exog: array_like
        A nobs x k array where nobs is the number of observations and k is the number of regressors. An intercept is not included by default and should be added by the user. See statsmodels.tools.add_constant.

    missing: str
        Available options are ‘none’, ‘drop’, and ‘raise’. If ‘none’, no nan checking is done. If ‘drop’, any observations with nans are dropped. If ‘raise’, an error is raised. Default is ‘none’.
    
    hasconst: None or bool
        Indicates whether the RHS includes a user-supplied constant. If True, a constant is not checked for and k_constant is set to 1 and all result statistics are calculated as if a constant is present. If False, a constant is not checked for and k_constant is set to 0.

指示 RHS 是否包含用户提供的常量。
如果为 True，则不检查常量并将 k_constant 设置为 1，并且所有结果统计信息都按照存在常量的方式进行计算。
如果为 False，则不检查常量并将 k_constant 设置为 0。

statsmodels.regression.linear_model.OLS

2、官网代码

>>> import statsmodels.api as sm
>>> import numpy as np
>>> duncan_prestige = sm.datasets.get_rdataset("Duncan", "carData")
>>> Y = duncan_prestige.data['income']
>>> X = duncan_prestige.data['education']
>>> X = sm.add_constant(X)
>>> model = sm.OLS(Y,X)
>>> results = model.fit()
>>> results.params
const        10.603498
education     0.594859
dtype: float64

>>> results.tvalues
const        2.039813
education    6.892802
dtype: float64

>>> print(results.t_test([1, 0]))
                             Test for Constraints
==============================================================================
                 coef    std err          t      P>|t|      [0.025      0.975]
------------------------------------------------------------------------------
c0            10.6035      5.198      2.040      0.048       0.120      21.087
==============================================================================

>>> print(results.f_test(np.identity(2)))
<F test: F=array([[159.63031026]]), p=1.2607168903696672e-20, df_denom=43, df_num=2>