similar to: Exactness of ppois

Full_Name: Morten Welinder Version: 2 OS: Solaris/space/gcc2.95.2 Submission from: (NULL) (65.213.85.217) I changed src/nmath/standalone/test.c to read: --------------------------------------------------------------------------------- #define MATHLIB_STANDALONE 1 #include <Rmath.h> #include <stdio.h> int main() { double x; for (x = 99990; x <= 100009; x++) printf

Hello, I am sure I am not the only person with this problem. I have a list with n elements, each consisting of a single column matrix with different row lengths. Each row has a name ranging from A to E. Here is an example: alph[[1]] A 1 B 2 C 3 D 4 alph[[2]] A 1 C 3 D 4 alph[[3]] A 1 D 4 E 5 I would like to create a matrix from the elements in the list with n columns such that the row names

function ppois is a function calculate the CDF of Poisson distribution, it should generate a non-decreasing result, but what I got is: > any(diff(ppois(0:19,lambda=0.9))<0) [1] TRUE Actually, > ppois(19,lambda=0.9)<ppois(18,lambda=0.9) [1] TRUE Which could not be TRUE. Code is tested in both R 3.5.1 and Microsoft R Open 3.5.1. _

>>>>> Serguei Sokol >>>>> on Tue, 4 Dec 2018 11:46:32 +0100 writes: > Le 04/12/2018 ? 11:27, I?aki Ucar a ?crit?: >> On Tue, 4 Dec 2018 at 11:12, <qweytr1 at mail.ustc.edu.cn> wrote: >>> function ppois is a function calculate the CDF of Poisson distribution, it should generate a non-decreasing result, but what I got is:

A quick question! The number of episodes per year of otitis media follows a Possion distribution with lambda = 1.6 episodes per year. Wouldn't the probability of getting 3 or more episodes of otitis media in the first 2 years of life be: > ppois(q=3, lambda=1.6*2, lower.tail = TRUE, log.p = FALSE) [1] 0.6025197 I am confused with the lambda and 3 or more.. thx much

Le 04/12/2018 ? 11:27, I?aki Ucar a ?crit?: > On Tue, 4 Dec 2018 at 11:12, <qweytr1 at mail.ustc.edu.cn> wrote: >> function ppois is a function calculate the CDF of Poisson distribution, it should generate a non-decreasing result, but what I got is: >> >>> any(diff(ppois(0:19,lambda=0.9))<0) >> [1] TRUE >> >> Actually, >> >>>

Full_Name: Murray H Smith Version: 0.63.3 OS: Windows NT Submission from: (NULL) (130.216.5.57) The ppois function is displaced by -0.5. Try: > ppois(-0.5,1) [1] 0.3678794 > ppois(-0.51,1) [1] 0 > ppois(0,1) [1] 0.3678794 and > par(mfrow=c(2,1)) > x<-seq(-1,5,0.01) > plot(x,ppois(x,1),type="s",ylab="F(x)",main="Poisson CDF?") >

Here's a contribution from Ian Smith that got bounced from the list. -------- Original Message -------- Subject: Re: [Rd] pt inaccurate when x is close to 0 (PR#9945) Date: Thu, 11 Oct 2007 06:02:43 -0400 From: iandjmsmith at aol.com To: murdoch at stats.uwo.ca Duncan, I tried sending the rest of this to R-devel but it was rejected as spam, hence the personal e-mail. R calculates the pt

