作者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