[問題]R 多維常態的cdf運算
在R的mvtnorm封包可運算多維常態的cdf，他的求算方法
不是向量的觀念。
在單維裡面
x = runif(100)
pnorm(x)
他是向量的觀念，不需要用loop的方式求數cdf，
但這招在mvtnorm裡卻行不通，我只能透過loop的方式求cdf。
require(mvtnorm)
obs=1000
xb1 = runif(obs)
xb2 = runif(obs)
rho = 0.2
ll = numeric(obs)
for (i in 1:obs){
ll[i]=
pmvnorm(lower= -Inf,upper=c(xb1[i],-xb2[i]),
corr=matrix(c(1,-rho,-rho,1),2,2))
}
請問是否有其它方法可省去loop(apply?)，更有效率求算數值。
謝謝^^"
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.86.136
→
07/09 02:01, , 1F
07/09 02:01, 1F