Help with the plot function

Dear all,

I found the last example of this link (https://sites.ualberta.ca/~
lkgray/uploads/7/3/6/2/7362679/6c_-_line_plots_with_error_bars.pdf) very
similar to the one I need to make for my paper, and I think I got what I
wanted by applying some of the suggestions of this mail list.

Here it is the code I devised (maybe there will be further improvements
from the list):


YEAR <- c(1996 , 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004 )
T_MAR <- c(2.8, 6.5, 5.4,2.4, 4, 4.1, 3, 4.4, 4.5)
T_APR <- c(5.7, 7.8, 7.7, 4.6, 4.7, 6.2,5.7, 5.9, 7)
T_MAY <- c(7, 8.8, 10, 6, 5.5, 7.6, 8.5, 7.3, 10.2)
BUD <- c(87, 98, 93, 85, 89, 91, 87, 92, 92)
BUD_SE <- c(3.6, 2, 2.4, 4, 2.4, 2.4, 4, 2.4, 3)
g1 <- data.frame(YEAR, T_MAR, T_APR, T_MAY, BUD, BUD_SE)

par(mar=c(5,5,5,5))
plot(YEAR,T_MAR,pch=0,type="b",col="green",yaxt="n",ylim=c(-
15,12),ylab="")
lines(T_MAR~YEAR, type="o", pch=19, col="green")
lines(T_APR~YEAR, type="o", pch=19, col="red")
lines(T_MAY~YEAR, type="o", pch=19, col="blue")
axis(side=2, at=c(0,6,12))
mtext("Temperature (?C)", side = 2, line=2.5, at=6)
legend(1999,14, bty="n",
       lty = c(1,1,1),
       lwd = c(2,2,2),
       pch = c(19,19,19),
       col = c("green","red","blue"), legend = c("March","April","May"))

par(new=TRUE)

plot(YEAR,BUD, pch=1,type="b",col="blue",yaxt="n",ylim=c(105,50), ylab="")
arrows(g1$YEAR,g1$BUD, g1$YEAR,g1$BUD + g1$BUD_SE, length=0.05, angle=90)
arrows(g1$YEAR,g1$BUD, g1$YEAR,g1$BUD-g1$BUD_SE, length=0.05, angle=90)
axis(side=4, at=c(100,90,80), labels=c("80","90","100"),las=1)
mtext("Bud Break (Julian Day)", side = 4, line=2.5, at=90, padj=0)
abline(h=78)


Thanks for the insight and help as to the functions I needed to take a
look at.

Andre


On Tue, Jun 20, 2017 at 12:38 AM, PIKAL Petr <petr.pikal at precheza.cz>
wrote:

Hi

You are quite close. With slight modification of your code:

par(mfrow = c(2, 1))
par(cex = 0.6)
par(mar = c(0, 0, 0, 0), oma = c(4, 4, 0.5, 0.5))
par(tcl = -0.25)
par(mgp = c(2, 0.6, 0))

plot(BUD~YEAR, type="o", ann=F, axes=F, pch=19, ylim=c(60,100),data=g1)
axis(4, las=2)
mtext("Bud Break (Julian Day)", side=4, padj=4)
arrows(g1$YEAR,g1$BUD, g1$YEAR,g1$BUD + g1$BUD_SE, length=0.05, angle=90)
arrows(g1$YEAR,g1$BUD, g1$YEAR,g1$BUD-g1$BUD_SE, length=0.05, angle=90)

plot(T_MAR~YEAR, type="l", pch=19, ann=F, axes=F, xlim=c(1996,2004),
ylim=c(0,12), data=g1)

title(ylab="Temperature (?C)",xlab="Year")
axis(1, at=seq(1996, 2004, 2))
axis(2, at=c(0,3,6,9,12), las=2)

I am quite close to what you probably expect. You need modify axes and
their annotation, which is left for your training.

Cheers
Petr

-----Original Message-----
From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Andr?

Luis

Neves
Sent: Tuesday, June 20, 2017 6:29 AM
To: Bert Gunter <bgunter.4567 at gmail.com>
Cc: R mailing list <r-help at r-project.org>
Subject: Re: [R] Help with the plot function

Hi, Bert:
Yes, I studied the functions you suggested, but I didn't get to adapt

it to

my example whose reproducible code I sent in my first email.

Here it is the code of the functions I studies:

## par
par(mfrow = c(2, 3))
par(cex = 0.6)
par(mar = c(3, 3, 0, 0), oma = c(1, 1, 1, 1))
for (i in 1:6) {
 plot(1, 1, type = "n")
mtext(letters[i], side = 3, line = -1, adj = 0.1, cex = 0.6)}


