similar to: Fwd: Rcpp sugar dpois

Dear r developers, I'm writing a set of new functions for an existing R package on R-forge (called COGARCH). I wrote the new R code (but still no documentation), and updated the NAMESPACE file. I installed it from my local repository and everything seems to work. I checked the local repository and the check did not produces errors, but 2 warnings and 2 notes (I attach the log) one of which

I would like to write a fast version of a function which contains slow loops that I cannot vectorize. I would like to rewrite it in C, but in the loops I need to call a few (two) R functions. Which is the best way to do so? I have read that one can call a C library which import all the R functions (R.h, is it?). Is that efficient? I tried to have a look into the internal C code to find just the

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130402/fd10647d/attachment.pl>

Chuck and Greg, Thanks a lot for your help! I have a much better understanding now of what is happening ?under the hood?. Kind Regards, Jason > On 31 May 2018, at 20:08, Greg Minshall <minshall at acm.org> wrote: > > Jason, > > as Chuck Berry (to whom, *thanks* for 'do {...} while(0)'!) suggested, > using grep, or even grep executed from find, such as >

Hi - I have a likelihood function that involves sums of two possions: L = a*dpois(Xi,theta1)*dpois(Yi,theta2)+b*(1-c)*a*dpois(Xi,theta1+theta3)*dpois(Yi,theta2) where a,b,c,theta1,theta2,theta3 are parameters to be estimated. (Xi,Yi) are observations. However, Xi and Yi are usually big (> 20000). This causes dpois to returns 0 depending on values of theta1, theta2 and theta3. My first

The nice new log=TRUE option in dpois appears to mess up the case where lambda=0 (I was trying to calculate the likelihood of a saturated model). Because the behavior is now always to calculate the probability in terms of exp(log(prob)), there's a test for lambda<=0 which really needs to be only lambda<0. dpois(0:5,0) ought to give 1 0 0 0 0 but gives NaNs instead. Here's

Hello all, I am trying to get a better understanding of the underlying code for the stats::dpois function in R and, specifically, what happens under the hood when it is called. I finally managed to track down the C course at: https://github.com/wch/r-source/blob/trunk/src/nmath/dpois.c. It would seem that the dpois C function is taking a double for each of the x and lambda arguments so I am a bit

Dear Masters, I have a question to submit consider the following script m<-4.95 obs<-rpois(36,m) # i generate 36 realization from a poisson(m) hist(obs,freq=F) curve(dpois(x,m),add=T,col="red") #i wish to overlay on the histogram the theorical poisson density function errors are returned saing the x vector doesn't contain integers.... really bizarre i can't give

depends on "autoconf: regenerate configure scripts with 4.4 version" This series removes some of the really old deadwood from the tools build and makes some other things which are on their way out configurable at build time with a default depending on how far down the slope I judge them to be. * nuke in tree copy of libaio * nuke obsolete tools: xsview, miniterm, lomount & sv *