[R] a question of substitute

Gabor Grothendieck ggrothendieck at gmail.com
Wed Jan 10 18:03:22 CET 2007


Looks like oneway.test has been changed for R 2.5.0.

Paste the code in this file:

   https://svn.r-project.org/R/trunk/src/library/stats/R/oneway.test.R

into your session.  Then fun.2 from your post will work without
the workarounds I posted:

   fun.2(values ~ group)


On 1/9/07, Adrian Dusa <dusa.adrian at gmail.com> wrote:
> On Tuesday 09 January 2007 15:14, Gabor Grothendieck wrote:
> > oneway.test is using substitute on its arguments so its literally
> > getting formula rather than the value of formula.
>
> Ah-haa... I understand now. Thanks for the tips, they both work as expected.
> Best,
> Adrian
>
> > Try these:
> >
> > fun.3 <- function(formula) {
> >       mc <- match.call()
> >       mc[[1]] <- as.name("oneway.test")
> >       eval.parent(mc)
> > }
> > fun.3(values ~ group)
> >
> > fun.4 <- function(formula) {
> >       do.call(oneway.test, list(formula))
> > }
> > fun.4(values ~ group)
> >
> > On 1/9/07, Adrian Dusa <dusa.adrian at gmail.com> wrote:
> > > Hi all,
> > >
> > > I want to write a wrapper for an analysis of variance and I face a
> > > curious problem. Here are two different wrappers:
> > >
> > > fun.1 <- function(formula) {
> > >    summary(aov(formula))
> > > }
> > >
> > > fun.2 <- function(formula) {
> > >    oneway.test(formula)
> > > }
> > >
> > > values <- c(15, 8, 17, 7, 26, 12, 8, 11, 16, 9, 16,
> > >            24, 20, 19, 9, 17, 11, 8, 15, 6, 14)
> > > group <- rep(1:3, each=7)
> > >
> > > # While the first wrapper works just fine:
> > > fun.1(values ~ group)
> > >
> > > # the second throws an error:
> > > fun.2(values ~ group)
> > > Error in substitute(formula)[[2]] : object is not subsettable
> > >
> > > ###
> > >
> > > I also tried binding the two vectors in a data.frame, with no avail.
> > > I did find a hack, creating two new vectors inside the function and
> > > creating a fresh formula, so I presume this has something to do with
> > > environments.
> > >
> > > Could anybody give me a hint on this?
> > > Thank you,
> > > Adrian
> > >
> > > --
> > > Adrian Dusa
> > > Romanian Social Data Archive
> > > 1, Schitu Magureanu Bd
> > > 050025 Bucharest sector 5
> > > Romania
> > > Tel./Fax: +40 21 3126618 \
> > >          +40 21 3120210 / int.101
> > >
> > > ______________________________________________
> > > 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.
>
> --
> Adrian Dusa
> Arhiva Romana de Date Sociale
> Bd. Schitu Magureanu nr.1
> 050025 Bucuresti sectorul 5
> Romania
> Tel./Fax: +40 21 3126618 \
>          +40 21 3120210 / int.101
>



More information about the R-help mailing list