[問題] 積分計算值
大家好,有個問題想請教
描述如下:
函數為兩個常態的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
→
12/12 12:39, , 1F
12/12 12:39, 1F
→
12/12 12:39, , 2F
12/12 12:39, 2F
→
12/12 12:39, , 3F
12/12 12:39, 3F
→
12/12 12:41, , 4F
12/12 12:41, 4F
→
12/12 12:43, , 5F
12/12 12:43, 5F
→
12/12 13:35, , 6F
12/12 13:35, 6F
→
12/12 16:19, , 7F
12/12 16:19, 7F
R_Language 近期熱門文章
PTT數位生活區 即時熱門文章
