I have a set of 52 weekly values, one is the desired high lake stage for
a week, and the other is the desired low lake stage for each week. It
looks like this:
week
High
Low
1
16
14.5
2
16
14.5
3
15.95
14.45
4
15.84
14.34
5
15.73
14.23
6
15.61
14.11
7
15.5
14
8
15.38
13.88
9
15.25
13.75
10
15.13
13.63
....
52
15
13.5
I can plot out two nice curves using the following code, but I can't
figure out how to shade in the area lying between the two curves. I
also can't figure out how to adjust the line weight as lwd doesn't seem
to affect scatter.smooth. Thanks in advance.
l=11 #lower ylim
h=17 #upper ylim
par(new=F)
good <- complete.cases(week, high) #get rid of the NA error and only
using non-missing pairs
par(new=T)
scatter.smooth(week[good], high[good],family = "gaussian",span = .2,
ylim=c(l,h), type='n', lwd=2, color="blue", ylab="",
xlab="")
good <- complete.cases(week, low)
par(new=T)
scatter.smooth(week[good], low[good],family = "gaussian",span = .2,
ylim=c(l,h), type='n', lwd=2, color="red", ylab="Stage,
feet",
xlab="Week of Year") #loess, lower span makes more irregular
par(new=F)
Gregory A. Graves
Lead Scientist
Everglades REstoration COoordination and VERification (RECOVER)
Watershed Division
South Florida Water Management District
Phones: DESK: 561 / 682 - 2429
CELL: 561 / 719 - 8157
[[alternative HTML version deleted]]
?polygon
Bert Gunter
Genentech Nonclinical Biostatisics
-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
On
Behalf Of Graves, Gregory
Sent: Thursday, August 13, 2009 11:17 AM
To: r-help at r-project.org
Subject: [R] shading between two smoothed curves
I have a set of 52 weekly values, one is the desired high lake stage for
a week, and the other is the desired low lake stage for each week. It
looks like this:
week
High
Low
1
16
14.5
2
16
14.5
3
15.95
14.45
4
15.84
14.34
5
15.73
14.23
6
15.61
14.11
7
15.5
14
8
15.38
13.88
9
15.25
13.75
10
15.13
13.63
....
52
15
13.5
I can plot out two nice curves using the following code, but I can't
figure out how to shade in the area lying between the two curves. I
also can't figure out how to adjust the line weight as lwd doesn't seem
to affect scatter.smooth. Thanks in advance.
l=11 #lower ylim
h=17 #upper ylim
par(new=F)
good <- complete.cases(week, high) #get rid of the NA error and only
using non-missing pairs
par(new=T)
scatter.smooth(week[good], high[good],family = "gaussian",span = .2,
ylim=c(l,h), type='n', lwd=2, color="blue", ylab="",
xlab="")
good <- complete.cases(week, low)
par(new=T)
scatter.smooth(week[good], low[good],family = "gaussian",span = .2,
ylim=c(l,h), type='n', lwd=2, color="red", ylab="Stage,
feet",
xlab="Week of Year") #loess, lower span makes more irregular
par(new=F)
Gregory A. Graves
Lead Scientist
Everglades REstoration COoordination and VERification (RECOVER)
Watershed Division
South Florida Water Management District
Phones: DESK: 561 / 682 - 2429
CELL: 561 / 719 - 8157
[[alternative HTML version deleted]]
______________________________________________
R-help at r-project.org mailing list
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.
Hi Gregory,> I can plot out two nice curves using the following code, but I can'tfigure out how to shade in the area lying between the two curves. Use polygon.> ?polygon>>> "Graves, Gregory" <ggraves@sfwmd.gov> 14/08/2009 4:17 am >>>I have a set of 52 weekly values, one is the desired high lake stage for a week, and the other is the desired low lake stage for each week. It looks like this: week High Low 1 16 14.5 2 16 14.5 3 15.95 14.45 4 15.84 14.34 5 15.73 14.23 6 15.61 14.11 7 15.5 14 8 15.38 13.88 9 15.25 13.75 10 15.13 13.63 .... 52 15 13.5 I can plot out two nice curves using the following code, but I can't figure out how to shade in the area lying between the two curves. I also can't figure out how to adjust the line weight as lwd doesn't seem to affect scatter.smooth. Thanks in advance. l=11 #lower ylim h=17 #upper ylim par(new=F) good <- complete.cases(week, high) #get rid of the NA error and only using non-missing pairs par(new=T) scatter.smooth(week[good], high[good],family = "gaussian",span = .2, ylim=c(l,h), type='n', lwd=2, color="blue", ylab="", xlab="") good <- complete.cases(week, low) par(new=T) scatter.smooth(week[good], low[good],family = "gaussian",span = .2, ylim=c(l,h), type='n', lwd=2, color="red", ylab="Stage, feet", xlab="Week of Year") #loess, lower span makes more irregular par(new=F) Gregory A. Graves Lead Scientist Everglades REstoration COoordination and VERification (RECOVER) Watershed Division South Florida Water Management District Phones: DESK: 561 / 682 - 2429 CELL: 561 / 719 - 8157 [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list 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. ****************************************************************************** This email message, including any attached files, is confidential and intended solely for the use of the individual or entity to whom it is addressed. The NSW Department of Commerce prohibits the right to publish, copy, distribute or disclose any information contained in this email, or its attachments, by any party other than the intended recipient. If you have received this email in error please notify the sender and delete it from your system. No employee or agent is authorised to conclude any binding agreement on behalf of the NSW Department of Commerce by email. The views or opinions presented in this email are solely those of the author and do not necessarily represent those of the Department, except where the sender expressly, and with authority, states them to be the views of NSW Department of Commerce. The NSW Department of Commerce accepts no liability for any loss or damage arising from the use of this email and recommends that the recipient check this email and any attached files for the presence of viruses. ****************************************************************************** [[alternative HTML version deleted]]
Graves, Gregory wrote:> I have a set of 52 weekly values, one is the desired high lake stage for > a week, and the other is the desired low lake stage for each week. It > looks like this: > ... > I can plot out two nice curves using the following code, but I can't > figure out how to shade in the area lying between the two curves. I > also can't figure out how to adjust the line weight as lwd doesn't seem > to affect scatter.smooth. Thanks in advance. >Hi Gregory, You can use "polygon" by running across either the upper or lower values, then backwards across the other set. For an example of how to do this, look at the source code of the dispersion function in the plotrix package. Jim