[R] Plotting a broken line?

Greg Snow Greg.Snow at intermountainmail.org
Wed Mar 7 18:30:38 CET 2007


If you insert an NA (or row of NA's) into the data at each place you
want a break (after indicator increases), then the regular plot with
type='l' will break the line for you.

Is this what you want?

-- 
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.snow at intermountainmail.org
(801) 408-8111
 
 

> -----Original Message-----
> From: r-help-bounces at stat.math.ethz.ch 
> [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Aldi Kraja
> Sent: Wednesday, March 07, 2007 9:21 AM
> To: r-help at stat.math.ethz.ch
> Subject: [R] Plotting a broken line?
> 
> Hi,
> 
> Is there a smart way in the R graphs to create a line that is 
> broken in intervals based on the indicator given below.
> following is a small test graph
> 
> Location,indicator,otherinfo
> 1.2,1,2.2
> 2.5,1,2.5
> 3.7,1,2.3
> 20.1,2,4.3
> 
> 22.5,2,5.2
> 25.0,2,3.4
> 27.3,2,2.2
> 
> 35.1,3,3.4
> 37.0,3,7.2
> 38.0,3,6.1
> 40.1,3,5.4
> 52.9,3,3.3
> 
> Right now in the plot the line is continuous, but I would 
> like to have it broken based on the indicator. If the line of 
> the plot reaches the last observation of indicator=1 then the 
> line needs to stop; the next line will start at location 22.5 
> and continue up top 27.3; the next line goes from 35.1 up to 52.9.
> 
>  > x<-read.table(file='c:\\aldi\\testgraph.csv',sep=',',header=T)
>  > x
>    Location indicator otherinfo
> 1       1.2         1       2.2
> 2       2.5         1       2.5
> 3       3.7         1       2.3
> 4      20.1         2       4.3
> 5      22.5         2       5.2
> 6      25.0         2       3.4
> 7      27.3         2       2.2
> 8      35.1         3       3.4
> 9      37.0         3       7.2
> 10     38.0         3       6.1
> 11     40.1         3       5.4
> 12     52.9         3       3.3
> 
>  >
> plot(x$Location,x$indicator,type='l',xlim=c(0,max(x$Location))
> ,ylim=c(0,max(x$indicator,x$otherinfo)))
>  > points(x$Location,x$otherinfo)
> 
> TIA,
> Aldi
> 
> --
> 
> ______________________________________________
> R-help at stat.math.ethz.ch 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.
>



More information about the R-help mailing list