#this is what I came up with
diff.temp <- function(x, a, b,use="pairwise.complete.obs")
{
library(zoo)
library(chron)
na.method <- pmatch(use, c("all.obs", "complete.obs",
"pairwise.complete.obs"))
par(mfrow=c(2,1))
f <- colnames(x)
d <- (x[,a]-x[,b])
plot(d, main= paste(f[a], "-", f[b]))
plot(density(na.omit(coredata(d))), main="")
}
On Sat, Jul 19, 2008 at 12:04 PM, stephen sefick <[EMAIL PROTECTED]> wrote:
> #this is my little function that I would like to use the column names of
> the x and y arguments in the function. I would like it to read
> # site1-site2 how would I do this
> diff.temp <- function(x, y ,use="pairwise.complete.obs")
> {
> na.method <- pmatch(use, c("all.obs", "complete.obs",
> "pairwise.complete.obs"))
> par(mfrow=c(2,1))
> d <- (x-y)
> plot(d, main="paste(colnames(x))-paste(colnames(y))")
> plot(density(na.omit(coredata(d))))
> }
>
> --
> Let's not spend our time and resources thinking about things that are so
> little or so large that all they really do for us is puff us up and make us
> feel like gods. We are mammals, and have not exhausted the annoying little
> problems of being mammals.
>
> -K. Mullis
>
--
Let's not spend our time and resources thinking about things that are so
little or so large that all they really do for us is puff us up and make us
feel like gods. We are mammals, and have not exhausted the annoying little
problems of being mammals.
-K. Mullis
[[alternative HTML version deleted]]
______________________________________________
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.
