similar to: Return value from a C program is different from the R console

Hello, I have been working with various probability distributions in R, and it seems the gamma distribution is inaccurate for some inputs. For example, qgamma(1e-100, 5e-101, lower.tail=FALSE) gives: 1.0. However, it seems this is incorrect; I think the correct answer should be 0.082372029620717283. When I check these numbers using pgamma, I get: pgamma(1,5e-101, lower.tail=FALSE) =

Full_Name: Version: 2.7.1 (2008-06-23) OS: windows vista Submission from: (NULL) (216.82.144.137) Hello, I have been working with various probability distributions in R, and it seems the gamma distribution is inaccurate for some inputs. For example, qgamma(1e-100, 5e-101, lower.tail=FALSE) gives: 1.0. However, it seems this is incorrect; I think the correct answer should be

I haven't looked inside to see what is causing this, but there's a big discontinuity in qgamma: curve(qgamma(x, shape=19), from=1e-10, to=2e-10) This appears in both R-patched and R-devel. Duncan Murdoch

David, Duncan, > Hi > > Following on David's rate argument > > try (with modifications of pch and grid) > > rate <- 1/4 > shape = 8 > rate = c(rep(1/4,100),rep(1/3,100)) I don't think the problem is related to the rate argument, which can well be vectorized, as is the case for a number of arguments in distrib functions in R (note that you are redefining it

dgamma and friends in S are documented as dgamma(x, shape, rate=1) pgamma(q, shape, rate=1) qgamma(p, shape, rate=1) rgamma(n, shape, rate=1) whereas R has dgamma(x, shape, scale=1, log = FALSE) pgamma(q, shape, scale=1, lower.tail = TRUE, log.p = FALSE) qgamma(p, shape, scale=1, lower.tail = TRUE, log.p = FALSE) rgamma(n, shape, scale=1) Note the use of rate vs scale. Indeed, as both S and

Hi I need to find the Hessian matrix for a complicated function from a certain kind of data but i keep getting this error Error in f1 - f2 : non-numeric argument to binary operator the data is given by U<-runif(n) Us<-sort(U) tau1<- 2 F1tau<- pgamma((tau1/theta1),shape,1) N1<-sum(Us<F1tau) X1<- Us[1:N1]

Hello! I am trying with this question again: I would like to test few distributional assumptions for some behavioral response data. There are few theories about true distribution of those data, like: normal, lognormal, gamma, ex-Gaussian (exponential-Gaussian), Wald (inverse Gaussian) etc. The best way would be via qq-plot, to show to students differences. First two are trivial: qqnorm(dat$X)

Hi all ! Thanks for the very good support of Dell R320 hardware, perc H310 is well supported, BCM5720 seems to work correctly and performances are great. The only problem i have found is very strange. The FreeBSD bootloader take many times to load, 30sec-2minutes to boot the kernel and show the bootloader menu. After that, the system boots properly, at a normal speed. Is there any issue or

Hello ALL! Can anyone tell me how to specify qqmath() function with distribution qgamma? In help it is bit vague how to pass shape and scale parameters for gamma. Thank you in advance. Best, PM

Folks, I'm fairly sure that I'm doing something stupid, but I'm getting a few really strange results from *some* of the distributions, but by no means all, when I link directly to the R shared library. I've tried this on both Windows with the precompiled Mingw binary of R-2.5.0 (compiling my code with MinGW-3.4.2), and by building R-2.5.0 on Mandriva Linux with gcc-3.4.4 and

Hi, I just made a custom fact to display a status about mysql. The fact variable is ''foo''. In a module I test this fact variable to apply or not a ressource : class bar { if $::foo == ''toto'' { notice("Execute ressource") # Reload fact variable ? if $::foo == ''toto'' { fail("Ressource not

R-2.10.1.tar.gz was built a short while ago. This is a maintenance release and fixes a number of mostly minor issues. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.10.1.tar.gz or wait for it to be mirrored at a CRAN site nearer to you. Binaries for various platforms will appear in due course. For the R Core Team Peter

R-2.10.1.tar.gz was built a short while ago. This is a maintenance release and fixes a number of mostly minor issues. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.10.1.tar.gz or wait for it to be mirrored at a CRAN site nearer to you. Binaries for various platforms will appear in due course. For the R Core Team Peter

(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

Hello, I try to migrate a VM and I get this error : loic at loic-VirtualBox:~/workspace/src/build/bin$ ./virsh --connect qemu:///system migrate --live vmUbuntuPartage2CacheWriteBack qemu+ssh://192.168.10x.2x/system --unsafe loic at 192.168.10x.2x's password: error: invalid argument: qemuDomainMigratePrepare3: unsupported flags (0x200) Does somebody know this error? I am working on windows

https://bugs.freedesktop.org/show_bug.cgi?id=83957 Priority: medium Bug ID: 83957 Assignee: nouveau at lists.freedesktop.org Summary: [NV50 - NVAF][MCP89] Major screen corruption and X crash soon after log in QA Contact: xorg-team at lists.x.org Severity: critical Classification: Unclassified OS:

Full_Name: Madeleine Yeh Version: 1.9.1 OS: AIX 5.2 Submission from: (NULL) (151.121.225.1) After compiling R-1.9.1 on AIX 5.2 using the IBM cc compiler, I ran the checks. One of them failed. Here is the output from running the check solo. root@svweb:/fsapps/test/build/R/1.9.1/R-1.9.1/tests/Examples: ># ../../bin/R --vanilla < stats-Ex.R R : Copyright 2004, The R

Full_Name: Morten Welinder Version: 1.5.1 OS: Solaris Submission from: (NULL) (65.213.85.136) I was having problems with qgamma's precision in the sense that pgamma(qpgamma(x)) was not as close to the identity function as I would like. I was seeing relative errors with random input of about 1e-8. This fits nicely with the code'd EPS2 value of 5e-7. To solve this I added a newton step

Dear R mailing list, Greetings!!! I'm doing a graphical test of suitability for the aircondit data. I actually just followed the exercise in Davison and Hinkley's book of Graphical Test for model checking. The program for testing if it comes from an exponential distribution seems to be running okay but the program for the gamma distribution distribution is running away :-). I hope

Hi there. I am receiving an unexpected error message when creating a new distribution for the survreg() function in the survival package. I understand the survival.distributions() function and have been following the Cauchy example provided in the help file. My goal is to use survreg to fit a gamma distribution to interval censored data. Here is a simple example of what I'm trying to do.