Displaying 20 results from an estimated 400 matches similar to: "hypergeometric function in ‘ mvtnorm’"
2008 Mar 22
1
Vectorization Problem
I have the code for the bivariate Gaussian copula. It is written with
for-loops, it works, but I wonder if there is a way to vectorize the
function.
I don't see how outer() can be used in this case, but maybe one can
use mapply() or Vectorize() in some way? Could anyone help me, please?
## Density of Gauss Copula
rho <- 0.5 #corr
R <- rbind(c(1,rho),c(rho,1)) #vcov matrix
id <-
2000 Nov 14
1
mvtnorm
Announcement: mvtnorm
Multivariate Normal and T Distribution
mvtnorm implements two R functions for the computation of the multivariate t
and normal distribution:
pmvt: Computes the the distribution function of the multivariate t
distribution for arbitary limits, degrees of freedom and
correlation matrices based on algorithms by Genz and Bretz.
pmvnorm: Computes the distribution
2000 Nov 14
1
mvtnorm
Announcement: mvtnorm
Multivariate Normal and T Distribution
mvtnorm implements two R functions for the computation of the multivariate t
and normal distribution:
pmvt: Computes the the distribution function of the multivariate t
distribution for arbitary limits, degrees of freedom and
correlation matrices based on algorithms by Genz and Bretz.
pmvnorm: Computes the distribution
2004 May 20
4
pmvt problem in multcomp
Hi, all:
Two examples are shown below.
I want to use the multiple comparison of Dunnett.
It succeeded in upper case "example 1".
However, the lower case "example 2" went wrong.
In "example 2", the function pmvt return NaN, so I cannot show
this simtest result. Is there any solution?
(I changed the variable "maxpts" to a large number in front of
the
2007 May 09
2
pvmnorm, error message
Hello there!
My operating system is Windows XP, my version of R is the latest (R-2.5.0). Recently I have downloaded the package "mvtnorm" and a problem with the command "pmvnorm" occured. Trying to enter the lines ...
A <- diag(3)
A[1,2] <-0.5
A[1,3] <- 0.25
A[2,3] <- 0.5
pvmnorm(lower=c(-Inf,-Inf,-Inf), upper=c(2,2,2),mean = c(0,0,0), corr=A)
I got the following
2005 Feb 24
1
Density of the Multivariate T Distribution
Hi,
I am looking for an efficient way to compute the
values of the density function of a multivariate T
distribution - something like "dmvnorm", but for T
distr. Does this exist somewhere?
Many thanks,
Jan Bulla
Goettingen University
2002 Jul 09
1
RE: mvtnorm package installation failure
Hi,
Thank you for the tip.
I tried to re-intall R from Debian "stable", in which R's version is 1.4.0. And the installation
of "mvtnorm" works.
I then re-installed R yet again from Debian "unstable" (woody), in which R's version is 1.5.1.
The installation of "mvtnorm" fails again with the same error message. Another package that
failed with the
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}}
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
2001 Sep 08
0
R-function available for noncentral hypergeometric distribution
For those who are interested, I have made available a R function for
noncentral hypergeometric distribution at
http://www.geocities.com/jg_liao/software/Hypergeometric/hypergeometric_in_R.txt
The paper that describes the algorithm will appear in The American
Statistician.
The function does not run on S-plus as the R's scoping rule is used.
Here is how the function can be used:
> n1
2010 Jun 08
1
hypergeometric series in R
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.
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
2011 Aug 31
0
Fitting the negative hypergeometric distribution
I'd like to fit the (1-displaced) negative hypergeometric distribution
to data samples such as the following:-
x <- c(370, 311, 299, 266, 265, 232, 197, 198, 170, 154, 133, 123, 120,
103, 80, 72, 69, 67, 67, 50, 36, 35, 26, 23, 15, 11, 9, 6, 5, 3, 3, 2,
2, 2)
i.e., I want to estimate the parameter values of K and M (with my data,
n would usually be the same as the length of the data
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:
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.:
2003 May 02
1
Does R implement Hypergeometric test?
[This email is either empty or too large to be displayed at this time]
2008 Jan 15
0
FDR for hypergeometric tests
Dear list,
I have performed several tests for the hypergeometric distribution
using phyper() for some gene annotation categories as follows
>phyper(26,830,31042,337, lower.tail=F)
>phyper(16,387,31042,337, lower.tail=F)
.
.
.
I am only running some selected categories but I would like to correct
this value for multiple testing since I
have 3121 possible tests according to 3121
1998 Feb 23
0
R-beta: Hypergeometric Probabilities
In both versions of R to which I currently have access (R-0.16.1 and
R-0.61.1), "phyper" stops returning correct cumulative probabilities as the
parameters of the hypergeometric distribution get large. For example, when
N1=1345, N2=1055, and n=1330, phyper returns either 0 or 1, and nothing in
between.
Looking at phyper.c, it's clear what's happening. First a term (called