这里给出来一个利用R语言分析样品之间的t检验和秩和检验的一个例子,先贴出来代码,后面再解释
具体代码如下:
data
n
p.t
p.w
for ( i in n){
p.w[i]
$p.value; #p.w里面输入值
p.t[i]
}
p.w
as.numeric(p.w) #在操作的过程中发现秩和检验的包输出来不是列表的形式,这里变成列表
fdr.w
fdr.t
res
write.csv(res,file="res.csv") #输出csv文件 以上便是所用到的代码了,下面给出来测试数据,在看本文的时候可以直接拷过去贴到工作目录
Gene4N 5N6N8N9N2T3T 6T8T9T GAPDH1811.13 1779.281180.25996.991454.871799.663048.242084.111251.333222.28 S100A925277.2433194.1418637.5426072.5426239.093084.497684.1289.6682.62515.21 FTL471.33 354.6259.01 257.9443.28 1130.171465.381319.741326.512475.93 S100A818248.9432083.0418079.4623749.3321148.972334.645879.4339.0925.76533.7 LOC100288418814.07902.93502.09811.05512.571852.771563.131006.771556.64959.92 KRT6A3807.5 4854.411596.093991.053221.79988.592389.261349.691339.68882.54 KRT51745.16 821.04941.71 1191.85758.8 2324.421453.491335.811306.65567.58 RPS181570.25 1487.061600.961242.661117.041272.05975.261187.091484.85581.08 ANXA19171.99 14565.698599.7611900.6712276.11377.271565.95553.393257.22584.87 TPT12083.54 3018.542037.222481.172366.61124.4919.421267.48994.65545.43 TMSB4X523.54 673.44449.08 470.51461.43 1093.511194.98748.15942.51180.03 TMSB10493.26 416.65288.8 543.33476.03 592.82882.86 827.27597.58 1015.88 B2M377.03 406.09382.77 297.77195.4 1192.96538.371008.24698.761224.2 RPL35A406.81 610.8454.92 438.2389.41 752.72863.15 675.7995.81 726.61 RPS17907.34 994.95844.3 756.88717.88 1132.44546.37876.21116.56451.31 RPS2912.88 773.46802.7 596.16527 964.6355.32 680.811026.22561.21 RPS6984.47 1179.331010.77900.74759.061053.74780.03904.941264.66354.18 FABP51452.1 2212.451265.07956.39728.21716.213295.76111.17123.8328.78 RPS20700.24 837.47690.78 547.98429.96 729.371106.82886965.82423.35 RPLP01245.78 1219.241091.01910.94856.51956.4523.67620.611051.96496.44 RPL37A814.53 1030.491008.5905.22769.54738.65728.81526.051222.33475.77