similar to: vectorizing a function

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). I'm not sure which function should be used- either phyper or qhyper or dhyper Where

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.

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 R help mailing list archive at Nabble.com.

Dear All, I am looking for an implementation in R of the Appell Hypergeometric function. Any suggestions will be more than appreciated! GP -- dr. Giovanni Parrinello External Lecturer Medical Statistics Unit Department of Biomedical Sciences Viale Europa, 11 - 25123 Brescia Italy Tel: +390303717528 Fax: +390303717488 email: parrinel at med.unibs.it

hi ?integer says: Note that on almost all implementations of R the range of representable integers is restricted to about +/-2*10^9: 'double's can hold much larger integers exactly. I am getting very confused as to when to use integers and when not to. In my line I need exact comparisons of large integer-valued arrays, so I often use as.integer(), but the above

Hello everyone I am trying to get one of my packages through R's QC. The package is clean for me under R-2.4.1, R-2.5.0, and R-devel, but Kurt gets > > * checking R code for possible problems ... WARNING > hypercube: no visible binding for global variable ?f? Function hypercube() [cut-&-pasted below] is intended to return an adjacency matrix for an n-dimensional

Hi R-users,   I have this code here: library(numDeriv)   fprime <- function(z) { alp  <- 2.0165;   rho  <- 0.868;   # simplified expressions   a      <- alp-0.5   c1     <- sqrt(pi)/(gamma(alp)*(1-rho)^alp)   c2     <- sqrt(rho)/(1-rho)   t1     <- exp(-z/(1-rho))   t2     <- (z/(2*c2))^a   bes1   <- besselI(z*c2,a)   t1bes1 <- t1*bes1   c1*t1bes1*t2 }   ## Newton

Jacobi's theta functions crop up here and there in various branches of mathematics, such as unimodular and elliptic functions. They have Taylor expansions, but the powers increase as the square of n, as in 1 + z + z^4 + z^9 + z^16 + z^25 + . . . so they converge very quickly, provided |z|<1 The following toy function shows how I'm implementing these objects. I just add terms

Dear R Users, There exists a non-central hypergeometric distribution function in the (MCMCpack) package, and a hypergeometric distribution function in the (stats) package. Is there a function for sampling from an extended hypergeometric distribution? Thanks, Narcyz This message is intended for the addressee named and may con...{{dropped}}

Hello, I'm looking for a package in R that performs, analysis of correlation matrices: cross-classified by 2 factors. The orginal reference to this method is by CJ Brien Biometrica (1998) 75(3):469-76. THank you, Agnes ------------------------------------- E. Agnes Richards, Ph.D. Candidate Department of Zoology University of Toronto at Mississauga 3359 Mississauga Rd. North

Hi everyone. I have a problem that I have been unable to determine either the best way to proceed and why the methods I'm trying to use sometimes fail. I'm using the pf() function in an optimization function to find a noncentrality parameter that leads to a specific value at a specified quantile. My goal is to have a general function that returns the noncentrality parameter that

Dear R Users, I employed the phyper() function to estimate the likelihood that the number of genes overlapping between 2 different lists of genes is due to chance. This appears to work appropriately. Now i want to try this with 3 lists of genes which phyper() does not appear to support. Some googling suggests i can utilize the Multivariate hypergeometric distribution to achieve this. eg.:

R Help on 'chisq.test' states that "if 'simulate.p.value' is 'TRUE', the p-value is computed by Monte Carlo simulation with 'B' replicates. In the contingency table case this is done by random sampling from the set of all contingency tables with given marginals, and works only if the marginals are positive... In the

R-users E-mail: r-help@r-project.org I have a quenstion on "gam()" in "gam" package. The help of gam() says: 'gam' uses the _backfitting algorithm_ to combine different smoothing or fitting methods. On the other hand, lm.wfit(), which is a routine of gam.fit() contains: z <- .Fortran("dqrls", qr = x * wts, n = n, p = p, y = y *

I have been reading, in various sources, that a poisson distribution is related to binomial, extending the idea to include numbers of events in a given period of time. In my case, the hypergeometric distribution seems more appropriate, but I need a temporal dimension to the distribution. I have weekly samples of two kinds of events: call them A and B. I have a count of A events. These change

Dear all, Can anyone please help me about any of the following questions: -------------------------------------------- 1. How can i find "factorial" of any number in R? I tried > prod(170:1) # to find factorial of 170 or 170! Is it the only procedure - or R has any better process / operational character to calculate factorial? Also, is it possible to calculate factorial of 500? Or

On Fri, 2017-01-06 at 12:54 -0500, m.roth at 5-cent.us wrote: > James B. Byrne wrote: > > > > On Thu, January 5, 2017 17:23, Always Learning wrote: > >> > >> Cyber attacks are gradually replacing armed conflicts. > > > > Better fight with bits than blood. > > Yes, but... attacks on the friggin' IoT could result in lots of blood. Or, >

Dear XpeRts, I prepared a no qoute Character string by the following command s<-noquote(paste (b1, collapse=",")) where, b1 is the vector of 24 intergers. > dput(b1) c(1L, 2L, 6L, 7L, 12L, 16L, 17L, 20L, 21L, 23L, 25L, 34L, 46L, 48L, 58L, 64L, 65L, 68L, 82L, 97L, 98L, 101L, 113L, 115L) > dput(s)

I have a question. I can find a function about "Gaussian hypergeometric function" in sofeware Maple V. I need this function to calculate. Is there a similar function in R ? Or I have to design such a function. Thanks! Jimmy Cheng

Hello This is a simple question but I couldn't google an answer. In the procrustes function of the vegan package, one uses plot(procrustes_object, kind=2) to obtain a plot of the residual differences. For instance: data(varespec) vare.dist <- vegdist(wisconsin(varespec)) library(MASS) mds.null <- isoMDS(vare.dist, tol=1e-7) mds.alt <- isoMDS(vare.dist,