micfalb-r wrote:
>
>
> > p.fun <- function(arg) {
> two_way_anova <- aov(arg ~ age * treatment, data = example.df)
> two_way_sum <- summary(two_way_anova)
> p_values <- two_way_sum[[1]]$"Pr(>F)"[1:3]
> return(p_values)
> }
>
> Unfortunately my setup seems to be flawed as I'm not capable to call my
> function:
> > p.fun(gene1)
> Error in eval(expr, envir, enclos) : object 'gene1' not found
> > p.fun("gene1")
> Error in model.frame.default(formula = arg ~ age * treatment, data =
> example.df, :
> variable lengths differ (found for 'age')
> and provide commented, minimal, self-contained, reproducible code.
>
Try this for your function
p.fun <- function(arg) {
lhs <- get(arg, envir=as.environment(example.df))
two_way_anova <- aov(lhs ~ age * treatment, data = example.df)
two_way_sum <- summary(two_way_anova)
p_values <- two_way_sum[[1]]$"Pr(>F)"[1:3]
return(p_values)
}
Use with
p.fun("gene1")
There are surely more elegant methods. E.g. no quotes in the argument of the
p.fun call.
/Berend
--
View this message in context:
http://r.789695.n4.nabble.com/Function-to-crop-p-values-from-multiple-Anovas-tp3400271p3400434.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
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.
