p.adjustの使い方

Question

以下のコードは、2つの行列を入力として使用してp-valueを返すt-test（un-paired）を実行しますが、p.adjusted値を生成するために必要な正確なコードを私に表示できます。FDR修正/調整p.values（p.adjust）

m1 <-go_samp_matrix_data[grep(paste(input$mychooser2$left, collapse='|'), rownames(go_samp_matrix_data), ignore.case=TRUE),] 
m2 <-go_samp_matrix_data[grep(paste(input$mychooser2$right, collapse='|'), rownames(go_samp_matrix_data), ignore.case=TRUE),] 

ttestmat1<-sapply(seq(ncol(m1)), function(x) f(m1[,x], m2[,x])) 

f <- function(x,y){ 
    test <- t.test(x,y, paired=FALSE) 
    out <- data.frame(pval = sprintf("%.3f", test$p.value)) 
    return(out) 
} 

Answer 1

あなたはp値のベクトルを必要とする、p.vectを言います。

p.adjust(p.vect, "fdr")