Re: [R] Ellipse in PCA with parameters "a" and "b"defined.

2013-01-22 Thread mary
thank you David, it was my first idea, but i don't know if it is right statistically speaking!!! -- View this message in context: http://r.789695.n4.nabble.com/Ellipse-in-PCA-with-parameters-a-and-b-defined-tp4656215p4656274.html Sent from the R help mailing list archive at Nabble.com. _

Re: [R] Ellipse in PCA with parameters "a" and "b"defined.

2013-01-22 Thread David L Carlson
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of mary > Sent: Tuesday, January 22, 2013 3:40 AM > To: r-help@r-project.org > Subject: Re: [R] Ellipse in PCA with parameters "a" and "b"defined. &

Re: [R] Ellipse in PCA with parameters "a" and "b"defined.

2013-01-22 Thread mary
Ok... so, in my model my "a" is built using the standard deviation of the first principal component and "b" with the second, so my "x" and "Y" should be : PCA $ scores [, 1], PCA $ scores [, 2] but in this way I do not get out a confidence interval set on my parameters but many ellipses. Thanks

Re: [R] Ellipse in PCA with parameters "a" and "b"defined.

2013-01-21 Thread Jim Lemon
On 01/22/2013 11:12 AM, mary wrote: Hi, I have to construct an ellipse interval region on a PCAbiplot, I have my parameters "a" and "b" and I would apply the formula: draw.ellipse(x, y, a = , b = ) I have done a PCA on my data so I have my scores and loading for the first and second component,

[R] Ellipse in PCA with parameters "a" and "b"defined.

2013-01-21 Thread mary
Hi, I have to construct an ellipse interval region on a PCAbiplot, I have my parameters "a" and "b" and I would apply the formula: draw.ellipse(x, y, a = , b = ) I have done a PCA on my data so I have my scores and loading for the first and second component, but my answer is what I have to choose

[R] ellipse

2012-06-06 Thread a.lugu...@gmail.com
Hi, I would like to draw elliptic confidence interval. I know I can use a package ellpise, but can I do it in this way: I want to draw an ellipse countour i.e. all points in R^2 such that t(x*-x)%*%A%*%(x*-x)<=c, for some x* in R^2 and constant c. A is 2x2 matrix. -- View this message in context

Re: [R] Ellipse that Contains 95% of the Observed Data

2010-03-29 Thread Bert Gunter
Typo: "**Paul_i** Exclusion Principle" Bert Gunter Genentech Nonclinical Biostatistics __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and pro

Re: [R] Ellipse that Contains 95% of the Observed Data

2010-03-29 Thread Bert Gunter
oun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Tom La Bone Sent: Monday, March 29, 2010 9:56 AM To: r-help@r-project.org Subject: Re: [R] Ellipse that Contains 95% of the Observed Data I know what "get a bigger sample means". I have no clue what "ask a

Re: [R] Ellipse that Contains 95% of the Observed Data

2010-03-29 Thread Tom La Bone
I know what "get a bigger sample means". I have no clue what "ask a more statistically meaningful question" means. Can you elaborate a bit? Tom -- View this message in context: http://n4.nabble.com/Ellipse-that-Contains-95-of-the-Observed-Data-tp1694538p1695357.html Sent from the R help mailing

Re: [R] Ellipse that Contains 95% of the Observed Data

2010-03-29 Thread Bert Gunter
Easy. See below. Bert Gunter Genentech Nonclinical Biostatistics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Tom La Bone Sent: Monday, March 29, 2010 6:56 AM To: r-help@r-project.org Subject: Re: [R] Ellipse that Contains 95

Re: [R] Ellipse that Contains 95% of the Observed Data

2010-03-29 Thread David Freedman
for a picture of the bagplot, try going to http://www.statmethods.net/graphs/boxplot.html -- View this message in context: http://n4.nabble.com/Ellipse-that-Contains-95-of-the-Observed-Data-tp1694538p1695236.html Sent from the R help mailing list archive at Nabble.com. _

Re: [R] Ellipse that Contains 95% of the Observed Data

2010-03-29 Thread Tom La Bone
Concisely, here is what I am trying to do: #I take a random sample of 300 measurements. After I have the measurements #I post stratify them to 80 type A measurements and 220 type B measurements. #These measurements tend to be lognormally distributed so I fit them to #determine the geometric me

Re: [R] Ellipse that Contains 95% of the Observed Data

2010-03-29 Thread S Ellison
The bagplot at http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=112 gives a nonparametric 2-d view analagous to a boxplot. S Ellison > > I can take the results of a simulation with one random variable and generate > an empirical interval that contains 95% of the observations, e.g.

Re: [R] Ellipse that Contains 95% of the Observed Data

2010-03-29 Thread Jim Lemon
On 03/29/2010 07:17 PM, Barry Rowlingson wrote: ... I think the problem as posed doesn't produce a unique ellipse. You could start with a circle of radius 0 centered on mean(x),mean(y) and then increase the radius until it has 95% of the points in it. As long as your points are in continuous sp

