I have a dataset about gene expression, which has n=62 observations and p=2000 variables (human genes). I also have one Y binary variable, which denotes as 2==Tumor and 1==Normal.

I need to create one-way ANOVA between y and each X. So there are 2000 ANOVAs. I need to extract each p-value. So I will create a vector that has 2000 components, here are 2000 p values. Then I need to order these p-values. The reason that I do this is I will do false discovery rate (FDR) (Benjamini and Hochberg) later.

```
one.way <- aov(response ~ X1, data = df)
summary(one.way)[[1]][["Pr(>F)"]][1]
one.way <- aov(response ~ X2, data = df)
summary(one.way)[[1]][["Pr(>F)"]][1]
one.way <- aov(response ~ X3, data = df)
summary(one.way)[[1]][["Pr(>F)"]][1]
```

I need to do response ~ X1 until response ~ X2000. How to do that and order these p-values and its corresponding X?

By the way, I omit uploading my dataset. Because otherwise the website will be full of my data.

```
dput(head(df,4))
```