[R] shading between two smoothed curves

[Bert Gunter]

?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.