Re: [R] Ellipse that Contains 95% of the Observed Data

2010-03-29 Thread Barry Rowlingson
On Mon, Mar 29, 2010 at 4:02 AM, Ben Bolker wrote: >  I'll be interested to hear what others come up with. >  I'm not sure the problem as you have stated it is well-posed, or > necessarily possible. Suppose there is a true unknown > bivariate probability distribution with a non-elliptical 95% > q

Re: [R] Ellipse that Contains 95% of the Observed Data

2010-03-28 Thread Ben Bolker
Tom La Bone gforcecable.com> writes: > > > I can take the results of a simulation with one random variable and generate > an empirical interval that contains 95% of the observations, e.g., > > x <- rnorm(1) > quantile(x,probs=c(0.025,0.975)) > > Is there an R function that can take the re

[R] Ellipse that Contains 95% of the Observed Data

2010-03-28 Thread Tom La Bone
I can take the results of a simulation with one random variable and generate an empirical interval that contains 95% of the observations, e.g., x <- rnorm(1) quantile(x,probs=c(0.025,0.975)) Is there an R function that can take the results from two random variables and generate an empirical

Re: [R] Ellipse: Major Axis Minor Axis

2009-09-08 Thread Jari Oksanen
Vishal gmail.com> writes: > > Jari, thanks for the quick answer. > > > sqrt(eigen(cov.trob(mydataforellipse)$cov)$values) > > what will this return? > > For my data, I get: > > > sqrt(eigen(cov.trob(r)$cov)$values) > [1] 1.857733e-05 4.953181e-06 > > Is this Left hand value the major or the

Re: [R] Ellipse: Major Axis Minor Axis

2009-09-02 Thread Vishal
I forgot to mention that the ellipse has a rotation. It's horizontal axis is not parallel to the x-axis. Just repeating, when I try to draw a circuit keeping the center same as the ellipse, and radius equal to the first value returned by > sqrt(eigen(cov.trob(mydataforellipse)$cov)$values) taking

Re: [R] Ellipse: Major Axis Minor Axis

2009-09-02 Thread Vishal
Jari, thanks for the quick answer. > sqrt(eigen(cov.trob(mydataforellipse)$cov)$values) what will this return? For my data, I get: > sqrt(eigen(cov.trob(r)$cov)$values) [1] 1.857733e-05 4.953181e-06 Is this Left hand value the major or the semi major length? I also try to plot a circuit keepi

Re: [R] Ellipse: Major Axis Minor Axis

2009-09-01 Thread Jari Oksanen
Vishal gmail.com> writes: > > I have a matrix(3000x2) of numbers and I have plotted a scatterplot > (defined in the ``car'' library.) > > scatterplot(r$V1,r$V2,ellipse=TRUE) > > The ellipse plotted is an error ellipse. > > I want the find the major(semi), minor(semi) minor axis length of this

[R] Ellipse: Major Axis Minor Axis

2009-09-01 Thread Vishal
I have a matrix(3000x2) of numbers and I have plotted a scatterplot (defined in the ``car'' library.) scatterplot(r$V1,r$V2,ellipse=TRUE) The ellipse plotted is an error ellipse. I want the find the major(semi), minor(semi) minor axis length of this ellipse. Does anybody know how to do it? I ca

Re: [R] ellipse

2009-05-05 Thread Duncan Murdoch
On 05/05/2009 6:28 PM, Antonio Lucadamo wrote: Dear all, I'm using the ellipse package and I would like to verify if the confidence region that I build with this package can be equivalent to an Union Intersection Test. I used different value for the t-statistic but I can not find the right eq

[R] ellipse

2009-05-05 Thread Antonio Lucadamo
Dear all, I'm using the ellipse package and I would like to verify if the confidence region that I build with this package can be equivalent to an Union Intersection Test. I used different value for the t-statistic but I can not find the right equivalence. Does someone know how to choose the

[R] ellipse

2009-05-04 Thread Antonio Lucadamo
Dear all, I'm using the ellipse package and I would like to verify if the confidence region that I build with this package can be equivalent to an Union Intersection Test. I used different value for the t-statistic but I can not find the right equivalence. Does someone know how to choose the

Re: [R] ellipse in pairs plot

2008-11-14 Thread Kaveh Vakili
Dear list moderator, I'm sending you a solution i found to the question i asked the list yesterday (i have not been able to figure out how to post it to the list directly). elly<-function(w,col=array(1,dim(w)[1]),pch=array(1,dim(w)[1])){ p<-dim(w)[2] op<-par(mfrow=c(3,3),mar=c

[R] ellipse in pairs plot

2008-11-13 Thread Kaveh Vakili
Hi dear list, I'm scrambling with this: i have a p-dimensional data frame for which i want to have a pairs() plot, with ellipse added in each lower.panel. I want the ellipse's parameter (mu, sigma) to be extracted from the appropriate entries in a pre-computed matrix ("dab" in the code below)