问题标签 [gauss]

我得到了一个用 GAUSS 编写的 MLE 估计器，我正试图将它重新编码到 R 中。我不使用，也从未使用过 GAUSS 本身（也无权访问它）。在代码中，有一行让我有点困惑。

在经过验证的 GAUSS 代码中，“输入”（部分注释）下有一行说明：

我正在努力让代码一块一块地工作，但我的第一个问题是相对简单的。

这是让我感到困惑的GAUSS片段：

我知道这是单位矩阵乘以“输入”invsig，但是在有效的模拟示例中（从附加到代码的日志文件中），可以使用 invsig 的标量值启动程序。IE：将初始值设置为 invsig = 1

在 R 中，这是行不通的。这是尝试获取此信息的简单“测试”代码：

我了解我在 R 中遇到的不符合错误。问题是我在 GAUSS 和 R 之间的转换中遗漏了什么？在阅读在线帮助手册时，GAUSS 通过使用单个符号 (*/+-) 进行矩阵运算，并且要按元素进行操作，请添加“。” 每次操作前。所以对我来说，GAUSS 代码说要进行矩阵乘法（%*%在 R 中），这就是简单函数试图做的事情。

非常感谢任何意见或建议！

我想使用在高斯统计软件中实现的面板数据计量经济学检验。使用 R 包MASS::write.matrix，我设法生成了一个 ASCII 文件并从 Gauss 中读取该文件。这适用于 txn 矩阵。但我想知道如何导出 atx nk 矩阵。nk 列会简单地相互附加吗？

我是这里的新手，也是 python 编程的新手，所以我需要一些帮助。我需要为教师制作关于 K-means 聚类和二维高斯分布的项目。我已经对三个集群进行了完整的编程，但我不知道必须为这三个集群制作 2D Gaussian。

这是我用来展示 K-means 算法的算法[我来自克罗地亚，所以我的语言中有一些词]：

我需要在这个算法之后制作二维高斯分布。你们能帮帮我吗？

此代码用于与矩阵一起使用的高斯消除，我不明白为什么屏幕上会出现以下错误。

错误是：

该程序应考虑以下功能： 将要求解的系统的二维数组作为现有变量存储器传递给程序。- 数据是双浮点类型。- 程序应将结果保存在另一个变量中。-程序应该在被解决之前显示线性方程组，最后应该显示解决方案。

我做了这个显示矩阵的代码：

我不知道如何应用高斯方法。你们能帮帮我吗？

在 Gauss-legendre 积分中，我们需要找到 legendre 函数的零点，但我找不到编写代码的方法来让我做到这一点？我知道那里有一个“xi”列表，这个函数等于零，但是我们可以编写程序自己找到这些“xi”吗？

有没有通用的解决方案？您必须将它们相乘，但很难实现。

对于二维情况，您可以使用表示单个正态分布的两个向量的外积。

我正在尝试编写一个通过高斯和其他一些过滤器过滤位图的代码。我开发了一个代码，它根据内核大小和标准偏差等输入参数生成内核。我还编写了一个解析位图图像并将其放入指针的代码。我的问题是，如何将生成的内核矩阵应用于指针，或者是否有其他方法可以将高斯内核值应用于该位图指针。

这里的示例是我使用不同过滤器（darkenImage）的代码，以便您可以看到我的代码样式：

这里是主函数的输入、输出部分

那么，如何制作 GaussFilter 函数呢？

我想知道高斯拉盖尔如何在大范围内工作。例如，我有一个从 (0, +inf) 开始的二维函数。当我在 python 中使用 gauss laguerre 通过对权重和横坐标进行汇总来对函数进行采样时，我并没有得到接近我使用的东西，比如 dblquad。下面是集成的示例代码。lgw 输出权重和横坐标，然后通过使用两个 for 循环将其用于双重积分。我没有看到像 x, y = 1e8, 1e8 这样的样本点是如何被捕获的。增加 n 不会给出高横坐标（至少不是要求很高）。

有人可以解释如何捕获更高的样本点吗？我没有正确使用正交吗？我可以集成具有小限制的功能，例如 1e2 左右。如果限制很高，比如 1e15，该怎么办？我从理论中看到了定义，但我看不到如何捕获更高的权重和横坐标。

谢谢

编辑：不可能进一步减少我的功能。被积函数的不同部分是用数值计算的，所以我没有任何解析表达式。我只能说函数是平滑的并且具有正弦曲线的行为。

“求 1000 以下的 3 或 5 的所有倍数之和”

我无法理解为什么下面的解决方案仍然返回正确的结果，因为 x3、x5 和 x15 在除法后使用 int。这意味着除法的结果总是四舍五入，小数点被忽略。

当我尝试用双打替换所有 3 个整数时，我得到了错误的结果。

该解决方案基于以下观察：

1 + 2 + ... + n = n*(n+1)/2

参考