我有一个函数func
(在R中)有两个像这样定义的变量(只是一个例子):
func <- function(t,gam){return(exp(-t/2)*(t+1)^gam)}
第二个参数gam
将从2个不同的样本估计,因此基本上gam
是长度为2的向量。我想要接收的是t函数的集成:0-&gt;知道gam
时的Inf。
我尝试了几种方法但无法从集成中获取值向量。任何建议将非常感谢。非常感谢。
答案 0 :(得分:0)
需要使用sapply
将两个项目向量(单独)传递给integrate
:
sapply( gam, function(x) integrate( func, lower=0,upper=Inf, gam=x)$value )
示例:
> sapply( c(.5,1), function(x) integrate( func, lower=0,upper=Inf, gam=x)$value )
[1] 3.311359 6.000000