uniroot
位于 stats
包(package)。 说明
函数uniroot
在从lower
到upper
的区间中搜索函数f
相对于其第一个参数的根(即零)。
将extendInt
设置为非"no"
字符串,意味着如果sign(f(x))
在区间终点不满足要求,则寻找正确的interval = c(lower,upper)
;请参阅“详细信息”部分。
用法
uniroot(f, interval, ...,
lower = min(interval), upper = max(interval),
f.lower = f(lower, ...), f.upper = f(upper, ...),
extendInt = c("no", "yes", "downX", "upX"), check.conv = FALSE,
tol = .Machine$double.eps^0.25, maxiter = 1000, trace = 0)
参数
f |
求根的函数。 |
interval |
包含要搜索根的区间的 end-points 的向量。 |
... |
要传递给 |
lower , upper |
要搜索的区间的下端点和上端点。 |
f.lower , f.upper |
分别与 |
extendInt |
字符串,指定当 |
check.conv |
逻辑指示底层 |
tol |
所需的精度(收敛容差)。 |
maxiter |
最大迭代次数。 |
trace |
整数;如果是肯定的,则产生追踪信息。更高的值提供更多细节。 |
细节
请注意,...
之后的参数必须完全匹配。
必须指定interval
或同时指定lower
和upper
:上端点必须严格大于下端点。对于 extendInt="no"
来说,端点处的函数值必须具有相反的符号(或零),这是默认值。否则,如果 extendInt="yes"
,则在两侧扩展间隔,以搜索符号更改,即,直到搜索间隔 满足 。
如果知道extendInt
可以(通常应该)指定为"upX"
(对于“upward crossing”)或"downX"
,分别。同样,定义 ,以在解决方案中要求 。在这种情况下,搜索间隔 可能被扩展为使得 和 。 如何在根 处改变符号,也就是说,如果函数在那里增加或减少,则
uniroot()
使用基于下面参考文献中给出的算法的 Fortran 子例程 zeroin
(来自 Netlib)。它们假设一个连续函数(已知该函数在区间内至少有一个根)。
如果 f(x) == 0
或算法一步的 x
变化小于 tol
(加上 x
中的表示错误容限),则声明收敛。
如果算法在 maxiter
步骤中未收敛,则会打印警告并返回当前近似值。
f
将被称为f(x, ...)
对于数值x.
传递给 f
的参数具有特殊语义,用于在调用之间共享。该函数不应该复制它。
值
至少包含五个组件的列表:root
和 f.root
给出根的位置以及在该点计算的函数的值。 iter
和 estim.prec
给出使用的迭代次数以及 root
的近似估计精度。 (如果根出现在端点之一,则估计精度为 NA
。) init.it
包含初始 extendInt
迭代次数(如果有),否则为 NA
。在这种extendInt
迭代的情况下,iter
包含其中的sum
和zeroin
迭代。
将来可能会添加更多组件。
例子
require(utils) # for str
## some platforms hit zero exactly on the first step:
## if so the estimated precision is 2/3.
f <- function (x, a) x - a
str(xmin <- uniroot(f, c(0, 1), tol = 0.0001, a = 1/3))
## handheld calculator example: fixed point of cos(.):
uniroot(function(x) cos(x) - x, lower = -pi, upper = pi, tol = 1e-9)$root
str(uniroot(function(x) x*(x^2-1) + .5, lower = -2, upper = 2,
tol = 0.0001))
str(uniroot(function(x) x*(x^2-1) + .5, lower = -2, upper = 2,
tol = 1e-10))
## Find the smallest value x for which exp(x) > 0 (numerically):
r <- uniroot(function(x) 1e80*exp(x) - 1e-300, c(-1000, 0), tol = 1e-15)
str(r, digits.d = 15) # around -745, depending on the platform.
exp(r$root) # = 0, but not for r$root * 0.999...
minexp <- r$root * (1 - 10*.Machine$double.eps)
exp(minexp) # typically denormalized
##--- uniroot() with new interval extension + checking features: --------------
f1 <- function(x) (121 - x^2)/(x^2+1)
f2 <- function(x) exp(-x)*(x - 12)
try(uniroot(f1, c(0,10)))
try(uniroot(f2, c(0, 2)))
##--> error: f() .. end points not of opposite sign
## where as 'extendInt="yes"' simply first enlarges the search interval:
u1 <- uniroot(f1, c(0,10),extendInt="yes", trace=1)
u2 <- uniroot(f2, c(0,2), extendInt="yes", trace=2)
stopifnot(all.equal(u1$root, 11, tolerance = 1e-5),
all.equal(u2$root, 12, tolerance = 6e-6))
## The *danger* of interval extension:
## No way to find a zero of a positive function, but
## numerically, f(-|M|) becomes zero :
u3 <- uniroot(exp, c(0,2), extendInt="yes", trace=TRUE)
## Nonsense example (must give an error):
tools::assertCondition( uniroot(function(x) 1, 0:1, extendInt="yes"),
"error", verbose=TRUE)
## Convergence checking :
sinc <- function(x) ifelse(x == 0, 1, sin(x)/x)
curve(sinc, -6,18); abline(h=0,v=0, lty=3, col=adjustcolor("gray", 0.8))
uniroot(sinc, c(0,5), extendInt="yes", maxiter=4) #-> "just" a warning
## now with check.conv=TRUE, must signal a convergence error :
uniroot(sinc, c(0,5), extendInt="yes", maxiter=4, check.conv=TRUE)
### Weibull cumulative hazard (example origin, Ravi Varadhan):
cumhaz <- function(t, a, b) b * (t/b)^a
froot <- function(x, u, a, b) cumhaz(x, a, b) - u
n <- 1000
u <- -log(runif(n))
a <- 1/2
b <- 1
## Find failure times
ru <- sapply(u, function(x)
uniroot(froot, u=x, a=a, b=b, interval= c(1.e-14, 1e04),
extendInt="yes")$root)
ru2 <- sapply(u, function(x)
uniroot(froot, u=x, a=a, b=b, interval= c(0.01, 10),
extendInt="yes")$root)
stopifnot(all.equal(ru, ru2, tolerance = 6e-6))
r1 <- uniroot(froot, u= 0.99, a=a, b=b, interval= c(0.01, 10),
extendInt="up")
stopifnot(all.equal(0.99, cumhaz(r1$root, a=a, b=b)))
## An error if 'extendInt' assumes "wrong zero-crossing direction":
uniroot(froot, u= 0.99, a=a, b=b, interval= c(0.1, 10), extendInt="down")
来源
基于 'zeroin.c' 在https://netlib.org/c/brent.shar.
参考
Brent, R. (1973) Algorithms for Minimization without Derivatives. Englewood Cliffs, NJ: Prentice-Hall.
也可以看看
相关用法
- R update 更新并重新拟合模型调用
- R update.formula 模型更新
- R stlmethods STL 对象的方法
- R medpolish 矩阵的中值波兰(稳健双向分解)
- R naprint 调整缺失值
- R summary.nls 总结非线性最小二乘模型拟合
- R summary.manova 多元方差分析的汇总方法
- R formula 模型公式
- R nls.control 控制 nls 中的迭代
- R aggregate 计算数据子集的汇总统计
- R deriv 简单表达式的符号和算法导数
- R kruskal.test Kruskal-Wallis 秩和检验
- R quade.test 四方测试
- R decompose 移动平均线的经典季节性分解
- R plot.stepfun 绘制阶跃函数
- R alias 查找模型中的别名(依赖项)
- R qqnorm 分位数-分位数图
- R eff.aovlist 多层方差分析的计算效率
- R pairwise.t.test 成对 t 检验
- R loglin 拟合对数线性模型
- R predict.smooth.spline 通过平滑样条拟合进行预测
- R bartlett.test 方差齐性的 Bartlett 检验
- R influence.measures 回归删除诊断
- R loess.control 设置黄土参数
- R Normal 正态分布
注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 One Dimensional Root (Zero) Finding。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。