问题标签 [coefficients]

我已将此答案中的代码移植到新的 C++ api 。但是，虽然我了解大部分代码，但我无法理解 calc_shift() 函数背后的想法以及如何提取像素移位。如果有人能给我一个解释，我将不胜感激。功能如下：

上述函数的调用方式如下所示：

我想理解上面的代码，因为我还想将它用于 k < 0 的枕形失真情况。如果我现在使用它，代码将陷入 k<0 给定值的无限循环。

我正在使用标准 ode45 方法求解具有复系数的 Odes 系统，但没有得到理想的结果。下面是我的 M.file

和命令 t=0:500;

[T,Y] = ode45(@tarak,t,[0 0 0 0 0 ]);

结果我得到一些 Y(1) 和 y(2) 的负值。但它应该总是积极的

我的方法有什么问题吗？

谢谢

我正在使用相当大的数据集，并且想估计每个 ID 的两个事件窗口的 F​​ama French 系数。我正在使用以下代码（dummy_reg将观察结果分配给相应的事件窗口）：

此外，我想使用系数来计算事件窗口之间的增量；但是，我不知道如何在估计过程中包含保存/生成新变量。

我尝试了以下方法，但没有奏效：

在 SymPy 0.7.5（Python 2.7.8-64 位）上。当常数项为零或不存在时，我无法从表达式中获取常数项（不存在的常数与零常数相同，对吗？）。

我可以从带有 的表达式中得到常数项expr.coeff(x, 0)。例如：

现在，在这些情况下我怎样才能得到 0？

我敢肯定这已经被问过并记录在案，但至少在没有一些尴尬的解决方法的情况下我找不到答案。谢谢达里奥

编辑交互式会话的完整输出：

MATLAB 表示：“fir2 函数还设计了加窗 FIR 滤波器，但具有任意形状的分段线性频率响应。这与 fir1 形成对比，fir1 仅设计标准低通、高通、带通和带阻配置的滤波器。”

我用 fir1 命令找到了 fir 滤波器系数，并通过使用 freqz 获得频率响应，如下所示

然后我将这个频率响应（hd）传递给 fir2，如下所示

根据 MA​​TLAB，b1 必须等于 b2，但对于 13 阶 FIR 滤波器，
结果如下：

b1=0.0042 0.0063 -0.0000 -0.0403 -0.1221 -0.2103 0.7470 -0.2103 -0.1221 -0.0403 -0.0000 0.0063 0.0042

b2=0.0017 -0.0044 0.0180 -0.0937 0.2075 -0.1097 -0.0012 0.0105 -0.0081 0.0050 -0.0025 0.0010 -0.0005

b1 不等于 b2。理论上这应该是正确的。我不明白出了什么问题。

我使用该包rpart对分类/回归树进行建模。我有变量 x,y,s 其中 x 在 {-1,1} 中，y 在 [0,1] 中是连续的，s$ 是一个具有 3 个级别的因子。

我用

最终模型非常合理，我可以使用fancyRpartPlot(fit). 我也可以修剪树，只剩下一个决定：

它说 $y$ 需要小于某个阈值。如何提取此阈值？我想在数据集上进行引导运行并将阈值保存在数组中。

我正在使用glmnet包（R），我需要提取用于岭回归的帽子矩阵以测试系数的显着性。我查看了函数glmnet的文档和源代码，但我没有找到有关如何获取该矩阵的任何线索。有谁知道是否有可能，使用这个特定的包？

我正在尝试从 qplot 中提取 500 个变量的截距和斜率。我正在使用的代码：

有人可以帮我提取每条回归线（500 条线/变量）的截距和斜率，如附图中所示。

我想从数据框中构建一个系数图。但是使用下面的代码，我遇到了以下错误。

示例代码

这是我试图将任何数学方程变量的系数固定为相同变量的总和系数，如“x^2+x^2+2x-x-25”为“+1x^2+1x^2+2x- 1x-25”，然后求和为“2x^2+x-25”，注意我已经用另一种方法完成了求和过程。

输入案例：

1. 一个
2. X
3. x^2
4. x^2+x^2

输出案例：

1. +1a
2. +1x
3. +1x^2
4. +1x^2+1x^2