在 R 中获取反向函数
魔术师大人
我试图在 R 中绘制给定函数的反向函数：
f<-function(x){
if( x>1 || x< -1)
{
0
}else{
0.75*(1-x^2)
}} #densityfunction f
fVec <- Vectorize(f)
F<-function(t){
integrate(fVec, lower=-1, upper=t)$value
}#integral of f over interval -1,1
FVec<-Vectorize(F) #vectorize F
inverse <- function (f, lower = -1, upper = 1) {
function (y) uniroot((function (x) f(x) - y), lower = lower, upper = upper)$root
}
Finv = inverse(F, -1, 1)
FinvVec<-Vectorize(Finv) # Vectorize
#plot(FVec, xlim=c(-2, 2)) #plot F
plot(FinvVec, xlim=c(-2, 2)) #plot F inv
我的问题是我收到错误：
Error in uniroot((function(x) f(x) - y), lower = lower, upper = upper) :
f() values at end points not of opposite sign
AFAIK 这意味着我的函数没有符号变化。我问了我的导师，他告诉我有一个标志变化，我在这里做错了。但我就是不知道是什么。你们能帮忙吗？
安德鲁·莱克
如果您将调用包装uniroot在 a 中try，那么它会在发现错误时继续运行
inverse <- function (f, lower = -1, upper = 1) {
function (y) try(uniroot((function (x) f(x) - y), lower = lower, upper = upper)$root)
}
结果图表明它仅适用于 (0,1) 中的 x。
plot(FinvVec, xlim=c(-2, 2))
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
编辑于
相关文章
Related 相关文章
- 1
在C ++中获取反向列表
- 2
如何在python2中获取反向枚举？
- 3
在rest框架中获取反向关系对象
- 4
反向cbind（）函数R
- 5
反向cbind（）函数R
- 6
在R中的咖喱函数中查找（反向递归？）
- 7
在R中的咖喱函数中查找（反向递归？）
- 8
在R中的gsub中将函数应用于反向引用
- 9
R中正态分布函数的反向/反向
- 10
R中的反向cumprod
- 11
R中的反向引用
- 12
Django ORM：获取反向指示
- 13
在R中提取反事实ARIMA预测的函数
- 14
无法从Django中的第三方软件包获取反向URL
- 15
无法从Django中的第三方软件包获取反向URL
- 16
如何在R中对dydyverse中的分数取反分位数？
- 17
获取R中的函数名称
- 18
影像r中的反向轴
- 19
R中矩阵的反向索引
- 20
在R中反向子集向量
- 21
在JavaScript中取反布尔函数
- 22
给定R env中定义的函数，获取函数参数
- 23
给定R env中定义的函数，获取函数参数
- 24
Queryset，仅获取反向关系之一
- 25
从完整路径获取反向URL名称
- 26
从完整路径获取反向URL名称
- 27
Queryset，仅获取反向关系之一
- 28
使用 netcat 通过 UDP 获取反向 shell
- 29
如何在R中的调用函数环境中获取对象？
我来说两句