similar to: Title for a group of plots?

I would like the text plotted with 'mtext' to be alighned like it is for printing on the console. Here is what I have: > print(emt) ME RMSE MAE MPE MAPE MASE original -1.034568e+07 1.097695e+08 2.433160e+07 -31.30554 37.47713 1.5100050 xreg 1.561235e+01 2.008599e+03 9.089473e+02 267.05490 280.66734

I want to place four plots on a page, and I would like to have all four plots share a common title. I have tried the following code, but the title is centered over the fourth graph and not centered across all four plots. Does anyone have any suggestions? R 2.1.1 windows xp oldpar<-par(mfcol =c(1,4),ask=TRUE) plot(p,varp) plot(p,SEp) plot(p,CVp) plot(p,ppval) title(paste("P and 95%CI

R 2.15.1 OS X (MLion) Colleagues, I am aware that changes in mfrow / mfcol in par() affect cex (from help: In a layout with exactly two rows and columns the base value of ?"cex"? is reduced by a factor of 0.83: if there are three or more of either rows or columns, the reduction factor is 0.66). I generate a multipage PDF in which mfrow varies such that cex is impacted. This affect

I am trying to upgrade to R 2.14 from R 2.13.1 I have compied all the libraries from the 'library' directory in my existing installation (2.13.1) to the installed R 2.14. Now I want to uninstall the old installation (R 2.13.1) and I get the error: Internal Error: Cannot find utCompiledCode record for this version of the uninstaller. Any ideas? Kevin [[alternative HTML

I noticed that when I fit a linear model using 'lm' there is an attribute called "factors" that is added to the "term". It doesn't seem to appear for 'model.matrix', just 'lm'. I have been unable to find where it gets constructed or what it means? It looks like a two dimensional array that I may be able to use so I would just like to get some

I can get an array of strings for the data that I want using 'paste()' as follows: paste('ma', 1:am$arma[2], '=', coef(am)[1:am$arma[2] + am$arma[1]], sep='') This results in a vector of strings like: [1] "ma1=1.17760133668255" "ma2=0.649795570407939" "ma3=0.329456750858276" What I would like is fixed.pars <-

Can someone enlighten me on why the following doesn't work? setwd('C:/Temp/R') d <- rep(1:53,2) (s <- ts(d, frequency=53, start=c(2000,10))) n <- length(s) k <- n%/%3 for(i in (n-k):n) { st <- c(start(s)[1] + (start(s)[2] + i)%/%frequency(s), (start(s)[2] + i) %% frequency(s)) ed <- c(start(s)[1] +

I noticed that the text() command adds text to a plot. Is there a way to either make the plot blank or add text to a "blank sheet". I would like to "plot" a page that contains just text, no plot lines, labels, etc. Suggestions? Kevin [[alternative HTML version deleted]]

Dear all, I'm a brazilian medical doctor that I research use of likelihood in the medical research. I need to do a plot with this layout: --------------------------------------- | Major Title | | | | Title Title | | -------- -------- | | | | | | | | | |

In browsing the source I see the following construct: res <- switch(type, working = , response = r, deviance = , pearson = if (is.null(object$weights)) r else r * sqrt(object$weights), partial = r) I understand that 'switch' will execute the code that is matched by its corresponding string value (in this case 'type'). What I don't

I have two time series a <- ts(1:10, start=c(1,6), end=c(2,5), frequency=10) b <- ts(1:5, start=c(2,1), end=c(2,5), frequency=10) Obviously 'b' is a subset of 'a'. I want a single index value indicating where that start of 'b' lines up with the start of 'a'. So in this simple example I would expect an index of 5. I was playing with 'merge'.

Hello! I want to display four xyplots in the same window (sorry if I don't use the propers words, but my English is rude) I know that if I use plot(y~x) then x=(1:100) y=(1:100) par(mfrow=c(2,2)) plot(y~x) plot(y~x) plot(y~x) plot(y~x) would give what I want, but this is not the case because I need the functions of xyplot. So I did the following: plot.new()

Dear all, I have a 3x2 plot and in addition to the title of the individual plots I would like to have an overall title for each row. I managed to get an overall title for the whole plot matrix with mtext: par(mfrow=(c(3,2)), mar=c(6.4,4.5,4.2, 1.8), oma=c(0,0,3,0)) for (i in 1:6) barplot(sample(1:10,5), main=paste("Plot #",i)) mtext("Overall Title Row 1", outer=TRUE) but I

What is wrong with the following? x <- 1:2 if(x[1] > 0) { if(x[2] > 0) { print("1 & 2 > 0") } else { print("1 > 0") } } else { if(x[2] > 0) {

Dear helpers, In an attempt to use a loop to generate graphs in a for loop in run into a problem. The plan is to fill each page with eight graphs (mfrow = c(4,2)) in to two columns. Only the buttom graphs ( meaning every fourth graph) have tick labels on the x axis to preserve space. I used an if .... Else statement to achieve that. The problem is that the first eight graphs are skipped

I have some data where the frequency is heavily weighted on the lower end. So I have lots of low values with very few higher values. I would like to find breakpoints that cover the data with as much detail as possible. I find that if I use hist() to automatically find the breaks for me it finds breaks that are too coarse for the low values. I have tried the other algorithms (like 'Scott'

I can put a legend on a plot with something like: legend('bottom', leg.txt, horiz = TRUE, fill = colors) But what if the arrays leg.txt and colors are too big? I would still like to provide a legend but to save space I would like to just show small boxes with the color filled in so it will still fit on the plot. If I could adjust the size of the boxes (at least in one dimension)

I'm trying to use mtext to create a main title over multiple plots. Below is a simple self-contained example and my sessionInfo (I should note I've also tried this with R-2.8.1 with the same results). When I execute the code chunk below, I get the plots, but no title. I've tried this using the screen driver, pdf, and postscript. I've used different sizes of paper. I suspect I am

Hi all, I'm trying to get the min and max of a sequence of number using a loop like the folllowing. Can anyone point me to why it doesn't work. Thanks. type <- c("min","max") n <- 1:10 for (a in 1:2) { print(type[a](n)) } -- Muhammad

I started to check what I thought I knew with autocovariance and it doesn’t jive with the the calculations given by ‘R’. I was wondering if there is some scaling or something that I am not aware of. Take the example Ø d <- 1:10 Ø (a <- acf(d, type="covariance", demean=FALSE, plot=FALSE)) Autocovariances of series ‘d’, by lag 0 1 2 3 4 5 6