similar to: weird error from MASS::eqcsplot with postscript driver

I found the following in eqscplot, library MASS: When given just one point, it fails: > eqscplot(x=1, y=1) Error in plot.window(xlim, ylim, log, asp, ...) : need finite xlim values After inserting browser() just before the last line in eqscplot() (in which plot() gets called), I found the following: > eqscplot(x=1, y=1) Called from: eqscplot(x = 1, y = 1) Browse[1]> xlim [1] NaN

This is not a bug. There is no way one can scale the plot axes equally when no indication has been given of the scaling expected for either! It has always seemd to be that plot(x=1, y=1) shoul dnot work: it has not supplied enough information! On Thu, 8 Nov 2001 pflugshaupt@geobot.umnw.ethz.ch wrote: > I found the following in eqscplot, library MASS: > > > When given just one

Dear all I was writing some code that needed functions from packages 'MASS' and 'CircStats', and I received a warning saying that the function 'eqscplot' in one of the packages was masked by another 'eqscplot' from the other package (in fact 'eqscplot' from 'CircStats' seems a short version of 'eqscplot' from 'MASS'). This masking of

Hi !! I am trying to analyze some of my data using linear discriminant analysis. I worked out the following example code in Venables and Ripley It does not seem to be happy with it. ============================ library(MASS) library(stats) data(iris3) ir<-rbind(iris3[,,1],iris3[,,2],iris3[,,3]) ir.species<-factor(c(rep("s",50),rep("c",50),rep("v",50)))

Please, could someone help me figure out what seems to be a very simple problem (and is still taking me hours...). I want to draw a simple scatterplot but with 'equal' axes, i.e. I want both axes to go from -3 to 3. Values for x lie between -2 and 0.5, values for y between -2.2 and 3. I have tried 'usr' and 'eqscplot' and a few other options, but it doesn't give me the

Hi, I am using the lda function from the MASS library. I would to find the decision boundaries of each class and subsequently plot them. I wonder if anybody can offer any help on this topic? Below I applied the lda function on a small dataset of mine. Any help will be much appreciated. > library(MASS) >

Full_Name: Laurent Gautier Version: 1.3.0-patched OS: IRIX 6.5 Submission from: (NULL) (130.225.67.199) Hello, The function La.eigen, called by cmdscale in the package mva behaves an unexplicable way (for me). The following lines show what happened. I tried the very same on linux, and it worked fine. >a <- matrix(c(1,2,3,2),3,3) >a [,1] [,2] [,3] [1,] 1 2 3 [2,]

Dear R-experts, Searching the help archives I found a recommendation to do multivariate outlier identification by mahalanobis distances based on a robustly estimated covariance matrix and compare the resulting distances to a chi^2-distribution with p (number of your variables) degrees of freedom. I understand that compared to euclidean distances this has the advantage of being scale-invariant.

I've got a small problem. I have some observational data (environmental samples: abiotic explanatory variable and biological response) to which I've fitted both a multiple linear regression model and also a gam (mgcv) using smooths for each term. The gam clearly fits far better than the lm model based on AIC (difference in AIC ~ 8), in addition the adjusted R squared for the gam is

Dear useRs, I'd like to produce some scatter plots where N units on the X axis are equal to N units on the Y axis (as measured with a ruler, on screen or paper). This approach x <- sample(10:200,40) ; y <- sample(20:100,40) windows(width=max(x),height=max(y)) plot(x,y) is better than plot(x,y) but doesn't solve the problem because of the other parameters (margins etc). Is

Greetings, I am working on a stylometric analysis of some latin texts; one of the latest stylometric techniques involves using principal components analysis. Not being a statistician, I can't really fully rely on PCA as my primary tool, since I don't really understand the statistics behind the PCA technique. Nevertheless, the ability to use PCA and graph the results has been

Last week I asked about data ellipses with rgl:::ellipse3d() with lines showing the principal axes. (The goal is a visual demonstration of PCA as a rotation of variable space to component space.) I was trying, unsuccessfully, to use princomp() to generate the PCA axes and plot them using segments3d: > > PC <- princomp(trees) > > sdev <- PC$sdev # component standard

Hi, I have a question about computing shortest Euclidean distances between two data frames of spatial points... I have 2 dataframes (not yet converted to spatial class) >Sewers<-data.frame(x=c(10,20,35,50),y=c(100,150,200,300)) >transect <- data.frame(x=seq(from=0, to=50, by=1),y=seq(from=100, to=150, by=1)) I would like to be able to compute the shortest distance from the

Hello: What I want to do is quite simple, but I can't find a way. I have a data frame with several points (x and y coords). I want to add another column with cluster membership. For example aggregate all the points that stand within a distance of 40 from each other. I've tried using "nncluster" from the package nnclust, but the results are not correct, for some

I want to plot something (eg a circle) with a fixed ratio of the x and y axis, or (even better) with a fixed size when I print it. Output should then be a circle (actually it'll be someting more complicated) with radius 5cm and not an ellipse. I'm _sure_ this is not new, but after looking 45min for a solution, I post here... Thanks for help Thomas

Am Dienstag, 25. April 2017, 11:21:31 schrieb Dirk Eddelbuettel: > On 25 April 2017 at 16:11, Johannes Ranke wrote: > | This looks similar to what I got this morning when I tested my > | (unreleased) > | backport of R 3.4.0 to Debian jessie. My test was > | > | library(MASS) > | example(rlm) > | > | and there was an object that was not found. I am on a train on the way

Greetings. I'm about to embark on my first big (to me at least!) R project, which will be to write a function to plot three-dimensional multiple correspondence analysis (mca) plots in a manner similar to scatterplot3d(). (plot.mca() plots only two dimensions, even though mca() will calculate more.) Before I do so, however, I would love to know that I'm not reinventing the wheel or any

Hi, Just downloaded the latest cvs from zaptel on my sparking new Athlon64 Centos4.2 system, but hitting a stumbling block... (sorry for the long post) #make cc -I. -O4 -g -Wall -DBUILDING_TONEZONE -m64 -DSTANDALONE_ZAPATA -DZAPTEL_CONFIG=\"/etc/zaptel.conf\" -c -o gendigits.o gendigits.c cc -o gendigits gendigits.o -lm ./gendigits cc -I. -O4 -g -Wall -DBUILDING_TONEZONE -m64

Hi, I was wondering if anyone knows of a way by which one can estimate values from a contour plot created by using the loess function? I am hoping to use the loess contour plot as a means of interpolation to identify the loess created values at points at pre-defined (x,y) locations. Could anyone point me in the right direction please? Thanks. Laura Quinn Institute of Atmospheric Science School

Hello everybody, I am doing a plot and I would like to know which parameter I have to use in order to get an equal scale for the axes x and y: Means, if the intervals in the x-axis is 0,5,10,15.... I would like to have the same scale for the y-axis. Thank you so much for your help. Halima from Leiden where the sun is shining. Holland