Displaying 20 results from an estimated 800 matches similar to: "hypergeometric series in R"
2008 Dec 03
1
hypergeometric
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
2002 Oct 23
1
vectorizing a function
Dear R-xperts
I have just written a little hypergeometric function, included below
[the hypergeometric function crops up when solving a common type of
ODE].
It works fine on single values of the primary argument z, but
vectorizing it is getting confusing.
The best I have come up with so far just tests for z being longer than
1 and if so, uses sapply() recursively. This is fine, except that it
2005 Sep 19
3
Extended Hypergeometric Distribution
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}}
2007 Feb 20
2
analysis of correlation matrices
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
2012 Mar 19
2
hypergeometric function in ‘ mvtnorm’
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.
2010 Mar 30
1
Multivariate hypergeometric distribution version of phyper()
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.:
2008 Feb 07
1
Appell Hypergeometric function
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
2008 Sep 25
1
What distribution is related to hypergeometric?
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
2003 Oct 01
1
hypergeometric & population estimates
"help"
We want to estimate the number of caribou in Jasper. We recently conducted
an aerial survey and saw 70 uncollared caribou and 8 of 11 collared
caribou. We want to estimate the number of caribou in this population with
95% confidence limits. Gary White uses the hypergeometric distribution and
determines the population estimates using maximum likelihood and 95%CL as
2009 Jun 09
2
Sweave and accents
Hello.
I want to write my notes in Sweave in my own language (spanish). But
my language has accents and when I run Sweave in R to translate my Snw
file into the tex file the accents are translated into unrecognizable
characters.
For example, the word "cami?n" (truck) is translated into "cami??n"
Somebody knows how can I do it?
One solution I don't like is using a
2011 Nov 11
1
Jordan Form of a matrix
Hello.
Is it possible to find the Jordan Form of a matrix with R?
Arnau.
------------------------------------------------------------
Arnau Mir Torres
Edifici A. Turmeda
Campus UIB
Ctra. Valldemossa, km. 7,5
07122 Palma de Mca.
tel: (+34) 971172987
fax: (+34) 971173003
email: arnau.mir at uib.es
URL: http://dmi.uib.es/~arnau
2011 May 05
1
functions pandit and treebase in the package apTreeshape
Hello.
I'm trying to use the functions pandit and treebase. They are in the package apTreeshape. Once I've loaded the package, R responses:
- no function pandit/treebase.
Somebody knows why or what is the reason?
Thanks,
Arnau.
------------------------------------------------------------
Arnau Mir Torres
Edifici A. Turmeda
Campus UIB
Ctra. Valldemossa, km. 7,5
07122 Palma de Mca.
2008 Oct 14
2
help about how can R compute AIC?
Hello.
I need to know how can R compute AIC when I study a regression model?
For example, if I use these data:
growth tannin
1 12 0
2 10 1
3 8 2
4 11 3
5 6 4
6 7 5
7 2 6
8 3 7
9 3 8
and I do
model <- lm (growth ~ tannin)
AIC(model)
R responses:
38.75990
I know the following formula to compute AIC:
AIC=
2010 Jun 02
1
compute the associate vector of distances between leaves in a binary non-rooted tree
Hello.
I'd like to compute the associate vector of distances between leaves in a binary non-rooted tree. The definition of a distance between two leaves in a binary non-rooted tree is the number of edges in the path joining the two leaves.
I've tried the ape package but I'm unable to find this vector.
For example, using rtree(5,rooted=F) I've obtained the following tree:
$edge
1997 May 12
1
R-alpha: Hypergeometric Distribution
A cut and paste typo has crept in and is rendering all values returned
for the hypergeometric distribution incorrect. The problem is in
src/main/arithmetic.c in the function "math4". The lines
PROTECT(sy = allocVector(REALSXP, n));
a = REAL(sa);
b = REAL(sb);
c = REAL(sc);
d = REAL(sc); /* <-- change this line */
y = REAL(sy);
should
2010 Aug 13
1
hypergeometric vs fisher.test
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 ?
What is the difference ?
Thx a lot
--
View this message in context:
2011 Nov 14
2
How to compute eigenvectors and eigenvalues?
Hello.
Consider the following matrix:
mp <- matrix(c(0,1/4,1/4,3/4,0,1/4,1/4,3/4,1/2),3,3,byrow=T)
> mp
[,1] [,2] [,3]
[1,] 0.00 0.25 0.25
[2,] 0.75 0.00 0.25
[3,] 0.25 0.75 0.50
The eigenvectors of the previous matrix are 1, 0.25 and 0.25 and it is not a diagonalizable matrix.
When you try to find the eigenvalues and eigenvectors with R, R responses:
> eigen(mp)
$values
[1]
2012 Oct 01
1
Retrieve hypergeometric results in large scale
I'm going to use
dhyper(x, m, n, k)
to get a 95% coverage. Let me use an example to explain my problem:
Suppose I have a urn containing 90 red and 10 black balls.
Now I wanna remove 3 from the urn. By the following codes:
m<-90;n<-10;k<-3;
x<-0:3
dhyper(x,m,n,k)
I can obtain the probability that 0,1,2,3 red balls will be removed.
0.000742115 0.025046382 0.247680891
2011 Feb 16
1
ExpertFit
Hello.
Somebody knows if there exists some R package that it is capable of doing the same of the software ExpertFit (see http://www.averill-law.com/ExpertFit-distribution-fitting-software.htm)?
Thanks,
Arnau.
------------------------------------------------------------
Arnau Mir Torres
Edifici A. Turmeda
Campus UIB
Ctra. Valldemossa, km. 7,5
07122 Palma de Mca.
tel: (+34) 971172987
fax: (+34)
2000 Mar 24
3
quantiles of the hypergeometric distribution (PR#502)
Hello!
I use R-version 1.0.0
To get the 0.95 quantile of the hypergeometric
distribution with the parameters m=45000,n=5000 and
k=600 I use the R-command
> qhyper(0.95,45000,5000,600).
The value obtained is 600. However, the true value
is 552. The latter can be obtained for example by
calling the corresponding distribution function
with the R commands
> x<-540:580
>