Re: [R] volume of ellipsoid

2013-11-17 Thread yuanzhi
Michael Friendly wrote On 11/14/2013 9:35 AM, yuanzhi wrote: Hi, Carl Witthoft yes, it looks like a mathematical question. I will try based on your suggestion to calculate the volume of the intersection. But I still want to know whether there are some functions in R which can calculate the

Re: [R] volume of ellipsoid

2013-11-15 Thread Michael Friendly
On 11/14/2013 9:35 AM, yuanzhi wrote: Hi, Carl Witthoft yes, it looks like a mathematical question. I will try based on your suggestion to calculate the volume of the intersection. But I still want to know whether there are some functions in R which can calculate the volume of an ellipsoid(area

Re: [R] volume of ellipsoid

2013-11-15 Thread David Winsemius
On Nov 15, 2013, at 8:16 AM, Michael Friendly wrote: On 11/14/2013 9:35 AM, yuanzhi wrote: Hi, Carl Witthoft yes, it looks like a mathematical question. I will try based on your suggestion to calculate the volume of the intersection. But I still want to know whether there are some

Re: [R] volume of ellipsoid

2013-11-14 Thread Carl Witthoft
Well, given that an upvoted question on math.stackexchange got no answers, I'd say you're asking a very difficult question. Perhaps this paper http://www.geometrictools.com/Documentation/IntersectionOfEllipsoids.pdf will be of some help. It's possible that you could do: 1) find the ellipse

Re: [R] volume of ellipsoid

2013-11-14 Thread yuanzhi
Hi, Carl Witthoft yes, it looks like a mathematical question. I will try based on your suggestion to calculate the volume of the intersection. But I still want to know whether there are some functions in R which can calculate the volume of an ellipsoid(area for p=2, hypervolume for p3) containing

Re: [R] volume of ellipsoid

2013-11-14 Thread yuanzhi
Carl Witthoft wrote Well, given that an upvoted question on math.stackexchange got no answers, I'd say you're asking a very difficult question. Perhaps this paper http://www.geometrictools.com/Documentation/IntersectionOfEllipsoids.pdf will be of some help. It's possible that you could

[R] volume of ellipsoid

2013-11-13 Thread Yuanzhi Li
Hi, everyone! I have a matrix X(n*p) which is n samples from a p-dimensional normal distribution. Now I want to make the ellipsoid containing (1-α)% of the probability in the distribution based on Mahalanobis distance: μ:(x-μ)'Σ^(-1)(x-μ)≤χ2p(α) where x and Σ is the mean and