r - 为什么 "Uniroot"命令给出不同的结果

我有一个关于 uniroot 命令的问题。我无法以合乎逻辑的方式提出我的问题，因为我不知道为什么在下面示例的第二种情况下每次结果都不同。
在第一种情况下，我的 f 函数的输出总是相同的:

library(mvtnorm)
f <- function(y1_upr,y2_upr = -0.05453663,target = 25e-4,df=3) {
    pmvt(upper = c(y1_upr,y2_upr),df = df) - target
}
uniroot(f,c(-10000,10000))$root

但我不知道为什么当我在同一个函数中添加另一个参数时，结果每次都会改变(见下文):

g <- function(y1_upr,
              y2_upr = -0.05453663,
              y3_upr = -0.06236616,
              target = 25e-4,
              df = 3) {
    pmvt(upper = c(y1_upr,y2_upr,y3_upr),df = df) - target
}
uniroot(g,c(-10000,10000))$root

您会看到，每次使用 g 函数应用 uniroot 命令时，都会得到不同的结果。
有没有人对此有想法？我能找到一种方法来解决我的结果吗？

最佳答案

该算法可能涉及一些随机性措施，以便选择一个起点。尝试:

set.seed(1)
uniroot(g,c(-10000,10000))$root

关于r - 为什么 "Uniroot"命令给出不同的结果，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/14348438/