问题标签 [linearmodels]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
2 回答
32 浏览

r - 如何对某组受访者运行线性模型?

我正在尝试运行线性模型,但仅限于特定的受访者群体。我在网上找到了示例代码,上面说使用 gapminder 功能来隔离受访者组。在这里,我只想在party 列中包含值为 1 的受访者。我按照我在网上找到的示例代码,

但我不断收到此错误

我试着跑步

查看错误发生的位置,但它说的和以前一样。

0 投票
0 回答
123 浏览

python - 为什么 Python 中的 PooledOLS 在长度相同时会给我错误“依赖和 exog 必须具有相同数量的观察值”?

我正在使用线性模型中的面板 OLS 来运行以下回归

但收到一条错误消息,指出内生变量和外生变量的长度不同。

但是,如果我检查两者的形状,y_endog它们X_1是相同的形状。y_endog我还通过在创建和数据框之前删除任何行并重置索引来确保没有 NaN 值X_1

X_1 = 1如果我运行估计中列数有效的代码。然而,一旦X_1超过 1 列,估计就不再起作用(例如,向 PooledOLS 模型添加一个附加变量)。

数据框

y_endog

y_endog

X_1

X_1

0 投票
0 回答
14 浏览

linearmodels - 在线性模型中,*.f_statistic 零假设真的是“所有参数前常数不为零”吗?

https://bashtage.github.io/linearmodels/panel/panel/linearmodels.panel.results.RandomEffectsResults.f_statistic.html?highlight=f%20statistic#linearmodels.panel.results.RandomEffectsResults.f_statistic上的文档报告 *. f_statistic 是经典的“经典 F 统计量。非常量回归量的显着性联合检验”。在这种情况下,我阅读并且空值是没有自变量的模型适合数据以及模型。


但是,当调用代码的 null 元素的实际字符串值时,它会显示“H0:所有参数,例如常量不为零”……恰恰相反。

null 只是错误地记录在代码中吗?这是linearmodels 4.17 btw ...谢谢!

结果的屏幕截图。

在此处输入图像描述

0 投票
0 回答
26 浏览

python - 我对 python 很陌生,我想在 sns.lmplot 中增加我的图例的字体大小



###用于轴



###for 轴刻度线标签



###为了传说,

传奇已经出现在情节之外,但我想更改文本的字体大小(基本上是 5 个分类单元名称)和标题

预先感谢您的帮助。巴纳马利

0 投票
1 回答
1024 浏览

python - 使用 Python 进行具有固定效应的面板数据回归

我将以下面板存储在df

状态 是的 持续的 x1 x2 时间
0 01 01001 2009 12 1 0.956007 639673 1
1 01 01001 2010 20 1 0.972175 639673 2
2 01 01001 2011 22 1 0.988343 639673 3
3 01 01002 2009 0 1 0 33746 1
4 01 01002 2010 1 1 0.225071 33746 2
5 01 01002 2011 5 1 0.450142 33746 3
6 01 01003 2009 0 1 0 45196 1
7 01 01003 2010 5 1 0.427477 45196 2
8 01 01003 2011 9 1 0.854955 45196 3
  • y是每个地区的抗议数量
  • constant是一列充满的
  • x1是移动网络提供商覆盖的地区区域的比例
  • x2是各区的人口数(注意是时间固定的)

如何在 Python 中运行以下模型? 模型


这是我尝试过的

ValueError: exog 没有完整的列排名。如果您希望继续进行模型估计而不考虑系数估计的数值精度,您可以设置 rank_check=False。

我究竟做错了什么?

0 投票
0 回答
20 浏览

r - DESeq2中设计模型的差异及其解释

我正在使用 DEseq2 并试图了解使用不同模型获得的结果。我的数据设计为

我有 3 个模型。

模型1:~基因型

模型2:~时间

模型3:~基因型+时间

我使用模型 1 得到 200 个显着基因,从模型 2 中得到另外 198 个显着基因。但是,当我使用 model3 时,我没有任何重要的基因。我想知道使用模型 3 时正在比较哪些样本。

即,当使用 (obj3, name="genotype_KO_vs_WT") 获得结果时,在 model3 中比较了哪些样本?它与model1有什么不同?

另外,使用 results(obj2, name="time_T2_vs_T1") 获取结果时,model3 中正在比较哪些样本?和model2有什么区别?

0 投票
1 回答
71 浏览

python - 如何在 Statsmodels 或 Linearmodels 中测试 2SLS 的系数相等性?

因此,如果我对多个治疗组和一个对照组进行实验,我会使用 Statsmodel ols 分析结果,看看是否有任何治疗组与对照组有统计学差异:

y ~ C(treatment_group, Treatment('控制')

然后我会运行 results.t_test_pairwise() 来确定每个治疗组的系数是否相等。即了解每个治疗组的结果是否在统计学上显着不同。

在当前情况下,我不只是运行标准 ols,而是使用 Statsmodel/Linearmodel 的 2SLS 来分析工具变量。我可以很好地运行分析,并得到结果。但是现在我需要看看不同仪器(三个不同治疗组)的系数是否相同,所以我知道不同治疗组的效果是否不同。

statsmodel的代码:

或者对于线性模型:

Josef's response here建议您可以使用 wald t 检验,但我需要使用限制矩阵而不是公式。因此,如果有人对如何做到这一点有任何想法,那将不胜感激。

0 投票
0 回答
42 浏览

python - 具有固定效应的两阶段回归最小二乘回归

是否有一种内置方法可以在 Python 的两阶段最小二乘回归中使用固定效果?

linearmodels有一个用于固定效果的专用模块和一个用于两阶段最小二乘的单独模块。前者允许您传入列以用作固定效果(每个级别一个虚拟列)。但是,我似乎在后者中找不到这个选项。

使用PanelOLS,您可以为回归添加固定效果。

我想知道这是否可以通过IV2SLS模块完成。

对于没有使示例可重现,我提前道歉。

0 投票
1 回答
71 浏览

python - 将线性模型回归摘要输出为乳胶

如何将拟合的线性模型对象的汇总表打印为乳胶?

例如,如何打印res为乳胶代码?

0 投票
2 回答
22 浏览

python - 使用线性模型将摘要另存为图像

我正在尝试将我正在使用包 linearmodels 进行的一些回归的摘要保存为 PNG 文件

我试图这样做,但它不起作用:

你能帮我找出问题所在吗?提前致谢!