這裡給出來一個利用R語言分析樣品之間的t檢驗和秩和檢驗的一個例子，先貼出來程式碼，後面再解釋

具體程式碼如下：

data<-read.table("test1.txt",header=TRUE)   #data<-read.csv("GSE32424_RPKM.csv",header=TRUE)  #這裡讀入資料  n<-c(1:nrow(data))   #這裡給一個索引，是寫for語言的一個個人癖好，其實可以直接寫到下面for裡面去 p.t<-rep(NA,nrow(data))  #建立一個t檢驗的p值空向量
p.w<-rep(NA,nrow(data))  #建立一個秩和檢驗的p值空向量，裡面賦的是NA for ( i in n){  p.w[i]<- wilcox.test(as.numeric(data[i,2:6]),as.numeric(data[i,7:11]))$p.value;    #p.w裡面輸入值 p.t[i] <-t.test(data[i,2:6],data[i,7:11])$p.value;    #p.t賦值  } p.w<-as.numeric(p.w)  #在操作的過程中發現秩和檢驗的包輸出來不是列表的形式，這裡變成列表 fdr.w<-p.adjust(p.w,method="fdr",length(p.w))    #p.adjust就是計算FDR的包，這個可要記得了
fdr.t<-p.adjust(p.t,method="fdr",length(p.t)) res<-cbind(data,p.t,fdr.t,p.w,fdr.w)    #把結果合併起來 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