我正在用R编写一个林肯-彼得森MLE函数。我有两个版本的分布函数，一个使用gamma函数，另一个使用R内部dmultinom (使用lgamma)。当我插入样本值时，两者返回相同的结果，但是dmultinom版本无法提供可接受的MLE估计。我很好奇这是什么原因造成的。(18,0,0), method="L-BFGS-B",
control=list(parsca
浏览 5提问于2019-06-20得票数 1
回答已采纳
我正在尝试在模拟的i.i.d.上实现Fisher评分。Poisson数据，但我收到堆栈溢出错误。我对函数做了一些简单的打印，发现在第一次迭代之后，猜测值没有变化。
浏览 12提问于2020-03-03得票数 4
回答已采纳
实现这个简单的寻根算法。我一辈子都找不出我的实现出了什么问题。树根不断膨胀，没有收敛的迹象。有什么建议吗？{ q0 = q; s0 = s;
q = q0 - f(q0)/((q
浏览 3提问于2013-11-10得票数 0
回答已采纳
我正在尝试用optim()实现混合高斯混合的MLE，使用R的本地数据集(来自质量的Geyser)。我的密码在下面。问题是optim工作得很好，但是，返回我传递给它的原始参数，并且还说它已经收敛。如果你能指出我偏离轨道的地方，我将不胜感激。我的期望是，它至少会产生不同的结果，即使不是大不相同。library(ggplot2)data("geyser")
浏览 0提问于2015-08-19得票数 0
回答已采纳
1回答
我正在使用mle2来估计一个非线性模型的参数，我想要估计参数估计(std.错误)。同样，我想使用模型来预测新的数据，我在这个过程中的一些步骤中遇到了问题(错误)。我尝试了许多不同的优化器和方法，但都得到了无法计算std的相同错误。错误，通常伴随着无法反转hessian的警告。或者我对我的参数的估计非常不可靠，没
浏览 31提问于2019-10-31得票数 0
1回答
我希望您不要认为这个问题没有意义，因为我一直在努力处理代码，我也不明白为什么它不收敛到相同的值。我编写了Newtow码来估计一个多项式logistic模型，当它收敛时，我得到了与我的对数和R函数多项式中的构造不同的显着值。0.7070484 -0.4951407 0.1878713
2 1.845272 -0.26688254 0.6945098 -0.2749830 -2.03
浏览 2提问于2017-10-06得票数 0
回答已采纳
1回答
我在玩Wolfram元胞自动机的规则110。给定一行0和1，您可以使用以下规则计算下一行：从00000000开始.1最后得到以下序列：出于好奇，我决定用多项式近似这些规则，这样细胞不仅可以是0和1，还可以是灰色之间的颜色：我的问题是为什么会发生这种情况，以及我如何解决这个问题。尝试了不同的除草优化器，但它们的性能都比LBFGS差。尝试了不同的学习速度。有什么想法吗?为什么会发
浏览 0提问于2018-05-31得票数 1
我正在使用R中的optim运行一个优化例程，并告诉程序我想要返回什么。例如，如果我放入return(op1$par)，它将返回我的所有4个变量值。这很好，如果我运行return(op1)，我显然会从优化例程中获得所有信息(par、值、收敛性等)。然而，在这种格式中，无法在输出中访问par值，它只是简单地详细说明了有4个值。现在我需要的是同时获得参数值和收敛
浏览 0提问于2012-09-08得票数 1
回答已采纳
1回答
C语言中的二分法
、、、、
我写了一段代码，用二分法求一个四次多项式的根。我也为第三个多项式写了同样的代码，运行得很好。我只是复制和粘贴，并添加第四度术语，它不工作fine.Here是我的代码它给了我当我将间隔改为-50，51时，输出变为0.500000，我认为问题是关于边界的</em
浏览 30提问于2021-03-12得票数 0
1回答
我在Python中实现了梯度下降，使用MSE作为损失函数执行正则化多项式回归，但在线性数据上(以证明正则化的作用)。所以我的模型在表格下面：在我的损失函数中，R代表正则化项：让我们以L2-范数作为正则化，损失函数w.r.t的偏导数。妇女的情况如下：最后，使用恒定的学习速率更新了相关系数：问题是，我不
浏览 2提问于2019-12-27得票数 1
回答已采纳
我使用R {fExtremes}为我的数据(向量)寻找GEV分布的最佳参数。但是获得以下错误消息
我追溯到拟合$ hessian，发现我的gevFit()的源代码(gevFit())显示fit$hessian是由optim()计算的。输出参数与初始参数<em
浏览 0提问于2016-08-23得票数 0
回答已采纳
当使用R-package glmmTMB() (参见)时，我知道在处理模型的收敛性时，我有一些选择。更具体地说，有一个control =参数可以传递glmmTMBControl()参数，手册中的部分如下：此外，其中一个小片段--即 --明确地谈到了如何处理收敛问题。=1e3))，即增加迭代次数
glmmTMBControl(optimizer=optim, optArgs=list(method="BFGS
浏览 19提问于2022-07-27得票数 2
回答已采纳
2回答
问题的简要说明:我使用牛顿-拉夫森算法在多项式中求根，但在某些情况下不起作用。为什么？我从“c++中的数值配方”中提取了一个牛顿-拉夫森混合算法，如果新-拉斐尔算法不能正确收敛(导数值低或收敛速度不快)，它就会一分为二。
我用几个多项式检查了算法，它起作用了。现在我在我的软件中测试，我总是得到一个特定的多项式的
浏览 4提问于2012-11-13得票数 9
回答已采纳
2回答
在课程中，没有关于时代的东西，但在实践中，它们随处可见。如果优化器在一次传递中找到了最佳权重，我们为什么还需要它们。为什么模型会有所改进？
浏览 1提问于2017-03-10得票数 4
我想要对R中的伽马分布的最大似然进行编程；到目前为止，我已经完成了以下工作：x<-scan("http://www.cmc.edu/pages/faculty/MONeillloglike2<-function(LL){beta<-LL$b
mle(
浏览 0提问于2015-02-22得票数 1
2回答
我要攻击的函数如下所示：在此可视化：Eq1用弧度表示，Eq2我用度表示如果有人能给我指明正确的方向，我将不胜感激。
浏览 3提问于2016-04-05得票数 0
我正在训练一个非常深的神经网络( Panoptic )，它在Google上的ResNet34主干上用于CityScapes数据集的全景分割，并注意到，由于作物大小很大，批处理的大小必须减少到每批1张图像，否则就会出现内存不足的问题。虽然我知道这会在训练中造成偏斜，而且很可能很难达到很好的收敛性，但我能否向专家们一般地问这个问题:在基于图像的处理中，一般认为1的批次
浏览 0提问于2022-06-30得票数 2
回答已采纳
我使用下面的调优代码来为我的svn模型找到最好的情况和epsilon。., ranges = list(epsilon = seq(-5.0, 5, 0.1), cost = 2^(0:3)))
但令人惊讶的是，它暗示了不幸的是，模型不像没有这些参数的模型那样精确。
浏览 0提问于2018-08-15得票数 1
回答已采纳
在下面的链接中，解释了如何使用元组返回多个值：Returning multiple values from a C++ function 所以我试着按照指示进行编码，以返回商和多项式除法的余数，但它不起作用; delta = R.degre() - D.degre(); return {Q, R};
}然后，当我按照上面的链接开始编写
浏览 3提问于2021-09-02得票数 0
回答已采纳
1回答
我正试图在我的ubuntu机器上执行一次繁重的科学计算。代码是用Python编写的，并被标记为可执行文件。然而，大约一个小时后，计算被终止，除了" terminated“之外没有任何错误消息，就像系统刚刚运行命令"killall python”一样。谢谢!
浏览 3提问于2012-04-10得票数 0
回答已采纳
云服务器
ICP备案
对象存储
云直播
腾讯会议
腾讯云开发者
