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
Starting with data from a microarray experiment and I would like to analyse the
influence of two factors (age, treatment) on gene expression.
Looking through the r-help archives and the web I tried the following:
I put my data in a dataframe similar to this one:
> example.df <- as.data.frame(mat
2 matches
Mail list logo
