Displaying 1 result from an estimated 1 matches for "distr1".
Did you mean:
distro
2010 Nov 07
3
Integrate and mapply
...lt;-cbind(t1,t2,t3,t4)
nu<-rchisq(15,4)
gam<-matrix(rchisq(75,df=5),15,5)
undint<-function(u){
prob<-function(i) {
i1<-comb[i,1]
i2<-comb[i,2]
i3<-comb[i,3]
i4<-comb[i,4]
val1<-gam[((i3-1)*2+1):(i3*2),i4]
vc1<-matrix(12,length(val1),1)-matrix(u,length(val1),1)-val1
distr1<-prod(pchisq(vc1,df=2))
p10<-distr1*1/(exp(-nu[i1])*(1+gam[i2,i4]))*(1-plogis((u+log(gam[i3,i4])-log(1+gam[i2,i4])),0,1))
p10
}
val<-c(1:15)
q<-sapply(val,prob)
}