Thanks Ben for verifying the issue. It is always reassuring to hear when others can reproduce the problem. I wrote a small patch that fixes the issue (https://github.com/r-devel/r-svn/pull/11): diff --git a/src/nmath/qnbinom.c b/src/nmath/qnbinom.c index b313ce56b2..d2e8d98759 100644 --- a/src/nmath/qnbinom.c +++ b/src/nmath/qnbinom.c @@ -104,6 +104,7 @@ double qnbinom(double p, double size,

Full_Name: Morten Welinder Version: 2.1.0 OS: src only Submission from: (NULL) (216.223.241.212) Now that pcauchy has been fixed, it is becoming clear that qcauchy suffers from the same problems. qcauchy(pcauchy(1e100,0,1,FALSE,TRUE),0,1,FALSE,TRUE) should yield 1e100 back, but I get 1.633178e+16. The code below does much better. Notes: 1. p need not be finite. -Inf is ok in the log_p

Full_Name: M Welinder Version: 1.4 OS: (src) Submission from: (NULL) (192.5.35.38) It seems to me that pweibull can be improved in the lower_tail=TRUE and log_p=FALSE case by using expm1. Something like -expm1(-pow(x / scale, shape)), I think. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read

I have DHCP and TFTP on a Solaris machine, pxelinux.0 loads, but it can't load the configuration file. The tftp server on Solaris tries to send the file, but data transfers do not receive ACKs. Previous posts to this list have indicated that the suspect was the PXE BIOS (ref "PXELinux can't load the config file?" dated March 20, 2003). Does that still appear to be the case

Hi all, I recently noticed that `qnbinom()` can take a long time to calculate a result if the `size` argument is very small. For example qnbinom(0.5, mu = 3, size = 1e-10) takes ~30 seconds on my computer. I used gdb to step through the qnbinom.c implementation and noticed that in line 106 (https://github.com/wch/r-source/blob/f8d4d7d48051860cc695b99db9be9cf439aee743/src/nmath/qnbinom.c#L106)

Hello All, I have been trying to do bonpow for a while now.Bonpow works for few graphs and for few others it goes out of memory. I did see reply to one of the posting Alph Centrality Crashed in Igraph memory error. The solution in the posting was to use sparse matrix. This is the link of the message. http://lists.gnu.org/archive/html/igraph-help/2008-04/msg00071.html I was able to obtain

Dear developement team, I want to share a massive storage casted with Ceph by samba with windows workstations. All works well. My problem so is that in windows the ceph storage size statistics are wrong. Instead of seeing a 44TB hard drive I see a 176GB hard drive. Under linux that issue doesn't show. The size are properly reported. I investigated around and it seems that the problem

Hi! I fairly new to R, have only done pretty basic things so far, so this may be a very basic question..... But I did search the forums and didn't see a solution.... I'm trying to get going with kripp.alpha(). I'm loading data from a file, like this: > library(irr) Loading required package: lpSolve > x <- read.table("foo", comment.char="#", header=TRUE,

x <- c( "a1", "a10", "a2" ) y <- c( "b10", "b2", "a12", "ca1" ) DF <- expand.grid( x = x, y = y ) # randomize set.seed( 42 ) DF <- DF[ sample( nrow( DF ) ), ] # missing from gtools mixedrank <- function( x ) { seq.int( length( x ) )[ gtools::mixedorder(x) ] } o <- do.call( order, lapply( DF, mixedrank ) )

Hi I have started using rails eg: rform and typo, but mostly I do flash dev (I know you guys decided you did not want flash content and used the name for something else!). Anyway what experiments have people been doing, anyone been Ajaxing flash content within rails application or even just embedding in rails, I did not really find anything on using flash in Typo, not that you can''t

Dear expeRts, I would like to create a scatter plot matrix with splom(). The lower panel should contain some additional information about the samples shown in the upper panel plot, see the splom() call below. Now two questions came up: (1) The lower panels show "tau" and "alpha" on top of each other. How can I plot *three* expressions on top of each other? I tried several

Full_Name: Morten Welinder Version: snapshot OS: Submission from: (NULL) (65.213.85.129) Two things are wrong. 1. There is nan test outside IEEE_754. 2. The meat part of the function should really be something like... if (!lower_tail) x = -x; if (fabs (x) > 1) { double temp = atan (1 / x) / M_PI; return (x > 0) ? R_D_Clog (temp) : R_D_val (-temp); } else