Re: [R] Hypergeometric Function seems to give wrong results

2015-07-08 Thread Ben Bolker
Carlos Nasher via R-help r-project.org> writes: > [snip] > I need to evaluate the Hypergeometric Function of the 2nd kind (Tricomi > confluent hypergeometric function). Therefore I'm using the kummerU > function from the fAsianOptions package. It seems to me that kummerU is > giving wrong resul

[R] Hypergeometric Function seems to give wrong results

2015-07-06 Thread Carlos Nasher via R-help
Hello R helpers, I need to evaluate the Hypergeometric Function of the 2nd kind (Tricomi confluent hypergeometric function). Therefore I'm using the kummerU function from the fAsianOptions package. It seems to me that kummerU is giving wrong results. Here's an example: library("fAsianOptions") ku

Re: [R] hypergeometric integral

2014-06-17 Thread JLucke
r=0,upper=.51,H=0.7) Error in integrate(Gauss2F1, lower = 0, upper = 0.51, H = 0.7) : non-finite function value > jlu...@ria.buffalo.edu Sent by: r-help-boun...@r-project.org 06/17/2014 10:24 AM To Xuse Chuse , cc r-help , r-help-boun...@r-project.org Subject Re: [R] hypergeomet

Re: [R] hypergeometric integral

2014-06-17 Thread JLucke
1.028055 1.019653 1.013143 1.007920 1.003618 [11] 1.00 #whereupon you find your second problem. #Joe Xuse Chuse Sent by: r-help-boun...@r-project.org 06/17/2014 08:32 AM To r-help , cc Subject [R] hypergeometric integral Dear all, I am trying to make a numerical integral of a Hyp

[R] hypergeometric integral

2014-06-17 Thread Xuse Chuse
Dear all, I am trying to make a numerical integral of a Hypergeometrical function 2F1 but it is not working. I dont know how to pass the arguments to my function. Thank you beforehand. Chuse. Here is my code: H <- 0.7; t <- 1; s <- seq(0,t,0.1); x <- 1-t/s; library(gsl) Gauss2F1 <- function(H,

Re: [R] hypergeometric function in ‘ mvtnorm’

2012-03-19 Thread Berend Hasselman
On 19-03-2012, at 16:54, statfan wrote: > Thanks for your advice. I actually meant to ask about the "pmvt" for the > distribution function. Viewing the source code "pmvt" uses the function > "mvt" which uses the function "probval" which sources the fortran code: > No it doesn't "source". It

Re: [R] hypergeometric function in ‘ mvtnorm’

2012-03-19 Thread statfan
Thanks for your advice. I actually meant to ask about the "pmvt" for the distribution function. Viewing the source code "pmvt" uses the function "mvt" which uses the function "probval" which sources the fortran code: Fortran("mvtdst", N = as.integer(n), NU = as.integer(df), LOWER = as.d

Re: [R] hypergeometric function in ‘ mvtnorm’

2012-03-19 Thread R. Michael Weylandt
To view the source of (most) functions, simply type funcname without parentheses: here, you get > dmvt function (x, delta, sigma, df = 1, log = TRUE, type = "shifted") { if (df == 0) return(dmvnorm(x, mean = delta, sigma = sigma, log = log)) if (is.vector(x)) { x <- matrix

[R] hypergeometric function in ‘ mvtnorm’

2012-03-18 Thread statfan
Is there any way to know how the "dmvt" function computes the hypergeometric function needed in the calculation for the density of multivariate t distribution? -- View this message in context: http://r.789695.n4.nabble.com/hypergeometric-function-in-mvtnorm-tp4483730p4483730.html Sent from the

Re: [R] hypergeometric vs fisher.test

2010-08-13 Thread TGS
"or an incompetent" Such harsh words Peter. You're not making this a friendly environment for people to ask questions. Is there a less competent attribute mailing list so that some of us don't offend you with our questions? On Aug 13, 2010, at 2:11 PM, Peter Dalgaard wrote: Andrea Franceschini

Re: [R] hypergeometric vs fisher.test

2010-08-13 Thread Peter Dalgaard
Andrea Franceschini wrote: > I ask the question also because I found this line in Wikipedia: > "The test (see above) based on the hypergeometric distribution > (hypergeometric test) is identical to the corresponding one-tailed > version of Fisher's exact test". > > Is this wrong ? May I kindly as

Re: [R] hypergeometric vs fisher.test

2010-08-13 Thread Andrea Franceschini
I ask the question also because I found this line in Wikipedia: "The test (see above) based on the hypergeometric distribution (hypergeometric test) is identical to the corresponding one-tailed version of Fisher's exact test". Is this wrong ? May I kindly ask a friendly explanation for not-exper

Re: [R] hypergeometric vs fisher.test

2010-08-13 Thread peter dalgaard
On Aug 13, 2010, at 3:47 PM, Andrea Franceschini wrote: > > Dear R team, > I have a simple question. > > I tried this command: > phyper(17,449,19551,181, FALSE) > [1] 1.47295e-07 > > and then I tried this command: > (fisher.test(matrix(c(17,449,181,19551),2,2), > alternative='greater'))$p.valu

[R] hypergeometric vs fisher.test

2010-08-13 Thread Andrea Franceschini
Dear R team, I have a simple question. I tried this command: phyper(17,449,19551,181, FALSE) [1] 1.47295e-07 and then I tried this command: (fisher.test(matrix(c(17,449,181,19551),2,2), alternative='greater'))$p.value [1] 3.693347e-06 Shouldn't be identical the results of the two commands ? Wh

Re: [R] hypergeometric series in R

2010-06-08 Thread Ben Bolker
Arnau Mir uib.es> writes: > > Hello. > > Somebody knows how to compute generalized hypergeometric series in R? > (see > http://functions.wolfram.com/HypergeometricFunctions/HypergeometricPFQ/02/ > to understand what I mean) library(sos) findFn("generalized hypergeometric function") -> see

[R] hypergeometric series in R

2010-06-08 Thread Arnau Mir
Hello. Somebody knows how to compute generalized hypergeometric series in R? (see http://functions.wolfram.com/HypergeometricFunctions/HypergeometricPFQ/02/ to understand what I mean) Thanks in advance, Arnau. __ R-help@r-project.org mailing list

Re: [R] hypergeometric

2008-12-03 Thread Robin Hankin
The hypergeo package should be able to deal with this, although the function you specify below looks like a degenerate case (if I understand it correctly) so the convergence rate is likely to be slow. Let me know how you get on best wishes Robin (author of hypergeo) Jarle Brinchmann wrote:

Re: [R] hypergeometric

2008-12-03 Thread Jarle Brinchmann
The dhyper etc deal with the hypergeometric _distribution_ while what you appear to want have is the hypergeometric special function (the connection is that the regular hypergeometric function is the generating function for the hypergeometric distribution if I recall correctly). Anyway, what you

[R] hypergeometric

2008-12-03 Thread Zakaria, Roslinazairimah - zakry001
Hi, I hope somebody can help me on how to use the hypergeometric function. I did read through the R documentation on hypergeometric but not really sure what it means. I would like to evaluate the hypergeometric function as follows: F((2*alpha+1)/2, (2*alpha+2)/2 , alpha+1/2, betasq/etasq).

[R] hypergeometric function

2008-12-03 Thread Roslina Zakaria
Hi,   I hope somebody can help me on how to use the hypergeometric function.  I did read through the R documentation on hypergeometric but not really sure what it means.   I would like to evaluate the hypergeometric function as follows: F((2*alpha+1)/2, (2*alpha+2)/2 , alpha+1/2, betasq/etasq).