[R] multiple scale

[Jeanie (Jie) Na]

Hi, 
Search the list, you might be able to find the message but here it is.

Copied From Dr. Brian Ripley's post
> x <- 1:10
> y <- rnorm(10)
> z <- runif(10, 1000, 10000)
> par(mar=c(5,4,4,4) + 0.1)  # Leave space for z axis
> plot(x, y)
> par(new=T)
> plot(x, z, type="l", axes=F, bty="n", xlab="", ylab="")
> axis(4, at=pretty(range(z)))
> mtext(z, 4, 3)
change last line to mtext("z", 4, 3).

J.
On Tue, 2005-08-02 at 17:18 +0100, herodote at oreka.com wrote:
> Hi all
> 
> i need to put on one graph 2 functions who's x axis is the same and y not.
> 
> I mean on horizontal the time, and on vertical left: pressure, on vertical right: rpm of a motor, is R able to do that?
> 
> i've found this that i could adapt maybe (i don't need time series really?) :/ :
> 
> (http://tolstoy.newcastle.edu.au/R/help/04/03/1456.html)
> ## 
> ## Description: A simple function which plots two time series on one plot where 
> ## the series can have different value intervals over the same time interval. 
> ## Usage: ts.plot.2Axis(xleft, xright) 
> ## Arguments: xleft is the time series for the left vertical axis and xright 
> ## is for the right axis. xleft and xright are defined as time series with 
> ## the 'ts' function in package ts. 
> ## ts.plot function must be available, do library(ts) to ensure this if 
> ## necessary. 
> ## In addition the usual 'ts.plot' and 'plot' parameters can be set 
> ## directly (mar, main, xlab, ylab, lwd) or through gpars as in ts.plot. 
> ## Also parameter digits is the preferred number of decimal digits on right 
> ## axis and ticks is the preferred number of tick marks on right axis. 
> ## Details: The time series for the right vertical axis is scaled with a simple 
> ## rule of thumb scaling. 
> ## The ts.plot function is used to plot the series. 
> ## Value: None. 
> ## Note: When scaling is not acceptable try switching the series parameters. 
> ## If a ylabel is to be set it is here only possible for the left axis. 
> ## See also: 'ts.plot', 'ts', 'legend'. 
> ## Author and date: Hauksson, Bjorn Arnar. March 2004. 
> ## Example: 
> ## First paste this function into the R console or use 'source'. 
> #library(ts) 
> #data(UKLungDeaths) 
> #x <- ldeaths 
> #y <- fdeaths/mdeaths 
> #ts.plot.2Axis(x, y) 
> #legTxt <- c("UK lung deaths", "UK female/male deaths (rhs)") 
> #legend(1976.5, 3950, legTxt, lty=c(1:2), col=c(1:2), lwd=2, bty="n") 
> ## 
> 
> 
> ts.plot.2Axis <- function(xleft, xright, digits=1, ticks=5, 
>                           mar=(c(4,4,4,4)+0.1), main="", 
>                           xlab="", ylab="", lwd=2, gpars=list()) { 
>         # Settings for other parameters than those in the function parameter list 
>         par(mar=mar) # Margins 
>       k <- ncol(as.matrix(xleft)) # Number of time series on left vertical scale 
>         lty <- c(1:(k+1)) # Line types 
>         col <- c(1:(k+1)) # Line colors 
> 
> 
>         # Scale time series on right vertical axis 
>         scale <- (max(xleft)-min(xleft))/(max(xright)-min(xright)) 
>       xright2 <- xright*scale 
>       meanScale <- mean(xleft) - mean(xright2) 
>         xright2 <- xright2 + meanScale 
> 
> 
>         # Plot the series 
>         ts.plot(xleft, xright2, lty=lty, col=col, main=main, ylab=ylab, xlab=xlab, 
>         lwd=lwd, gpars=gpars) 
> 
> 
>         # Add the right vertical axis labels 
>         lab <- seq(round(min(xright), digits), round(max(xright), digits), 
>         length=ticks) 
>         labAt <- seq(min(xright2), max(xright2), length=ticks) 
>         axis(side=4, labels=lab, at=labAt) 
> } 
> 
> are there a better way to plot what i want?
> 
> thks.
> guillaume.
> 
> ////////////////////////////////////////////////////////////
> // Webmail Oreka : http://www.oreka.com
> ////////////////////////////////////////////////////////////
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
> 
-- 
Jeanie (Jie) Na

Programmer Analyst II                       _   _    ___    _   _    ____
Department of Quantitative Health Sciences [_]-[_]  / - \  | |_| |  (_(_
Cleveland Clinic Foundation                 |   |  ( |_| ) )  _  (   _| )
Tel: (216)4451369                          [_]-[_]  \_\_\  |_| |_|  (___/