Re: [R] "Groups" in XYPLOT
Much nicer, thanks, Thomas. I've made a small change to make the
differences more obvious for this example.
#First dataset
N <- 25
x <- round(rnorm(N),2)
y <- round(rnorm(N),2)
df.1 <- data.frame(x = x, y = y)
#Plot it
xyplot(x~y,xlab="Test 1 Data",ylab="P(A>A*)", data=df.1)
#Second Dataset
N <- 20
x <- round(rnorm(N),2) + 10
y <- round(rnorm(N),2)
df.2 <- data.frame(x = x, y = y)
#How to get this in the same panel as plot 1?
xyplot(x~y,xlab="Test 2 Data",ylab="P(A>A*)", data=df.2)
df <- make.groups(df.1, df.2) ## Thanks, Deepayan!
xyplot(x~y,
groups=which,
xlab="All Test Data",
ylab="P(A>A*)",
auto.key=list(space="right"),
data=df)
Cheers
Andrew
On Sat, Mar 17, 2007 at 11:48:05PM -0500, Thomas Colson wrote:
> Sorry 'bout that.
>
> Here's about as simple as I can get it:
>
> #First dataset
> N <- 25
> x <- round(rnorm(N),2)
> y <- round(rnorm(N),2)
> df <- data.frame(x = x, y = y)
> #Plot it
> xyplot(df$x~df$y,xlab="Test 1 Data",ylab="P(A>A*)")
>
> #Second Dataset
> N <- 20
> x <- round(rnorm(N),2)
> y <- round(rnorm(N),2)
> df <- data.frame(x = x, y = y)
> #How to get this in the same panel as plot 1?
> xyplot(df$x~df$y,xlab="Test 2 Data",ylab="P(A>A*)")
>
>
>
>
> As far as combining that data into one datasetI'm trying to plot both
> datasets into one panel to show the differences in the plots of the
> two...utilizing the coloring/symbology funtions...wouldn't combining them
> null that capability?
>
> Thanks for your replies!
>
>
> Thomas Colson, PhD
> North Carolina State University
> Department of Forestry and Environmental Resources
> (919)673-8023
> [EMAIL PROTECTED]
>
> Schedule: www4.ncsu.edu/~tpcolson
> -Original Message-
> From: Andrew Robinson [mailto:[EMAIL PROTECTED]
> Sent: Saturday, March 17, 2007 10:44 PM
> To: Thomas Colson
> Cc: r-help@stat.math.ethz.ch
> Subject: Re: [R] "Groups" in XYPLOT
>
> Hi again Thomas,
>
> ah, sorry, I should be more precise. Please construct a reproducible worked
> example that does not require us to download 7 Mb of data. You might also
> try the suggestions that I made and let us know if they worked for you.
>
> Cheers
>
> Andrew
>
> On Sat, Mar 17, 2007 at 10:37:46PM -0500, Thomas Colson wrote:
> > Thanks for the warning:
> > Here is the link to the datasets, rather large at 2 and 5 mb. Another
> > note is that one set has more datapoints than the other, don't know if
> > this can be done with xyplot.
> > http://www4.ncsu.edu/~tpcolson/coastcurvfreqs.txt
> > http://www4.ncsu.edu/~tpcolson/coastslopefreqs.txt
> >
> > Thomas Colson, PhD
> > North Carolina State University
> > Department of Forestry and Environmental Resources
> > (919)673-8023
> > [EMAIL PROTECTED]
> >
> > Schedule: www4.ncsu.edu/~tpcolson
> > -Original Message-
> > From: Andrew Robinson [mailto:[EMAIL PROTECTED]
> > Sent: Saturday, March 17, 2007 10:15 PM
> > To: Thomas Colson
> > Cc: r-help@stat.math.ethz.ch
> > Subject: Re: [R] "Groups" in XYPLOT
> >
> > Hi Thomas,
> >
> > sadly, the full code is not much help to us in the absence of the
> > data. Can I suggest that you construct a reproducible worked example
> > to help explain your question? For what it's worth I suspect that the
> > answer is that you need to join these datasets into one and theneitehr
> > use the groups argument, or the "+" protocol on the LHS of the plot
> formula.
> >
> > Cheers
> >
> > Andrew
> >
> > On Tue, Mar 27, 2007 at 04:51:55PM -0500, Thomas Colson wrote:
> > > I'm not sure I'm barking up the right tree here, but would I need to
> > > make use of groups to plot two separate datasets within ONE panel in
> > > xyplot? The desired end result is a single xy plot of two separate
> > > (but similar in values and ranges).
> > >
> > > Full code follows, xyplot code at bottom
> > >
> > >
> > >
> > >
> > >
> > > #Determine Frequencies
> > > ##coastal_slope
> > > #needs the maptools package to read ESRI grid
> > > require(maptools)
> > > #import the flow slope grid
> > > basin.map <- readAsciiGrid("C:/R_PLots/coastal_slp.asc",
> > > colname="slope") basin_slope <- (basin.map$slope)
Re: [R] "Groups" in XYPLOT
The make.groups function did the trick. Thank you so much for the seemingly
obvious solution!
piedfac <- read.table("C:/R_PLots/piedmontfacfreqs.txt", header=TRUE,
sep=",", na.strings="NA", dec=".", strip.white=TRUE)
coastfac <- read.table("C:/R_PLots/coastalfacfreqs.txt", header=TRUE,
sep=",", na.strings="NA", dec=".", strip.white=TRUE)
fac<-make.groups(piedfac,coastfac)
xyplot(fac$rank.PRank~fac$basin_area,groups=fac$which,scales=list(y=list(log
=TRUE,at=c(.0001,.001,.01,.1,1)),x=list(log=TRUE,at=c(10,100,1000,1,1000
00,100))),xlab="Drainage Area m^2",ylab="P(A>A*)")
Thomas Colson, PhD
North Carolina State University
Department of Forestry and Environmental Resources
(919)673-8023
[EMAIL PROTECTED]
Schedule: www4.ncsu.edu/~tpcolson
-Original Message-
From: Deepayan Sarkar [mailto:[EMAIL PROTECTED]
Sent: Saturday, March 17, 2007 11:19 PM
To: Thomas Colson
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] "Groups" in XYPLOT
On 3/17/07, Thomas Colson <[EMAIL PROTECTED]> wrote:
> Thanks for the warning:
> Here is the link to the datasets, rather large at 2 and 5 mb. Another
> note is that one set has more datapoints than the other, don't know if
> this can be done with xyplot.
As long as the two datasets have the same column name, you should be able to
use the 'make.groups' function to combine them. The resulting data frame
should have a column called 'which' identifying the origin.
Deepayan
__
R-help@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.
Re: [R] "Groups" in XYPLOT
Sorry 'bout that.
Here's about as simple as I can get it:
#First dataset
N <- 25
x <- round(rnorm(N),2)
y <- round(rnorm(N),2)
df <- data.frame(x = x, y = y)
#Plot it
xyplot(df$x~df$y,xlab="Test 1 Data",ylab="P(A>A*)")
#Second Dataset
N <- 20
x <- round(rnorm(N),2)
y <- round(rnorm(N),2)
df <- data.frame(x = x, y = y)
#How to get this in the same panel as plot 1?
xyplot(df$x~df$y,xlab="Test 2 Data",ylab="P(A>A*)")
As far as combining that data into one datasetI'm trying to plot both
datasets into one panel to show the differences in the plots of the
two...utilizing the coloring/symbology funtions...wouldn't combining them
null that capability?
Thanks for your replies!
Thomas Colson, PhD
North Carolina State University
Department of Forestry and Environmental Resources
(919)673-8023
[EMAIL PROTECTED]
Schedule: www4.ncsu.edu/~tpcolson
-Original Message-
From: Andrew Robinson [mailto:[EMAIL PROTECTED]
Sent: Saturday, March 17, 2007 10:44 PM
To: Thomas Colson
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] "Groups" in XYPLOT
Hi again Thomas,
ah, sorry, I should be more precise. Please construct a reproducible worked
example that does not require us to download 7 Mb of data. You might also
try the suggestions that I made and let us know if they worked for you.
Cheers
Andrew
On Sat, Mar 17, 2007 at 10:37:46PM -0500, Thomas Colson wrote:
> Thanks for the warning:
> Here is the link to the datasets, rather large at 2 and 5 mb. Another
> note is that one set has more datapoints than the other, don't know if
> this can be done with xyplot.
> http://www4.ncsu.edu/~tpcolson/coastcurvfreqs.txt
> http://www4.ncsu.edu/~tpcolson/coastslopefreqs.txt
>
> Thomas Colson, PhD
> North Carolina State University
> Department of Forestry and Environmental Resources
> (919)673-8023
> [EMAIL PROTECTED]
>
> Schedule: www4.ncsu.edu/~tpcolson
> -Original Message-
> From: Andrew Robinson [mailto:[EMAIL PROTECTED]
> Sent: Saturday, March 17, 2007 10:15 PM
> To: Thomas Colson
> Cc: r-help@stat.math.ethz.ch
> Subject: Re: [R] "Groups" in XYPLOT
>
> Hi Thomas,
>
> sadly, the full code is not much help to us in the absence of the
> data. Can I suggest that you construct a reproducible worked example
> to help explain your question? For what it's worth I suspect that the
> answer is that you need to join these datasets into one and theneitehr
> use the groups argument, or the "+" protocol on the LHS of the plot
formula.
>
> Cheers
>
> Andrew
>
> On Tue, Mar 27, 2007 at 04:51:55PM -0500, Thomas Colson wrote:
> > I'm not sure I'm barking up the right tree here, but would I need to
> > make use of groups to plot two separate datasets within ONE panel in
> > xyplot? The desired end result is a single xy plot of two separate
> > (but similar in values and ranges).
> >
> > Full code follows, xyplot code at bottom
> >
> >
> >
> >
> >
> > #Determine Frequencies
> > ##coastal_slope
> > #needs the maptools package to read ESRI grid
> > require(maptools)
> > #import the flow slope grid
> > basin.map <- readAsciiGrid("C:/R_PLots/coastal_slp.asc",
> > colname="slope") basin_slope <- (basin.map$slope) #read the slopes
> > into a dataframe
> > freqs<-as.data.frame(table(basin_slope))
> > #rank the frequencies based on each unique occerence, note, ranks
> > from
> > 1 to n
> > r<-rank(freqs$basin_slope)
> > n<-length(r)
> > #determing the probability, n+1 insures there is no 100%, 1-
> > reverses the order so #low slopes gets high probability of
> > exceedence z<-cbind(Rank = r, PRank = 1-(r/(n+1))) #attach the
> > probability to the table, result is high prob of exceed is in row
> > with low slope #and low probabibility is in row with high slope
> > freqs$rank<-z write.table(freqs, "C:/R_PLots/coastslopefreqs.txt",
> > sep=",", col.names=TRUE, row.names=TRUE, quote=TRUE, na="NA")
> >
> > ##coastal_curvature
> > #needs the maptools package to read ESRI grid
> > require(maptools)
> > #import the curvature grid
> > basin.map <- readAsciiGrid("C:/R_PLots/coastal_crv.asc",
> > colname="curv") basin_curv <- (basin.map$curv) #read the curvs into
> > a dataframe
> > freqs<-as.data.frame(table(basin_curv))
> > #rank the frequencies based on each unique occerence, note, ranks
> > from
> > 1 to n
> > r<-rank(freqs$basin_curv)
> > n<-length(r)
&g
Re: [R] "Groups" in XYPLOT
On 3/17/07, Thomas Colson <[EMAIL PROTECTED]> wrote: > Thanks for the warning: > Here is the link to the datasets, rather large at 2 and 5 mb. Another note > is that one set has more datapoints than the other, don't know if this can > be done with xyplot. As long as the two datasets have the same column name, you should be able to use the 'make.groups' function to combine them. The resulting data frame should have a column called 'which' identifying the origin. Deepayan __ R-help@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.
Re: [R] "Groups" in XYPLOT
Hi again Thomas,
ah, sorry, I should be more precise. Please construct a reproducible
worked example that does not require us to download 7 Mb of data. You
might also try the suggestions that I made and let us know if they
worked for you.
Cheers
Andrew
On Sat, Mar 17, 2007 at 10:37:46PM -0500, Thomas Colson wrote:
> Thanks for the warning:
> Here is the link to the datasets, rather large at 2 and 5 mb. Another note
> is that one set has more datapoints than the other, don't know if this can
> be done with xyplot.
> http://www4.ncsu.edu/~tpcolson/coastcurvfreqs.txt
> http://www4.ncsu.edu/~tpcolson/coastslopefreqs.txt
>
> Thomas Colson, PhD
> North Carolina State University
> Department of Forestry and Environmental Resources
> (919)673-8023
> [EMAIL PROTECTED]
>
> Schedule: www4.ncsu.edu/~tpcolson
> -Original Message-
> From: Andrew Robinson [mailto:[EMAIL PROTECTED]
> Sent: Saturday, March 17, 2007 10:15 PM
> To: Thomas Colson
> Cc: r-help@stat.math.ethz.ch
> Subject: Re: [R] "Groups" in XYPLOT
>
> Hi Thomas,
>
> sadly, the full code is not much help to us in the absence of the data. Can
> I suggest that you construct a reproducible worked example to help explain
> your question? For what it's worth I suspect that the answer is that you
> need to join these datasets into one and theneitehr use the groups argument,
> or the "+" protocol on the LHS of the plot formula.
>
> Cheers
>
> Andrew
>
> On Tue, Mar 27, 2007 at 04:51:55PM -0500, Thomas Colson wrote:
> > I'm not sure I'm barking up the right tree here, but would I need to
> > make use of groups to plot two separate datasets within ONE panel in
> > xyplot? The desired end result is a single xy plot of two separate
> > (but similar in values and ranges).
> >
> > Full code follows, xyplot code at bottom
> >
> >
> >
> >
> >
> > #Determine Frequencies
> > ##coastal_slope
> > #needs the maptools package to read ESRI grid
> > require(maptools)
> > #import the flow slope grid
> > basin.map <- readAsciiGrid("C:/R_PLots/coastal_slp.asc",
> > colname="slope") basin_slope <- (basin.map$slope) #read the slopes
> > into a dataframe
> > freqs<-as.data.frame(table(basin_slope))
> > #rank the frequencies based on each unique occerence, note, ranks from
> > 1 to n
> > r<-rank(freqs$basin_slope)
> > n<-length(r)
> > #determing the probability, n+1 insures there is no 100%, 1- reverses
> > the order so #low slopes gets high probability of exceedence
> > z<-cbind(Rank = r, PRank = 1-(r/(n+1))) #attach the probability to the
> > table, result is high prob of exceed is in row with low slope #and low
> > probabibility is in row with high slope freqs$rank<-z
> > write.table(freqs, "C:/R_PLots/coastslopefreqs.txt", sep=",",
> > col.names=TRUE, row.names=TRUE, quote=TRUE, na="NA")
> >
> > ##coastal_curvature
> > #needs the maptools package to read ESRI grid
> > require(maptools)
> > #import the curvature grid
> > basin.map <- readAsciiGrid("C:/R_PLots/coastal_crv.asc",
> > colname="curv") basin_curv <- (basin.map$curv) #read the curvs into a
> > dataframe
> > freqs<-as.data.frame(table(basin_curv))
> > #rank the frequencies based on each unique occerence, note, ranks from
> > 1 to n
> > r<-rank(freqs$basin_curv)
> > n<-length(r)
> > #determing the probability, n+1 insures there is no 100%, 1- reverses
> > the order so #low curvature gets high probability of exceedence
> > z<-cbind(Rank = r, PRank = 1-(r/(n+1))) #attach the probability to the
> > table, result is high prob of exceed is in row with low curv #and low
> > probabibility is in row with high curv freqs$rank<-z
> > write.table(freqs, "C:/R_PLots/coastcurvfreqs.txt", sep=",",
> > col.names=TRUE, row.names=TRUE, quote=TRUE, na="NA")
> >
> >
> >
> >
> >
> > ##Make XYPLOT and export to ps coastcurv <-
> > read.table("C:/R_PLots/coastcurvfreqs.txt", header=TRUE, sep=",",
> > na.strings="NA", dec=".", strip.white=TRUE)
> > xyplot(coastcurv$rank.PRank~coastcurv$basin_curv,scales=list(y=list(lo
> > g=TRUE
> > ,at=c(.0001,.001,.01,.1,1)),x=list(log=TRUE,at=c(0.0001,0.001,0.01,0.1
> > ,1,10)
> > )),xlab="Curvature",ylab="P(C>C*)")
> > dev.copy2eps(file="C:/R_PLots/coastcurv_cad.eps",
Re: [R] "Groups" in XYPLOT
Thanks for the warning:
Here is the link to the datasets, rather large at 2 and 5 mb. Another note
is that one set has more datapoints than the other, don't know if this can
be done with xyplot.
http://www4.ncsu.edu/~tpcolson/coastcurvfreqs.txt
http://www4.ncsu.edu/~tpcolson/coastslopefreqs.txt
Thomas Colson, PhD
North Carolina State University
Department of Forestry and Environmental Resources
(919)673-8023
[EMAIL PROTECTED]
Schedule: www4.ncsu.edu/~tpcolson
-Original Message-
From: Andrew Robinson [mailto:[EMAIL PROTECTED]
Sent: Saturday, March 17, 2007 10:15 PM
To: Thomas Colson
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] "Groups" in XYPLOT
Hi Thomas,
sadly, the full code is not much help to us in the absence of the data. Can
I suggest that you construct a reproducible worked example to help explain
your question? For what it's worth I suspect that the answer is that you
need to join these datasets into one and theneitehr use the groups argument,
or the "+" protocol on the LHS of the plot formula.
Cheers
Andrew
On Tue, Mar 27, 2007 at 04:51:55PM -0500, Thomas Colson wrote:
> I'm not sure I'm barking up the right tree here, but would I need to
> make use of groups to plot two separate datasets within ONE panel in
> xyplot? The desired end result is a single xy plot of two separate
> (but similar in values and ranges).
>
> Full code follows, xyplot code at bottom
>
>
>
>
>
> #Determine Frequencies
> ##coastal_slope
> #needs the maptools package to read ESRI grid
> require(maptools)
> #import the flow slope grid
> basin.map <- readAsciiGrid("C:/R_PLots/coastal_slp.asc",
> colname="slope") basin_slope <- (basin.map$slope) #read the slopes
> into a dataframe
> freqs<-as.data.frame(table(basin_slope))
> #rank the frequencies based on each unique occerence, note, ranks from
> 1 to n
> r<-rank(freqs$basin_slope)
> n<-length(r)
> #determing the probability, n+1 insures there is no 100%, 1- reverses
> the order so #low slopes gets high probability of exceedence
> z<-cbind(Rank = r, PRank = 1-(r/(n+1))) #attach the probability to the
> table, result is high prob of exceed is in row with low slope #and low
> probabibility is in row with high slope freqs$rank<-z
> write.table(freqs, "C:/R_PLots/coastslopefreqs.txt", sep=",",
> col.names=TRUE, row.names=TRUE, quote=TRUE, na="NA")
>
> ##coastal_curvature
> #needs the maptools package to read ESRI grid
> require(maptools)
> #import the curvature grid
> basin.map <- readAsciiGrid("C:/R_PLots/coastal_crv.asc",
> colname="curv") basin_curv <- (basin.map$curv) #read the curvs into a
> dataframe
> freqs<-as.data.frame(table(basin_curv))
> #rank the frequencies based on each unique occerence, note, ranks from
> 1 to n
> r<-rank(freqs$basin_curv)
> n<-length(r)
> #determing the probability, n+1 insures there is no 100%, 1- reverses
> the order so #low curvature gets high probability of exceedence
> z<-cbind(Rank = r, PRank = 1-(r/(n+1))) #attach the probability to the
> table, result is high prob of exceed is in row with low curv #and low
> probabibility is in row with high curv freqs$rank<-z
> write.table(freqs, "C:/R_PLots/coastcurvfreqs.txt", sep=",",
> col.names=TRUE, row.names=TRUE, quote=TRUE, na="NA")
>
>
>
>
>
> ##Make XYPLOT and export to ps coastcurv <-
> read.table("C:/R_PLots/coastcurvfreqs.txt", header=TRUE, sep=",",
> na.strings="NA", dec=".", strip.white=TRUE)
> xyplot(coastcurv$rank.PRank~coastcurv$basin_curv,scales=list(y=list(lo
> g=TRUE
> ,at=c(.0001,.001,.01,.1,1)),x=list(log=TRUE,at=c(0.0001,0.001,0.01,0.1
> ,1,10)
> )),xlab="Curvature",ylab="P(C>C*)")
> dev.copy2eps(file="C:/R_PLots/coastcurv_cad.eps", width=8.0,
> height=8.0,
> pointsize=10)
>
>
> How to get this in the first plot graphic?
>
> coastslope <- read.table("C:/R_PLots/coastslopefreqs.txt",
> header=TRUE, sep=",", na.strings="NA", dec=".", strip.white=TRUE)
> xyplot(coastslope$rank.PRank~coastslope$basin_slope,scales=list(y=list
> (log=T
> RUE,at=c(.0001,.001,.01,.1,1)),x=list(log=TRUE,at=c(0.0001,0.001,0.01,
> 0.1,1,
> 10))),xlab="Slope",ylab="P(S>S*)")
> dev.copy2eps(file="C:/R_PLots/coastslope_cad.eps", width=8.0,
> height=8.0,
> pointsize=10)
>
> Thomas Colson, PhD
> North Carolina State University
> Department of Forestry and Environmental Resources
>
> _
Re: [R] "Groups" in XYPLOT
Hi Thomas,
sadly, the full code is not much help to us in the absence of the
data. Can I suggest that you construct a reproducible worked example
to help explain your question? For what it's worth I suspect that the
answer is that you need to join these datasets into one and theneitehr
use the groups argument, or the "+" protocol on the LHS of the plot
formula.
Cheers
Andrew
On Tue, Mar 27, 2007 at 04:51:55PM -0500, Thomas Colson wrote:
> I'm not sure I'm barking up the right tree here, but would I need to make
> use of groups to plot two separate datasets within ONE panel in xyplot? The
> desired end result is a single xy plot of two separate (but similar in
> values and ranges).
>
> Full code follows, xyplot code at bottom
>
>
>
>
>
> #Determine Frequencies
> ##coastal_slope
> #needs the maptools package to read ESRI grid
> require(maptools)
> #import the flow slope grid
> basin.map <- readAsciiGrid("C:/R_PLots/coastal_slp.asc", colname="slope")
> basin_slope <- (basin.map$slope)
> #read the slopes into a dataframe
> freqs<-as.data.frame(table(basin_slope))
> #rank the frequencies based on each unique occerence, note, ranks from 1 to
> n
> r<-rank(freqs$basin_slope)
> n<-length(r)
> #determing the probability, n+1 insures there is no 100%, 1- reverses the
> order so
> #low slopes gets high probability of exceedence
> z<-cbind(Rank = r, PRank = 1-(r/(n+1)))
> #attach the probability to the table, result is high prob of exceed is in
> row with low slope
> #and low probabibility is in row with high slope
> freqs$rank<-z
> write.table(freqs, "C:/R_PLots/coastslopefreqs.txt", sep=",",
> col.names=TRUE, row.names=TRUE, quote=TRUE, na="NA")
>
> ##coastal_curvature
> #needs the maptools package to read ESRI grid
> require(maptools)
> #import the curvature grid
> basin.map <- readAsciiGrid("C:/R_PLots/coastal_crv.asc", colname="curv")
> basin_curv <- (basin.map$curv)
> #read the curvs into a dataframe
> freqs<-as.data.frame(table(basin_curv))
> #rank the frequencies based on each unique occerence, note, ranks from 1 to
> n
> r<-rank(freqs$basin_curv)
> n<-length(r)
> #determing the probability, n+1 insures there is no 100%, 1- reverses the
> order so
> #low curvature gets high probability of exceedence
> z<-cbind(Rank = r, PRank = 1-(r/(n+1)))
> #attach the probability to the table, result is high prob of exceed is in
> row with low curv
> #and low probabibility is in row with high curv
> freqs$rank<-z
> write.table(freqs, "C:/R_PLots/coastcurvfreqs.txt", sep=",", col.names=TRUE,
> row.names=TRUE, quote=TRUE, na="NA")
>
>
>
>
>
> ##Make XYPLOT and export to ps
> coastcurv <- read.table("C:/R_PLots/coastcurvfreqs.txt", header=TRUE,
> sep=",", na.strings="NA", dec=".", strip.white=TRUE)
> xyplot(coastcurv$rank.PRank~coastcurv$basin_curv,scales=list(y=list(log=TRUE
> ,at=c(.0001,.001,.01,.1,1)),x=list(log=TRUE,at=c(0.0001,0.001,0.01,0.1,1,10)
> )),xlab="Curvature",ylab="P(C>C*)")
> dev.copy2eps(file="C:/R_PLots/coastcurv_cad.eps", width=8.0, height=8.0,
> pointsize=10)
>
>
> How to get this in the first plot graphic?
>
> coastslope <- read.table("C:/R_PLots/coastslopefreqs.txt", header=TRUE,
> sep=",", na.strings="NA", dec=".", strip.white=TRUE)
> xyplot(coastslope$rank.PRank~coastslope$basin_slope,scales=list(y=list(log=T
> RUE,at=c(.0001,.001,.01,.1,1)),x=list(log=TRUE,at=c(0.0001,0.001,0.01,0.1,1,
> 10))),xlab="Slope",ylab="P(S>S*)")
> dev.copy2eps(file="C:/R_PLots/coastslope_cad.eps", width=8.0, height=8.0,
> pointsize=10)
>
> Thomas Colson, PhD
> North Carolina State University
> Department of Forestry and Environmental Resources
>
> __
> R-help@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.
--
Andrew Robinson
Department of Mathematics and StatisticsTel: +61-3-8344-9763
University of Melbourne, VIC 3010 Australia Fax: +61-3-8344-4599
http://www.ms.unimelb.edu.au/~andrewpr
http://blogs.mbs.edu/fishing-in-the-bay/
__
R-help@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.
