Re: [R] extracting colnames to label plots in a function

From: stephen sefick <ssefick_at_gmail.com>
Date: Sat, 19 Jul 2008 12:50:59 -0400

h <- structure(c(1L, 2L, 3L, 4L, NA, 6L, 7L, 8L, NA, 10L, 1L, NA, 8L, 7L, 6L, 5L, 4L, 3L, 2L, 1L), .Dim = c(10L, 2L), .Dimnames = list(

    c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10"), c("site1",     "site2")), index = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10), class = "zoo")

diff.temp <- function(x, a, b,use="pairwise.complete.obs") {

    library(zoo)
    library(chron)
    na.method <- match.arg(use, c("all.obs", "complete.obs",     "pairwise.complete.obs"))
    opar <- par
    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="") par(opar)
}

give a NULL after execution
what is the problem?

On Sat, Jul 19, 2008 at 12:42 PM, stephen sefick <ssefick_at_gmail.com> wrote:

> I can fix the par settings. I am new to function writing. I would like to
> use the na.method in the
> d<- (x-y)
>
> how?
>
> I don't know what a driver is sorry for my ignorance
>
>
> On Sat, Jul 19, 2008 at 12:36 PM, Gabor Grothendieck <
> ggrothendieck_at_gmail.com> wrote:
>
>> See ?paste and the collapse argument, in particular:
>>
>> plot(d, main = paste(paste(colnames(x), collapse = " "),
>> paste(colnames(y), collapse = " "), sep = " - "))
>>
>> Also your function sets na.method but never uses it and leaves the par
>> settings
>> changed afterwards. See ?par. It could also benefit from the use of
>> ?match.arg
>>
>> Also please include drivers that call the posted function so one can
>> run them in
>> a reproducible manner.
>>
>>
>> On Sat, Jul 19, 2008 at 12:04 PM, stephen sefick <ssefick_at_gmail.com>
>> 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
>> >
>> > [[alternative HTML version deleted]]
>> >
>> > ______________________________________________
>> > R-help_at_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.
>> >
>>
>
>
>
> --
> 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_at_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.
Received on Sat 19 Jul 2008 - 16:53:46 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sat 19 Jul 2008 - 17:32:00 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive