r - Homoscedastic 2 tailed paired t-test with replicates -

i need calculate p-values of 7 treatments vs control 3 replicates, i'm writing function calculating different things (regulation, ic50, ld50 , other parameters i'm interesting into) , merge df, , i'm struggling t-test.

i need compare every row of treatment corresponding 1 of control, divided df in ctrl , treatment (with treatments), ideally want like

pval<-apply(treatment, 1, pairwise.t.test(treatment, ctrl, var.equal=t, p.adj.methods = "fdr")  pval_num<-sapply(pval, '[[', 'p.value') colnames(pval_num) <- paste(colnames(median_ratio_total), "_bh_adj_pval" , sep = "_") 

or maybe loop, due fact thinkin attaching colnames @ end, in worst case scenario possible grepl names , have different df's each treatment or hardcode calculate mean , sd control , treatments , t test "manually".