similar to: Draw overlap ellipse in R

library (ellipse) shape1 = c (1, 0, 0,1) dim(shape1) = c(2,2) ellipse (center = c(0,0), shape = shape1, radius = 1) ============================= Error in plot.xy(xy.coords(x, y), type = type, col = col, lty = lty, ...) : plot.new has not been called yet It is really frustrating. Also what do the shape matrix, radius correspond to an ellipse function (x-x0)^2/a + (y-y0)^2/b = 1 ?

Hi, I drew an ellipse with the package ellipse. Now I would like to know if a point is inside the ellipse. Is any R functions to do it without computing the equation of the ellipse manually? Thanks. For example, if I do "plot(ellipse(0.8), type = 'l')", I would like to know if (0,1) belongs to the drawn ellipse. Regards, Alain -- Alain Guillet Statistician and Computer

Dear S-PLUS/R users: Do you know any default function or a user contributed function that can draw an ellipse with given axes and origin? Thanks for any help. Paul. -- Romance, like alcohol, should be enjoyed, but should not be allowed to become necessary. -- Edgar Friedenberg -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list --

Hello everyone I really need some help here. I made a confidence ellipse using the function ellipse from the package ellipse: ellipse(SD, centre=colMeans(pcsref),t=sqrt((p * (n-1)/(n-p))*qf(0.99, p,n-p)) Now, I want to write a function whom return TRUE or FALSE if a given observation is in the confidence ellipse. But I have no clue how to do it Can anyone help me? Best regards Jessica

Hello again, First I want to thank all the people who answered my question about line width in graphs. I promise I will learn the 'par' help page by heart for the end of the month ! I now want to trace some ellipses to emphasize groups of data. I found how to trace circles with 'symbols()', but no ellipse. I'm planning on writing my own function based on

That's probably a stupid question, but I'm looking for a low-level command which plots ellipse, specifying only center and deformation axes. The purpose is to illustrate bivariates gaussians with 2D .95 confident levels without using any specific library. Thanxs for your help, Guillaume

Hi: Does Ellipse and dataellipse function in R produce the same ellipse? I wanted to see how the radius for the Ellipse function in R calculated. Also what is the var-covariance matrix, if any, assumed for the dataellipse function? Heres an example of the code where I am generating Multivariate normal data and creating ellipse using the 2 functions: library(car) library(mvtnorm) mu =

Hi, I'm using the R library(car) to draw confidence/prediction ellipses in a scatterplot. >From what i understood the ellipse() function return an ellipse based parameters: shape, center, radius . If i read dataEllipse() function i can see how these parameters are calculated for a confidence ellipse. ibrary(car) a=c(12,12,4,5,63,63,23) b=c(13,15,7,10,73,83,43) v <-

Hello, SAS' point and click interface has the option of produce a scatterplot with a superimposed confidence ellipse for the correlation coefficient. Since I generally like R so much better, I would like to reproduce this in R. I've been playing with the ellipse package. In order to have the points and the ellipse on the same graph I've done the following. (Load ellipse

Hi, I created a principal component plot using the first two principal components. I used the function princomp() to calculate the scores. now, I would like to superimpose an ellipse representing the center and the 95% confidence interval of a series of points in my plot (as to illustrate the grouping of my samples). I looked at the ellipse() function in the ellipse package but can't get it

I have been using the ellipse function from the car package and the covariance matrix to draw an ellipse around a group of points to show the confidence limits. However, the points are actually represented by 3 variables so rather than plot each pairwise combination of variables in 2D I would like to plot the 'ellipse' in 3D using the djmrgl package. Can anyone offer advice on how I can

Hi all Arnaud, Some newly purchased Ellipse ASR-model UPS find their way to my office for a first evaluation test. It concerns the Ellipse 1500 and Ellipse 1000 models. As you know we are still working with hidups. Peripheral recognition, start of NUT, mains short interruptions management, all that work fine... as expected ;-) But changes occured on the delayed shutdown procedure after the

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(10000) 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 ellipse that contains 95% of the observations, e.g., x <- rnorm(10000) y <- rnorm(10000) ?

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 as X and Y into the formula? if "a" and "b" are scalars or

Dear list members, This might be a silly question but I just can't figure it out. I am using the "ellipse" package on covariance matrices. I would simply like to plot my ellipses WITH its two axis ploted as well. These axis represents the 2 eigen vectors of my matrix and it is important that I can graphically show them. Is there an easy way to do so? Many thanks, Antoine -- View

Full_Name: Volker Franz Version: Version 1.6.2 (2003-01-10) OS: Debian Submission from: (NULL) (192.124.28.104) Hi there, it seems to me that data.ellipse of package "car" (Version 1.0-1) produces confidence interval's which are too big. To see this, do: library(car) plot(c(-2,2),c(-2,2),pch=0) data.ellipse(rnorm(10000),rnorm(10000),levels=0.68,plot.points=F) abline(v=+1)

Dear All, I would be very appreciative of your help with the following 1). I am running multivariate multiple regression through the manova() function (kindly suggested by Professor Venables) and getting two different answers for test=c("Wilks","Roy","Pillai") and tests=c("Wilks","Roy",'"Pillai") as shown below. In the

Hi everyone,I need you all help. I want to create few image in same diagram. For example, I wan draw two ellipse in same diagram. plot(ellipsePoints(2,5, alpha = 60), asp=1) and plot(ellipsePoints(2,5, alpha = 30), asp=1) How to display both ellipse together? both is different in angle. Thank you. -- Hemavathi Ramulu [[alternative HTML version deleted]]

Hello. So as my english is not so good, i'll try to write my story in bash commands :o) miwa@zork3$ uname -mrs FreeBSD 5.5-STABLE i386 miwa@zork3$ su root@zork3# cd /usr/ports/sysutils/nat root@zork3# make && make install root@zork3# cd ../nut-usb root@zork3# make && make install root@zork3# cd /usr/local/etc/nut root@zork3# cat ups.conf user=root [MGE-Ellipse]

My UPS is a Eaton Ellipse 1000S (But its identified as MGE Ellipse 0463:ffff). I have been using the usb interface with the usbhid-ups driver, but I found that sometimes (too much times) the driver doesn't respond and NUT doesn't receive changes nor notifications from the UPS. The kernel lose the connections with USB device, in the kernel logs appears: Dec 30 11:24:12 rei kernel: