similar to: setting up a plot without plotting any data

Hi all: following the previous discussion, it looks like plot(x) with y=NULL still does not work correctly. If one tries for example plot(1:5) it works, but already for plot(runif(100)) it does not. I posted the proposed correction for plot.POSIXct and plot.POSIXlt before. Please voice your opinions whether the following fix for plot.default could be reasonable? I include the full function and

There seems to be a minor bug in barplot.default when used with log scale w= here one or more values is NA: dat <- matrix(1:25, 5) dat[2,3] <- NA barplot(dat, beside =3D T) #Plots and appropriate barplot with gaps for m= issing data barplot(dat, beside =3D T, log =3D "y") #Error in if (min(height + offset) <=3D 0) stop("log scale error: at least = one 'height +

Hi! The other day I asked a couple of questions about plotting some data. Firstly I wanted to specify the colours used for the text plotted using: > text(chem.predict[,1:2], labels = km$cluster) Matt Wiener suggested: > col.vec <- c("black", "red", "blue") > text(chem.predict[,1:2], labels = km$cluster, color = col.vec[km$clusters]) Cheers Matt, that

Hello, I want to plot 21 scatter plots in a 7*3 matrix on a single A4 page (using mfcol). Below is an example, which -by and large- looks very close to the desired result. But I do not get it to fill out the whole page. If I specify the a different page size, the plots still do not fill out the whole page. Anybody an idea what to do? Maybe adjusting indvidiual plot size (if so, how?). Or is there

http://www.ncbi.nlm.nih.gov/pubmed/21418051 for the full reference. I don't have an electronic copy, but I do have that issue of Biometrics in my office. I'll have a copy sent over. Terry On 07/10/2012 04:08 PM, r-help-request at r-project.org wrote: > Send R-help mailing list submissions to > r-help at r-project.org > > To subscribe or unsubscribe via the World Wide

hi R-users: I want to draw three plot into one figure by layout and the script has been shown below. But I find R does not show the xlab and ylab information completely as shown the figure attached. How can I midify the script.? thank you . xxlab<-paste(cpmd," (",ro,"%)",sep=" ") yylab<-paste(rfmd," (",co,"%)",sep=" ")

I've created a hacked version of barplot() which is more compatible with the S version, but currently ONLY FOR VECTORS. Differences && new features: * The `space' argument is interpreted as the fraction of the average bar width. (The current version has the width of the bar plus the space in between constrained to sum to 1.) * There is a new argument `width'. * There is a

Hi, I'm trying to use the date library and I think I'm getting an incorrect warning when I make a plot with plot.date(). For example: > library(date) > x <- as.date(1:100) > y <- 1:100 > plot(x, y, log = "y") Warning message: parameter "log" couldn't be set in high-level plot() function The problem appears to be in the last line of plot.date()

Thomas, Thank you for your reply about the for () loop. The as.character advice worked. Sorry for the delay in getting back to?I had to set the project aside for a few weeks. This didn?t work exactly as is for (patient in as.character(1:n)){ pt <- MRN == patient (rest of the function) } But this did for (patient in as.character(levels(MRN))){ pt <- MRN == patient (rest of the function)

Hi all, To my knowledge, the current version of contour.default() does not handle the 'asp' parameter. This can be embarassing when displaying eg geographical maps, etc... Submitted to the opinion of more experienced R programmers, contour.defaut() function should be changed according to the followings: line 7: add = FALSE,asp=NA,...) line 33: plot.window(xlim, ylim,

Plotting vector fields as we know it from Wolfram Mathematica's PlotVectorField[{f1(x1), f2(x2)}, {x1, min, max}, {x2, ymin, ymax}, PlotPoints -> ...] ### FUNCTION DEFINITION #! Outer product with vector function ### # expand.outer joins functionality of "expand.grid" and "outer", so that # a vector function can be applied, values of which are then stored in

Hi all, I've written some code to generate 4 maps per screen and write the output to a jpeg. The output is fairly quick at the start (about 5 jpegs per minute) and then slows down greatly (1-2 jpegs per minute). Is there some way to speed it up? One of my thoughts is to keep the base map static on the screen and just update the points that are being plotted on the map (with the exception

[Env: R 2.12.2, Win XP] I'm creating figures using MASS::eqcsplot to provide equal scaling of the axes. My figures work OK when I plot to the screen, but when I try to do the same plot as a postscript file, I get an unexplicable error, > figframe() Error in if (yuin > xuin * ratio) yuin <- xuin * ratio else xuin <- yuin/ratio : missing value where TRUE/FALSE needed >

I read in a thread in r-help today that the windows device in 2.6 supports transparency, so I tried an example and had some issues. The density plots should be filled with transparent color in the following example (similar to the points), however the color is "fully" transparent. This works in the Cairo device, but not in the windows device. Thanks, --Matt Matt Austin

Hi,   I would like to distinguish my plots using colors but I got error message.  How do I correct that?     plot(ecdf(z), main ="CDF for observed and simulated weighted sum",type="l",lwd=2,col="blue", xlab="Weighted sum (mm)", ylab="Cumulative Percent", xlim=c(0,15), xaxs ='i', yaxs ='i',ylim=c(0,1)) par(new=TRUE)

Dear all, I have done further modifications on the histogram functions that I reported earlier this year, and I hope this can be of use and perhaps included in the distribution. I have been using this stuff a couple of months myself, and while it is nothing sophisticated, it has it's applications. :-) I did a few small modifications today to make it a bit more compact. I have modified the

I have the following codes but can anyone make it shorter i.e making these FOR loop into one loop ... thanks... par(mfrow=c(2,4)) for(i in 16:23){ hist(data[,i],main=paste(colnames(data)[i],sep=""),ylab="Frequency",xlim=c(1,5),xlab="Score",ylim=c(0,100)) } png("histogram.png") dev.off() par(mfrow=c(2,4)) for(i in

###Dear R users ###I have been using SensoMineR package from CRAN for most of my work in sensory data analysis and from my usage experience, I encountered some areas for improvement and considered ###modifying the function in SensoMineR package for my personal use. I felt that it could be useful to share this to the community for enabling adoption by other users where they might require a

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

I'm dealing with a classification problem using ordinal logistic regression. In the case of binary logistic regression with unequal proportions of 0's and 1's, a threshold in the interval [0,1] has to be adapted to transform back the predicted probabilities into 0 and 1. This can be done quite straightforward using e.g. the Kappa statistics as accuracy criterion. With