[Thread Prev][Thread Next][Index]
shaded area between two lines
Hi,
I would like to shade the area between two lines.
Is there a straightforward way to do this?
I started using the 3 argument shade command, but I do not
quite get what I want.
I want to represent the shaded area +/- one standard
deviation around something along my time axis. So I calculated
the 2d-variables:
TT=(t-dt, t+dt) XX=(x-std, x+std) and made up a 2d-dummy FF
of the same size containing only 1's.
Then
shade FF,TT,XX
does look right, except that I do not have the time axis I had
before any more. Somehow, I can't get the right time axis
into that plot.
This would be nice, because one could them overlay some other
data one wants to compare to the indicated range.
However, I fail getting any overlay to work...
Here is what I do:
I have 15 min data:
define axis/t="15-AUG-1998 00:07:30":"31-DEC-1999
23:52:30":0.010416667/unit="days since 01-Jan-1998 00:00:00" timtax
define axis/x=1:2:1 timxax2
define axis/x=1:4:1 timxax4
define grid/t=timtax/x=timxax2 twgrd2
define grid/t=timtax/x=timxax4 twgrd4
! time series
set data/ez/var=w/grid=twgrd4/col=4 wrccrf_swc_199899_tim_d.dat
! data file with the data for the 3 argument shade
set data/ez/var=ff,xx,tt/grid=twgrd2/col=3 tt
shade ff,tt,xx
plot/over w[d=1,i=@ave]
This plot over now does not understand the axes from the previous shade.
Vice versa,
plot w[d=1,i=@ave]
then
shade/over/noaxis ff,tt,xx
draws a new plot.
Any ideas would be most appreciated.
Cheers,
Joerg
--
Jörg Kaduk Tel.: 1 650 325 1521 x 416
Carnegie Institution of Washington FAX: 1 650 325 6857
Dept. of Plant Biology
260 Panama Street joerg@jasper.stanford.edu
Stanford, CA 94305-1297 http://Jasper.Stanford.EDU/joerg/
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement