[R] Building a formula using paste
john seers (IFR)
john.seers at bbsrc.ac.uk
Fri Feb 1 13:06:35 CET 2008
Thank you very much. That problem has been niggling me for some time.
I slotted in your code and it worked. I just need to spend a bit of time
understanding it ...
Thanks again.
John Seers
---
-----Original Message-----
From: Bill.Venables at csiro.au [mailto:Bill.Venables at csiro.au]
Sent: 01 February 2008 09:57
To: Bill.Venables at csiro.au; john seers (IFR); r-help at stat.math.ethz.ch
Subject: RE: [R] Building a formula using paste
Please disregard my previous reply. Now that I have read your question,
all becomes clear.
To do what you want to do is a bit tricky. Here is one way
f1 <- as.name("F1")
fm <- eval(bquote(aov(.(f1) ~ sensoryTerm, data = vdata)))
Them fm is the fitted model object, suitably formed.
E&OE, of course!
Bill Venables
CSIRO Laboratories
PO Box 120, Cleveland, 4163
AUSTRALIA
Office Phone (email preferred): +61 7 3826 7251 Fax (if absolutely
necessary): +61 7 3826 7304
Mobile: +61 4 8819 4402
Home Phone: +61 7 3286 7700
mailto:Bill.Venables at csiro.au
http://www.cmis.csiro.au/bill.venables/
-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
On Behalf Of Venables, Bill (CMIS, Cleveland)
Sent: Friday, 1 February 2008 7:45 PM
To: john.seers at bbsrc.ac.uk; r-help at stat.math.ethz.ch
Subject: Re: [R] Building a formula using paste
use as.formula() rather than just formula().
Bill Venables
CSIRO Laboratories
PO Box 120, Cleveland, 4163
AUSTRALIA
Office Phone (email preferred): +61 7 3826 7251 Fax (if absolutely
necessary): +61 7 3826 7304
Mobile: +61 4 8819 4402
Home Phone: +61 7 3286 7700
mailto:Bill.Venables at csiro.au
http://www.cmis.csiro.au/bill.venables/
-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
On Behalf Of john seers (IFR)
Sent: Friday, 1 February 2008 7:31 PM
To: r-help at stat.math.ethz.ch
Subject: [R] Building a formula using paste
Hi All
I have looked but cannot find an answer to this.
I want to build a formula dynamically (in a function), using for example
paste, and using it in a model:
fr<-"F1"
f1<-formula(paste(fr, "~ SensoryTerm"))
m1<-aov(f1, data=vdata)
So this is the equivalent of m1<-aov(F1 ~ SensoryTerm, data=vdata)
This works fine but the problem is the formula appears as "f1" in the
summary(m1) of the model and not "F1 ~ SensoryTerm". With many models
this can be a bit confusing ....
Is there a way to code this so the formula appears in the model summary?
Regards
John Seers
> sessionInfo()
R version 2.6.1 (2007-11-26)
i386-pc-mingw32
locale:
LC_COLLATE=English_United Kingdom.1252;LC_CTYPE=English_United
Kingdom.1252;LC_MONETARY=English_United
Kingdom.1252;LC_NUMERIC=C;LC_TIME=English_United Kingdom.1252
attached base packages:
[1] stats graphics grDevices datasets utils methods base
other attached packages:
[1] pixmap_0.4-7 gplots_2.3.2 gdata_2.3.1 gtools_2.4.0
rcom_1.5-2.2
[6] nnet_7.2-38 e1071_1.5-17 class_7.2-38 tree_1.0-26 mgcv_1.3-29
[11] MASS_7.2-38 car_1.2-7 RODBC_1.2-2 RWinEdt_1.7-9
>
---
______________________________________________
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.
______________________________________________
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.
More information about the R-help
mailing list