作者bcs (= ="frailty..gggg XD)
看板Statistics
标题[问题]R 多维常态的cdf运算
时间Wed Jul 8 13:36:37 2009
在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
1F：→ clickhere:改它的 C/FORTRAN code. 07/09 02:01