作者slow0201 (海星)
看板R_Language
标题[问题] 积分计算值
时间Thu Dec 12 12:13:20 2013
大家好,有个问题想请教
描述如下:
函数为两个常态的CDF和一个常态的PDF的乘积
积分下限为负无穷大,上限为无穷大
请问直接使用下面程式积分是否容易产生错误结果?
我跑出来的结果比论文提供的值低很多
请问是否有别的写法可以使得我得到的结果较接近
论文提供的结果,谢谢!
--------------------------------------------
程式编写为
x=1
y=1.486607
z=1
w=1.486607
fun1<-function(u)
{
fun_ev=pnorm(x*u+y)*pnorm(-z*u+w)*dnorm(u)
return(fun_ev)
}
AS=integrate(fun1,lower=-Inf,upper=Inf)$value
----------------------------------------------
出来的结果AS是0.7096956
----------------------------------------------
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.109.73.41
1F:→ celestialgod:我帮你用MATLAB算也是0.7097 12/12 12:39
2F:→ celestialgod:你应该去看一下你的FUNCTION也没有写错 12/12 12:39
3F:→ celestialgod: 有 12/12 12:39
4F:→ celestialgod:其他FUNCTION如quadinf (package:pracma) 12/12 12:41
5F:→ celestialgod:AS2=quadinf(fun, -Inf, Inf, tol=1e-10);AS2=0.7097 12/12 12:43
6F:→ slow0201:感谢楼上的帮忙,实在是不知道文章中的0.87是怎麽算出来的 12/12 13:35
7F:→ psinqoo:sum(AS) [1] 0.7096957 summary(AS) 0.7097 12/12 16:19