差次的発現有意性テストを行う必要があるバルク RNA シーケンス データがいくつかあります。WT と KO の 2 つの条件があり、それぞれ 2 つのレプリケートがあり、次のようなデータフレームが得られます (列はカウントです)。
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
私の質問は、データの火山プロットを作成できるように、各遺伝子の p 値を含む右側の列を取得するにはどうすればよいですか? 基本的に、その列を生成するにはどの統計テストを使用する必要があり、そのために R でどの関数を使用すればよいでしょうか? これが技術的にここで質問すべき質問ではない場合は申し訳ありませんが、率直に言って、他のどこに投稿すればよいかわかりませんでした。前もって感謝します！