Displaying 1 result from an estimated 1 matches for "yinv".
Did you mean:
ying
2001 Oct 28
0
Summary: A speed improvement challenge
...umsum(x)/sum(x),rep(1,n))[1:n],n=n) #
make
pcum's all same length
apply(sweep(pcum,2,runif(2),">"),2,function(x,n)(1:n)[x][1],n=n)
Nick Ellis
CSIRO Marine Research mailto:Nick.Ellis at csiro.au
-------------------------------------------------------------------------------
try:
yinv<-sapply(1:m,function(i,x,y,freq,aty) ifelse(any(s<-abs(y-aty[i]) <
del),
sample(x[s], 1, replace=F, prob=freq[s]/sum(freq[s])),
approx(y, x, xout=aty[i], rule=2)$y),
x,y,freq,aty)
===================================================================
Vadim Kutsyy, PhD...