I have some bulk-RNA sequencing data that I need to do differential expression significance testing on. I have two conditions, WT and KO, with two replicates each, giving me a dataframe that looks like the following (the columns are in counts):
WT1 WT2 KO1 KO2
gene1 1.3 1.23 3.42 3.45
gene2 2.6 2.54 1.22 1.21
gene3 5.54 2.32 1.21 1.10
My questions are, how do I get a column on the right with a p-value for each gene so that I can construct a Volcano plot of the data? Basically, what statistical test do I need to use to generate that column, and what function do I use in R to do so? I'm sorry if this isn't technically a question that I'm supposed to ask here, but frankly I didn't know where else to post. Thanks in advance!