Re: [R] automating a script to read a file
Petr, R- Users
The scheme you provided yesterday worked very well. I am now trying to add
more information to the graph(s) to make them more informative for my
purposes.
Essentially, I am trying to use a double Y axis graphic that will include
the skew normal curves that I wrote about yesterday and a histogram.
The graph is very similar to the one produced by the following:
library(latticeExtra)
doubleYScale(histogram(x), densityplot(x), use.style = FALSE)
However, I'd like to substitute the following syntax:
doubleYScale(hist(mydata, breaks=20, prob=T, xlim=c(-100, 2000),
plot.spdf(x), use.style=FALSE)
This does not work as doubleYScale expects histogram and densityplot, and
I'd like to use the plot.spdf routine in its place.
Any Suggestions would be greatly appreciated.
I'm working with
R 2.15.0 (2012-03-12)
Platform i386-pc-mingw32/ie86 (32-bit)
Thanks
Steve
Steve Friedman Ph. D.
Ecologist / Spatial Statistical Analyst
Everglades and Dry Tortugas National Park
950 N Krome Ave (3rd Floor)
Homestead, Florida 33034
steve_fried...@nps.gov
Office (305) 224 - 4282
Fax (305) 224 - 4147
Petr Savicky
savi...@cs.cas.c
z To
Sent by: r-help@r-project.org
r-help-bounces@r- cc
project.org
Subject
Re: [R] automating a script to read
04/23/2012 04:42 a file
PM
On Mon, Apr 23, 2012 at 04:02:45PM -0400, steve_fried...@nps.gov wrote:
Hi,
The following script (which I did not develop) is used to calculate and
plot a skewed normal curve. The script currently requires the user to
input six parameters, rather than reading these directly from a file.
I've been spinning wheels here, trying to figure out how to modify the
script to automate it. I have four data sets, each in excess of 300
records that I need to process.
My initial thoughts were to use the lapply and use a pdf graphic device
to
capture the plots to do this, but my R programming skills are too limited
to determine how to best accomplish this.
Hi.
If you read the parameters from a file and put them to a matrix,
then all the plots may be produced using a loop like the following.
#some parameters
p - matrix(1:18, nrow=3, ncol=6)
for (i in 1:nrow(p)) {
plot.spdf(p[i, 1], p[i, 2], p[i, 3], p[i, 4], p[i, 5], p[i, 6])
readline(press Enter to continue)
}
If you use pdf() for sending the graphics to a file, then remove
the readline command.
Hope this helps.
Petr Savicky.
__
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.
__
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.
Re: [R] automating a script to read a file
On Tue, Apr 24, 2012 at 2:09 PM, steve_fried...@nps.gov wrote:
library(latticeExtra)
doubleYScale(hist(mydata, breaks=20, prob=T, xlim=c(-100, 2000),
plot.spdf(x), use.style=FALSE)
This does not work as doubleYScale expects histogram and densityplot, and
I'd like to use the plot.spdf routine in its place.
Maybe it doesn't work because ?hist is base graphics, not the same as
?histogram in lattice. You did not provide plot.spdf or any data, so
this just a place to start type response.
Cheers
Any Suggestions would be greatly appreciated.
I'm working with
R 2.15.0 (2012-03-12)
Platform i386-pc-mingw32/ie86 (32-bit)
Thanks
Steve
Steve Friedman Ph. D.
Ecologist / Spatial Statistical Analyst
Everglades and Dry Tortugas National Park
950 N Krome Ave (3rd Floor)
Homestead, Florida 33034
steve_fried...@nps.gov
Office (305) 224 - 4282
Fax (305) 224 - 4147
Petr Savicky
savi...@cs.cas.c
z To
Sent by: r-help@r-project.org
r-help-bounces@r- cc
project.org
Subject
Re: [R] automating a script to read
04/23/2012 04:42 a file
PM
On Mon, Apr 23, 2012 at 04:02:45PM -0400, steve_fried...@nps.gov wrote:
Hi,
The following script (which I did not develop) is used to calculate and
plot a skewed normal curve. The script currently requires the user to
input six parameters, rather than reading these directly from a file.
I've been spinning wheels here, trying to figure out how to modify the
script to automate it. I have four data sets, each in excess of 300
records that I need to process.
My initial thoughts were to use the lapply and use a pdf graphic device
to
capture the plots to do this, but my R programming skills are too limited
to determine how to best accomplish this.
Hi.
If you read the parameters from a file and put them to a matrix,
then all the plots may be produced using a loop like the following.
#some parameters
p - matrix(1:18, nrow=3, ncol=6)
for (i in 1:nrow(p)) {
plot.spdf(p[i, 1], p[i, 2], p[i, 3], p[i, 4], p[i, 5], p[i, 6])
readline(press Enter to continue)
}
If you use pdf() for sending the graphics to a file, then remove
the readline command.
Hope this helps.
Petr Savicky.
__
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.
__
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.
__
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.
Re: [R] automating a script to read a file
On Tue, Apr 24, 2012 at 6:00 PM, Steve Friedman skfgla...@gmail.com wrote:
The original post does contain the function plot.spdf and some sample data.
That's great. And who's responsibility is it to make sure we have
context for THIS post ? This is a mailing list, not a forum and MY
mail client doesn't keep track of YOUR postings.
BTW, you are the one seeking our help, so even if you think it was a
picky request, a proper response would be to simply post your
function+data set. Not send us looking for some original post.
I thought I also made it clear that I understood that hist and histogram are
not the same.
No you didn't. But now IT IS clear you don't realize grid and base
graphics don't mix.
I am looking for a way to use the alternatives in a function similar to but
not identical to doubleYScale.
LatticeExtra functions expect any lattice plots, including local (user
defined) panel/plot functions, if only you would have posted plot.spdf
maybe you would have gotten a solution by now...
If plot.spdf is actually a base graphics plotter and you just stumbled
upon DoubleYscale somehow, ?par ?axis hist(...,add=T) are all
available to you.
Best
Steve
On Apr 24, 2012 7:45 PM, ilai ke...@math.montana.edu wrote:
On Tue, Apr 24, 2012 at 2:09 PM, steve_fried...@nps.gov wrote:
library(latticeExtra)
doubleYScale(hist(mydata, breaks=20, prob=T, xlim=c(-100, 2000),
plot.spdf(x), use.style=FALSE)
This does not work as doubleYScale expects histogram and densityplot,
and
I'd like to use the plot.spdf routine in its place.
Maybe it doesn't work because ?hist is base graphics, not the same as
?histogram in lattice. You did not provide plot.spdf or any data, so
this just a place to start type response.
Cheers
Any Suggestions would be greatly appreciated.
I'm working with
R 2.15.0 (2012-03-12)
Platform i386-pc-mingw32/ie86 (32-bit)
Thanks
Steve
Steve Friedman Ph. D.
Ecologist / Spatial Statistical Analyst
Everglades and Dry Tortugas National Park
950 N Krome Ave (3rd Floor)
Homestead, Florida 33034
steve_fried...@nps.gov
Office (305) 224 - 4282
Fax (305) 224 - 4147
Petr Savicky
savi...@cs.cas.c
z
To
Sent by: r-help@r-project.org
r-help-bounces@r-
cc
project.org
Subject
Re: [R] automating a script to
read
04/23/2012 04:42 a file
PM
On Mon, Apr 23, 2012 at 04:02:45PM -0400, steve_fried...@nps.gov wrote:
Hi,
The following script (which I did not develop) is used to calculate and
plot a skewed normal curve. The script currently requires the user to
input six parameters, rather than reading these directly from a file.
I've been spinning wheels here, trying to figure out how to modify the
script to automate it. I have four data sets, each in excess of 300
records that I need to process.
My initial thoughts were to use the lapply and use a pdf graphic
device
to
capture the plots to do this, but my R programming skills are too
limited
to determine how to best accomplish this.
Hi.
If you read the parameters from a file and put them to a matrix,
then all the plots may be produced using a loop like the following.
#some parameters
p - matrix(1:18, nrow=3, ncol=6)
for (i in 1:nrow(p)) {
plot.spdf(p[i, 1], p[i, 2], p[i, 3], p[i, 4], p[i, 5], p[i, 6])
readline(press Enter to continue)
}
If you use pdf() for sending the graphics to a file, then remove
the readline command.
Hope this helps.
Petr Savicky.
__
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.
__
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.
__
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.
__
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.
Re: [R] automating a script to read a file
On Mon, Apr 23, 2012 at 04:02:45PM -0400, steve_fried...@nps.gov wrote:
Hi,
The following script (which I did not develop) is used to calculate and
plot a skewed normal curve. The script currently requires the user to
input six parameters, rather than reading these directly from a file.
I've been spinning wheels here, trying to figure out how to modify the
script to automate it. I have four data sets, each in excess of 300
records that I need to process.
My initial thoughts were to use the lapply and use a pdf graphic device to
capture the plots to do this, but my R programming skills are too limited
to determine how to best accomplish this.
Hi.
If you read the parameters from a file and put them to a matrix,
then all the plots may be produced using a loop like the following.
#some parameters
p - matrix(1:18, nrow=3, ncol=6)
for (i in 1:nrow(p)) {
plot.spdf(p[i, 1], p[i, 2], p[i, 3], p[i, 4], p[i, 5], p[i, 6])
readline(press Enter to continue)
}
If you use pdf() for sending the graphics to a file, then remove
the readline command.
Hope this helps.
Petr Savicky.
__
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.
Re: [R] automating a script to read a file
Petr,
Thank you very much this works. A little more tweaking and I'll have what
I need.
Thanks
Steve Friedman Ph. D.
Ecologist / Spatial Statistical Analyst
Everglades and Dry Tortugas National Park
950 N Krome Ave (3rd Floor)
Homestead, Florida 33034
steve_fried...@nps.gov
Office (305) 224 - 4282
Fax (305) 224 - 4147
Petr Savicky
savi...@cs.cas.c
z To
Sent by: r-help@r-project.org
r-help-bounces@r- cc
project.org
Subject
Re: [R] automating a script to read
04/23/2012 04:42 a file
PM
On Mon, Apr 23, 2012 at 04:02:45PM -0400, steve_fried...@nps.gov wrote:
Hi,
The following script (which I did not develop) is used to calculate and
plot a skewed normal curve. The script currently requires the user to
input six parameters, rather than reading these directly from a file.
I've been spinning wheels here, trying to figure out how to modify the
script to automate it. I have four data sets, each in excess of 300
records that I need to process.
My initial thoughts were to use the lapply and use a pdf graphic device
to
capture the plots to do this, but my R programming skills are too limited
to determine how to best accomplish this.
Hi.
If you read the parameters from a file and put them to a matrix,
then all the plots may be produced using a loop like the following.
#some parameters
p - matrix(1:18, nrow=3, ncol=6)
for (i in 1:nrow(p)) {
plot.spdf(p[i, 1], p[i, 2], p[i, 3], p[i, 4], p[i, 5], p[i, 6])
readline(press Enter to continue)
}
If you use pdf() for sending the graphics to a file, then remove
the readline command.
Hope this helps.
Petr Savicky.
__
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.
__
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.
