在R中使用optim和L-BFGS-B方法时的非有限差分误差
是指在优化问题中使用这两种方法时，计算梯度的近似值与真实梯度之间的误差。
优化问题是云计算领域中的一种常见任务，用于寻找某个函数的最优解。在R中，optim函数提供了多种优化算法，其中包括L-BFGS-B方法，它是一种拟牛顿方法，用于求解无约束非线性优化问题。
为了使用optim和L-BFGS-B方法，首先需要定义一个目标函数，然后使用optim函数来最小化该函数。在这个过程中，优化算法会尝试不同的参数值，并计算函数的梯度来指导搜索方向。
然而，由于计算真实梯度可能需要非常复杂的数学推导，或者计算量很大，所以在实际中往往使用数值方法来近似计算梯度。非有限差分方法是其中一种常用的数值方法，它通过在参数值附近采样多个点，计算函数在这些点上的差分来近似梯度。
非有限差分误差是指使用非有限差分方法计算的梯度与真实梯度之间的差异。由于非有限差分方法是一种近似方法，所以误差是难以完全避免的。误差的大小通常取决于参数的取值范围、差分步长和函数的性质等因素。
要降低非有限差分误差，可以尝试以下几种方法：
- 减小差分步长：通过减小差分步长，可以提高差分的精度，但同时会增加计算量。
- 增加采样点数：通过增加采样点的数量，可以更准确地估计梯度，但也会增加计算量。
- 检查函数实现：确保目标函数的实现正确，并且能够返回正确的梯度值。
- 考虑使用其他优化方法：optim函数还提供了其他一些优化方法，可以尝试使用这些方法，以获得更准确的结果。
在腾讯云中，可以使用腾讯云服务器（CVM）来进行云计算任务，并且可以使用腾讯云提供的弹性MapReduce（EMR）来进行大数据处理。此外，腾讯云还提供了一系列与云计算相关的产品和服务，如腾讯云函数（SCF）、腾讯云容器服务（TKE）等。
参考链接：
- optim函数文档：https://www.rdocumentation.org/packages/stats/versions/3.6.2/topics/optim
- 弹性MapReduce（EMR）产品介绍：https://cloud.tencent.com/product/emr
- 腾讯云函数（SCF）产品介绍：https://cloud.tencent.com/product/scf
- 腾讯云容器服务（TKE）产品介绍：https://cloud.tencent.com/product/tke
相关·内容
1回答
我尝试使用mle来估计负二项分布的参数。这是我的代码。= dnbinom(b, size, prob, log = TRUE)}stats4::mle(LL2, start = list(sizehessian = TRUE, lower = lower, : 当我移除边界时</
浏览 4提问于2021-06-25得票数 1
我正在尝试使用我的对数似然函数来获得五个参数的估计的置信区间(1pi，一个界限在0和1之间的比例，4个betas是指数下降速率)。", hessian = T) 在尝试估计hessian时，即使优化方法本身运行良好，并生成五个参数的估计值(即当hessian = FALSE时)，我仍然得到非有限差分错误。我需要
浏览 74提问于2021-11-03得票数 0
我有病人的生存数据。其中一些被审查了。我想用R中的flexsruv包对经验数据拟合一个广义伽玛函数，而我想要拟合的所有其他函数(如指数、威布尔等)。工作正常，我得到了dist = "gengamma"的以下错误代码 non-finite finite-differ
浏览 3提问于2017-11-03得票数 1
回答已采纳
2回答
我希望对optim()的输出参数进行限制。对于我的正弦函数(其中'designL‘是我的自变量，'ratio’是我的因变量数据，dfm是我的数据帧)，它不必要地收敛数以千计的异相。function(data, par) {}
resultt <- optim(par = c(lo_0, kc_
浏览 1提问于2018-02-06得票数 1
回答已采纳
我正在研究R(初学者)中的损失厌恶模型，并希望从具有3列(损失/增益值(包括连续值和决策编码为0或1(二进制))的数据集中估计一些参数dropbox.com/s/fpw3obrqcx8ld1q/GrandAverage.RData(lambda = 1，temp = temp_s)，nobs = length(set$Decision))系数: lambda 1.298023 1.041057
当我尝试在没有alpha参数的
浏览 0提问于2017-01-05得票数 0
我试图使用p包2.13版本计算rethinking的后验： alist( p ~ dnorm(1e-4, 0.5e-4) data = list(n = 10, N = 10e4),)误差的重新思考：：quap(n~ dbinom(N，p)，p~d范数(1e-04，5e-05))，：非有限差</
浏览 3提问于2020-12-20得票数 2
回答已采纳
在optim函数中，如何在MLE下设置par[1]、par[2]、par[3]的边界？getSymbols("GOOG", from = "2008-06-30", to = "2018-06-30", src = "yahoo")
goog.daily+ log(sigma) + 2*log(
浏览 0提问于2018-08-19得票数 0
我使用optim来估计自定义函数的参数。此自定义函数包括在参数输入无效时产生错误消息的函数(它包括多个分布，但假设如果出现正态分布的负标准差，则产生错误。现在，如果我在BFGS中使用optim选项，那么optim可能在某个时候尝试无效的参数输入，这将导致在计算自定义函数时出现错误，并停止optim过程。如果我<
浏览 2提问于2016-10-05得票数 3
我在R中使用qmap，我有obs数据(Obs)和模型输出(An)，它们都在数据帧中，具有单行和单列。当我使用公式fitQmapDIST(obs,an,dist="GammaDist")时，我得到以下错误：
optim中的错误(par= vstart，fn = fnobj，fix.arg = fix.arg，obs = data，\n values B需要‘fn’\n
浏览 9提问于2015-08-02得票数 0
1回答
这是一个更普遍的问题，有点独立于数据，所以我没有一个MWE。但是，我现在遇到了一个问题，需要我为fn的几个参数之一设置界限。据我所知，optim(method="L-BFGS-B",...)允许我设置参数的</em
浏览 8提问于2019-12-11得票数 1
回答已采纳
1回答
我正在写一篇论文，要求我找到Gumbel的I型指数分布的MLE。我已经证明了的似然函数和日志似然函数，但我很难在r中实现它，以便使用optimization执行优化。我的代码生成NA值。下面是我的密码。}
colnames(data) = c("X1", "X2
浏览 12提问于2022-05-26得票数 -2
回答已采纳
我正在尝试用R的“最优”来最大化可能性。我得到的错误是“非有限差值”。我使用L-BFGS-B是因为我必须约束第11个参数(Bernoulli "p")为0<=p<=1，因为我需要这个约束，所以我不能使用像“Nelder-Mead”这样的非梯度方法。它在模拟数据中工作得很好！
注意，我在这里使用的是floor函数，因为"Tr
浏览 2提问于2017-07-27得票数 1
在R中的BTYD包中执行以下命令时： params <- bgnbd.EstimateParameters(cal.cbs) 其中cal.cbs为下式： cal.cbs <- dc.BuildCBSFromCBTAndDates(cal.cbt, cal.cbs.dates,
per="week") 我收到以下错误: Error in optim(logparams，bg
浏览 6提问于2019-12-02得票数 0
1回答
我尝试通过输入alpha0、alpha1和beta 1的模型及其参数来模拟一个过程。 { xs[i]<-rpois(1,lambda = l[i]) s[j,1:200]<-xs并利用optim函数最小化负对数可能性，从而得到alpha 0、alpha 1和</
浏览 4提问于2022-06-14得票数 0
)/N)}
maxmim <- optimize(func_2, c(0,1), maximum = TRUE) 错误消息不断弹出，提示参数"N“丢失，我该怎么做才能找到p的最佳值
浏览 31提问于2021-04-12得票数 1
回答已采纳
1回答
我试图在R的optim中使用pgtol，而没有得到任何结果。* (y-x^2)*x return(gr)默认优化和检查，梯度高于我想要的。-05 -3.449572e-05
R stats包对optim文档的帮助说，“有助于控制”L B“方法的收敛性。它是对当前搜索方向的投影梯度的公差
浏览 5提问于2022-02-25得票数 3
我的函数有3个参数，x, y, and w。给定w，我需要找到最小化函数的最佳x and y对。但是，x and y被限制在间隔(-0.5 ,0.5]中。我试着使用optim函数，当我没有限制时，它就能工作。然而，当我应用这些限制并使用L B方法时，我会得到一个错误。我搞不懂我为什么会犯这个错误。
下面是我的函数和一些我尝试过的方法</e
浏览 9提问于2022-01-07得票数 1
回答已采纳
1回答
我有一个可以工作的optim函数，它适用于单个行。我想对dataframe中的每一行应用相同的函数(实际上，可能有大量行)。我尝试过使用dplyr中的行方式，但这是不正确的。", lower = c(0), upper = c(100), control = list(maxit = 1000))
但是，当尝
浏览 2提问于2021-04-14得票数 1
回答已采纳
我一直在使用Excel求解器来处理以下问题y = a*b*c*x/((1 - c*x)(1 - c*x + b*c*x))0 < a < 100
0 <我最小化abs的和(y- f(x))。我可以在每个x处限制函数的结果的系数和值的范围。我尝试了nls (没有试图施加约束)，虽然Excel提供了我想要提供<
浏览 6提问于2012-07-09得票数 6
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频
腾讯云开发者
