similar to: PR#8528

R versions 2.1.0 to present. Examples shown were computed under Windows R-devel, current SVN, but ix86 Linux shows similar behaviour (sometimes NaN or -Inf rather than Inf, depending on the compiler and optimization level used). The replacement pgamma algorithm used from R 2.1.0 has an inadequate design and no supporting documentation whatsoever. There is no reference given to support the

(1) ...need to speed up a monte-carlo sampling...any suggestions about how I can get R to use all 8 cores of a mac pro would be most useful and very appreciated... (2) spent the last few hours trying to get pnmath to compile under os- x 10.5.4... using gcc version 4.2.1 (Apple Inc. build 5553) as downloaded from CRAN, xcode 3.0... ...xcode 3.1 installed over top of above after

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

On a recent FreeBSD 8.0-CURRENT (i386) building R (any version) breaks with the following messages: ---------------------------------------------------------------------- [...snip...] gcc -std=gnu99 -I. -I../../src/include -I../../src/include -I/usr/local/include -DHAVE_CONFIG_H -g -O2 -c wilcox.c -o wilcox.o gcc -std=gnu99 -I. -I../../src/include -I../../src/include -I/usr/local/include

There is a problem with `rpois'. It does seem to take care about the order of the arguments. This is an example: > rpois(n=1,lambda=2) [1] 3 > rpois(lambda=2,n=1) [1] 2 0 It obviously uses the first argument as the number of samples to be drawn, which is wrong. I used Version 0.49 Beta (April 23, 1997). Fredrik

There is a problem with `rpois'. It does seem to take care about the order of the arguments. This is an example: > rpois(n=1,lambda=2) [1] 3 > rpois(lambda=2,n=1) [1] 2 0 It obviously uses the first argument as the number of samples to be drawn, which is wrong. I used Version 0.49 Beta (April 23, 1997). Fredrik

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,

I came across the message below and decided to respond. I don't use R but I do use other FSF products and I am a believer in making high quality software freely available, particularly for basic mathematical and statistical functions. I can thorougly recommend TOMS Algorithm 708 but I believe there are a couple of problems from you point of view. Firstly there may be restrictions on how it

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)

Hi! I have just downloaded the R-0.60.1 sources and have problems compiling R on a Sun Ultra 1 running Solaris 2.6 and gcc 2.7.2.1. I have not been able to find to find any compiling hints in the documentation or the FAQ. After ./configure I use make and get the output below. Any hints are welcome. I am not on the list, so please answer me directly too. Best regards Jens --- Jens Lund

I have been looking around for VoIP providers but have not found a good listing. Is there no "yellow pages" for VoIP providers? Google mostly returns services like Vonage, Packet8, NuFone, ect. None seam to be very reseller friendly and none offer LNP or local DID's for my area. Anyone know of a list (even a partial one) Jeromie Reeves

Hi Martin, thanks for verifying. I agree that the Cornish-Fisher seems to struggle with the small size parameters, but I also don't have a good idea how to replace it. But I think fixing do_search() is possible: I think the problem is that when searching to the left y is decremented only if `pnbinom(y - incr, n, pr, /*l._t.*/TRUE, /*log_p*/FALSE)) < p` is FALSE. I think the solution is

Hi, I am looking for a provider that can provide DIDs in Atlanta, GA (national in future) and who can port an existing number. LiveVoip is a perfect fit but doesn't do LNP :-(. Currently I only need a handful of numbers ported, but will have more requirement in future. I'll need unlimited simultaneous incoming/outgoing calls (can pay for incoming as well). Any pointers are much

Hello Mates, I am facing a really weird issue with my mail server, somehow I cannot connect using port 993. It works with 143 but not 993 nor 995. Here is a little bit of more information: # telnet localhost 143 Trying ::1... telnet: connect to address ::1: Connection refused Trying 127.0.0.1... Connected to localhost. Escape character is '^]'. * OK [CAPABILITY IMAP4rev1 LITERAL+

Hi, I am trying to obtain the heteroskedasticity consitent standard errors (HCSE) after apply 2SLS. I obtain 2SLS through tsls from package sem or systemfit: #### tsls #### library (sem) Reg2SLS <-tsls(LnP~Sc+Ag+Ag2+Var+R+D,~I2+Ag+Ag2+Var+R+D) summary (Reg2SLS) #### systemfit #### library (systemfit) RS <- LnP~Sc+Ag+Ag2+Var+R+D Inst <- ~I2+Ag+Ag2+Var+R+D labels

Hi again ! Thanks for the info on updating the config.site file which I have done. I have also added -lm in the Makeconf manually because this is needed explicitly for DEC cc. However, there are still a few problems when linking some of the files as you can see from the enclosed log. Ciao, Andreas ------------------------------------------------------- R-0.63.3>make make[1]: Entering

Hi, Just downloaded the R-0.62.4 of R and tried to install it. With the standard procedure : ./configure make At the end of the compilation I got the following message : ld: ../lib/libunix.a(system.o): main: multiply defined fort: Severe: Failed while trying to link. *** Exit 1 Stop. *** Exit 1 Stop. *** Exit 1 Stop. *** Exit 1 Stop. I attach the printout after the ./configure and make

Dear R Users, I am trying to users gsub to remove multiple cases of square brackets and their different contents in a character string. A sample of such a string is shown below. However, I am having great difficulty understanding regexp syntax. Any help is greatly appreciated. Ally "tree STATE_286000 [&lnP=-12708.453945423369] = [&R]

This is a multipart MIME message. --==_Exmh_981436288450 Content-Type: text/plain; charset=iso-8859-1 Content-Transfer-Encoding: quoted-printable Hi ! I have problems compiling R successfully on a DEC-UINX 4.0E. I have applied the recommended config.site, which I enclose. As can be seen from the compilation log there are linking errors... I did a 'make check' which fails for the R

Mike, No they have not. Calls are failing again today. They have offered to refund my money but that does not solve the problem. My asterisk server is only 4 to 12 ms away from their "network". I have had VERY good luck with nufone.(40 to 45ms away) Only have 1 or 2% fail rate. Going to be calling txlink.net on Monday. Seems that LiveVoIP does not care about asterisk users. They like