Hi Jake, My mistake, that should be: par(xpd=FALSE) abline(lm(...)) par(xpd=NA) I suspect that the odd looking regression lines are due to points off the plot as you have specified xlim and ylim for your plots. Jim On Wed, Jan 18, 2017 at 2:22 PM, Jake William Andrae
<jake.andrae at adelaide.edu.au> wrote:
If I run the script like this, with the commands you suggested, it looks like the figure attached. Scratching my head, the regression lines shouldn?t look like this.
#Growing season
#Construction of plots (without Rhmax)
# 12 figures arranged in 3 rows and 3 columns
attach(mtcars)
par(mfrow=c(3,3), mar = c(.6,.5,.6,1), oma = c(5,4,2,2), xpd = NA)
#Concentration
plot(Growing_season_precipitation, Concentration, xaxs = "i", yaxs = "i", xlim = c(0, 700), ylim = c(0,500), xlab = NA, ylab = "Concentration", xaxt='n', pch=21, bg='black', rect(par("usr")[1],par("usr")[3],par("usr")[2],par("usr")[4],col = "white"))
par(xpd=TRUE)
abline(lm(Growing_season_precipitation~Concentration))
par(xpd=NA)
plot(Growing_season_VPD, Concentration, xaxs = "i", yaxs = "i", xlim = c(0.6, 1.8), ylim = c(0,500), xlab = NA, ylab = "", xaxt='n', pch=21, bg='black', yaxt='n', rect(par("usr")[1],par("usr")[3],par("usr")[2],par("usr")[4],col = "white"))
par(xpd=TRUE)
abline(lm(Growing_season_VPD~Concentration))
par(xpd=NA)
plot(Growing_season_temperature, Concentration, xaxs = "i", yaxs = "i", xlim = c(20,34), ylim = c(0,500), xlab = NA, ylab = "", xaxt='n', pch=21, bg='black', yaxt='n', rect(par("usr")[1],par("usr")[3],par("usr")[2],par("usr")[4],col = "white"))
par(xpd=TRUE)
abline(lm(Growing_season_temperature~Concentration))
par(xpd=NA)
#ACL
plot(Growing_season_precipitation, ACLTotal, xaxs = "i", yaxs = "i", xlim = c(0, 700), ylim = c(28,32), xlab = "", ylab = "ACL", xaxt='n', pch=21, bg='black', rect(par("usr")[1],par("usr")[3],par("usr")[2],par("usr")[4],col = "white"))
par(xpd=TRUE)
abline(lm(Growing_season_precipitation~ACLTotal))
par(xpd=NA)
plot(Growing_season_VPD, ACLTotal, xaxs = "i", yaxs = "i", xlim = c(0.6, 1.8), ylim = c(28,32), xlab = "", ylab = "", xaxt='n', pch=21, bg='black', yaxt='n' ,rect(par("usr")[1],par("usr")[3],par("usr")[2],par("usr")[4],col = "white"))
par(xpd=TRUE)
abline(lm(Growing_season_VPD~ACLTotal))
par(xpd=NA)
plot(Growing_season_temperature, ACLTotal, xaxs = "i", yaxs = "i", xlim = c(20,34), ylim = c(28,32), xlab = "", ylab = "", xaxt='n', pch=21, bg='black', yaxt='n', rect(par("usr")[1],par("usr")[3],par("usr")[2],par("usr")[4],col = "white"))
par(xpd=TRUE)
abline(lm(Growing_season_temperature~ACLTotal))
par(xpd=NA)
#CPI
plot(Growing_season_precipitation, CPITotal, xaxs = "i", yaxs = "i", xlim = c(0, 700), ylim = c(0,30), xlab = "Total precipitation (mm)", ylab = "CPI", pch=21, bg='black', rect(par("usr")[1],par("usr")[3],par("usr")[2],par("usr")[4],col = "white"))
par(xpd=TRUE)
abline(lm(Growing_season_precipitation~CPITotal))
par(xpd=NA)
plot(Growing_season_VPD, CPITotal, xaxs = "i", yaxs = "i", xlim = c(0.6, 1.8), ylim = c(0,30), xlab = "Average daily VPD (kPa)", ylab = "", pch=21, bg='black', yaxt='n' ,rect(par("usr")[1],par("usr")[3],par("usr")[2],par("usr")[4],col = "white"))
par(xpd=TRUE)
abline(lm(Growing_season_VPD~CPITotal))
par(xpd=NA)
plot(Growing_season_temperature, CPITotal, xaxs = "i", yaxs = "i", xlim = c(20,34), ylim = c(0,30), xlab = "Average daily temperature (oC)", ylab = "", pch=21, bg='black', yaxt='n', rect(par("usr")[1],par("usr")[3],par("usr")[2],par("usr")[4],col = "white"))
par(xpd=TRUE)
abline(lm(Growing_season_temperature~CPITotal))
par(xpd=NA)
#Plot main title
title(main="Growing season (June-November, inclusive)",outer=T)
-----Original Message-----
From: Jim Lemon [mailto:drjimlemon at gmail.com]
Sent: Wednesday, 18 January 2017 1:40 PM
To: Jake William Andrae <jake.andrae at adelaide.edu.au>; r-help mailing list <r-help at r-project.org>
Subject: Re: [R] Adding regression line to each individual plot in a window with multiple plots
Hi Jake,
In the second line of your script, you set xpd=NA. That means that abline will draw a line across the entire display region rather than restricting it to the plot region. If you must set this at the beginning, then add these lines:
par(xpd=TRUE)
abline(lm(...))
par(xpd=NA)
for each regression line you display.
Jim
On Wed, Jan 18, 2017 at 1:41 PM, Jake William Andrae <jake.andrae at adelaide.edu.au> wrote:
Hi Jim, that works but this is what it then looks like. -----Original Message----- From: Jim Lemon [mailto:drjimlemon at gmail.com] Sent: Wednesday, 18 January 2017 1:06 PM To: Jake William Andrae <jake.andrae at adelaide.edu.au> Cc: r-help at r-project.org Subject: Re: [R] Adding regression line to each individual plot in a window with multiple plots Hi Jake, As I don't have your data set, try this: attach(mtcars) plot(mpg~disp,xaxs="i",yaxs="i") abline(lm(mpg~disp)) Jim On Wed, Jan 18, 2017 at 12:04 PM, Jake William Andrae <jake.andrae at adelaide.edu.au> wrote:
Hi Everyone,
I've constructed a script that adds multiple plots to the plot window, but I'm having a bit of trouble adding a regression line to each individual plot. Of course, the regression lines will vary depending on the variables plotted against one another. I've attached the script.
#Growing season
#Construction of plots
# 12 figures arranged in 3 rows and 4 columns
attach(mtcars)
par(mfrow=c(3,4), mar = c(.6,.5,.6,1), oma = c(15,4,2,2), xpd = NA)
#Concentration plot(Growing_season_precipitation, Concentration, xaxs
= "i", yaxs = "i", xlim = c(0, 700), ylim = c(0,500), xlab = NA, ylab
= "Concentration", xaxt='n', pch=21, bg='black',
abline(lm(Growing_season_precipitation~Concentration)))
plot(Growing_season_VPD, Concentration, xaxs = "i", yaxs = "i", xlim
= c(0.6, 1.8), ylim = c(0,500), xlab = NA, ylab = "", xaxt='n',
pch=21, bg='black', yaxt='n') plot(Growing_season_Rhmax,
Concentration, xaxs = "i", yaxs = "i", xlim = c(35, 60), ylim =
c(0,500), xlab = NA, ylab = "", xaxt='n', pch=21, bg='black',
yaxt='n') plot(Growing_season_temperature, Concentration, xaxs = "i",
yaxs = "i", xlim = c(20,34), ylim = c(0,500), xlab = NA, ylab = "",
xaxt='n', pch=21, bg='black', yaxt='n') #ACL
plot(Growing_season_precipitation,
ACLTotal, xaxs = "i", yaxs = "i", xlim = c(0, 700), ylim = c(28,32),
xlab = "", ylab = "ACL", xaxt='n', pch=21, bg='black')
plot(Growing_season_VPD, ACLTotal, xaxs = "i", yaxs = "i", xlim =
c(0.6, 1.8), ylim = c(28,32), xlab = "", ylab = "", xaxt='n',
pch=21, bg='black', yaxt='n') plot(Growing_season_Rhmax, ACLTotal,
xaxs = "i", yaxs = "i", xlim = c(35, 60), ylim = c(28,32), xlab = "",
ylab = "", xaxt='n', pch=21, bg='black', yaxt='n')
plot(Growing_season_temperature, ACLTotal, xaxs = "i", yaxs = "i",
xlim = c(20,34), ylim = c(28,32), xlab = "", ylab = "", xaxt='n',
pch=21, bg='black', yaxt='n') #CPI
plot(Growing_season_precipitation,
CPITotal, xaxs = "i", yaxs = "i", xlim = c(0, 700), ylim = c(0,30),
xlab = "Total precipitation (mm)", ylab = "CPI", pch=21,
bg='black') plot(Growing_season_VPD, CPITotal, xaxs = "i", yaxs =
"i", xlim = c(0.6, 1.8), ylim = c(0,30), xlab = "Average daily VPD
(kpa)", ylab = "", pch=21, bg='black', yaxt='n')
plot(Growing_season_Rhmax, CPITotal, xaxs = "i", yaxs = "i", xlim =
c(35, 60), ylim = c(0,30), xlab = "Average daily RHmax (%)", ylab =
"", pch=21, bg='black',
yaxt='n') plot(Growing_season_temperature, CPITotal, xaxs = "i", yaxs
= "i", xlim = c(20,34), ylim = c(0,30), xlab = "Average daily
temperature (oC)", ylab = "", pch=21, bg='black', yaxt='n') #Plot
main title title(main="Growing season (June-November,
inclusive)",outer=T)
Any help would be greatly appreciated!
[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.