par(mfrow = c(2, 3))
par(cex = 0.6)
par(mar = c(0, 0, 0, 0), oma = c(4, 4, 0.5, 0.5))
par(tcl = -0.25)
par(mgp = c(2, 0.6, 0))

for (i in 1:6) {
      plot(1, axes = FALSE, type = "n")
      mtext(letters[i], side = 3, line=-1,adj=0.1,cex=0.6, col =

"grey40")

  if (i %in% c(4, 5, 6))
    axis(1, col = "grey40", col.axis = "grey20", at = seq(0.6,
                                                          1.2, 0.2))
  if (i %in% c(1, 4))
    axis(2, col = "grey40", col.axis = "grey20", at = seq(0.6,
                                                          1.2, 0.2))
  box(col = "grey60")}

## Layout

m <- rbind(c(1, 1), c(2, 3))
m
layout(m)
layout.show(3)
layout(m)
par(mar = c(3, 3, 0, 0))
for (i in 1:3) plot(1, 1, type = "n")

Thank you,

Andre


On Mon, Jun 19, 2017 at 9:53 PM, Bert Gunter <bgunter.4567 at gmail.com>
wrote:

1. Did you study the functions (esp. ?layout) to which I referred you?

2. Show us your code! -- "to no avail" is meaningless!

-- Bert


Bert Gunter

"The trouble with having an open mind is that people keep coming along
and sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )


On Mon, Jun 19, 2017 at 7:01 PM, Andr? Luis Neves <

andrluis at ualberta.ca>

wrote:

I'm trying to recreate a graph similar to the last one found on this

link:

https://sites.ualberta.ca/~lkgray/uploads/7/3/6/2/

7362679/6c_-_line_plots_with_error_bars.pdf

The difference is that I want budbreak on the top, and the

temperatures

at

the bottom.

I tried to set par before each graph and include lines, with no

avail.

Thanks, Bert.

Andre

On Mon, Jun 19, 2017 at 7:41 PM, Bert Gunter <

bgunter.4567 at gmail.com>

wrote:

See

?layout
?split.screen
?par  (the mfrow and mfcol values)

depending exactly on what you want to do and how you want to do it.
Essentially, these all allow you to make separate plots at

different

regions of the device.


Cheers,
Bert




Bert Gunter

"The trouble with having an open mind is that people keep coming

along

and sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )


On Mon, Jun 19, 2017 at 6:30 PM, Andr? Luis Neves

<andrluis at ualberta.ca

wrote:

Dear friends,

I have the following dataframe:

YEAR <- c(1996 , 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004 )
T_MAR <- c(2.8, 6.5, 5.4,2.4, 4, 4.1, 3, 4.4, 4.5)
T_APR <- c(5.7, 7.8, 7.7, 4.6, 4.7, 6.2,5.7, 5.9, 7)
T_MAY <- c(7, 8.8, 10, 6, 5.5, 7.6, 8.5, 7.3, 10.2)
BUD <- c(87, 98, 93, 85, 89, 91, 87, 92, 92)
BUD_SE <- c(3.6, 2, 2.4, 4, 2.4, 2.4, 4, 2.4, 3)
g1 <- data.frame(YEAR, T_MAR, T_APR, T_MAY, BUD, BUD_SE)

###PLOT
dev.new(width=6.5, height=5)
par (cex=1, family="sans", mar=c(5,5,5,5.5))
plot(T_MAR~YEAR, type="l", pch=19, ann=F, axes=F,

xlim=c(1996,2004),

ylim=c(0,12), data=g1)

title(ylab="Temperature (?C)",xlab="Year")
axis(1, at=seq(1996, 2004, 2))
axis(2, at=c(0,3,6,9,12), las=2)
par(new=T)
plot(BUD~YEAR, type="o", ann=F, axes=F, pch=19,

ylim=c(60,100),data=g1)

axis(4, las=2)
mtext("Bud Break (Julian Day)", side=4, padj=4)
arrows(g1$YEAR,g1$BUD, g1$YEAR,g1$BUD + g1$BUD_SE, length=0.05,
angle=90)
arrows(g1$YEAR,g1$BUD, g1$YEAR,g1$BUD-g1$BUD_SE, length=0.05,

angle=90)

However, I'd like to draw a multi-panel graph with budbreak on

the top

(as
it is), and with the temperatures for March, April, and May on

the

bottom,
with their respective legends.

I was wondering if you could help me out with this.

Thanks a million for your help.

--
Andre

        [[alternative HTML version deleted]]

see

code.

--
Andre

--
Andre

      [[alternative HTML version deleted]]

ng-guide.html

--
Andre