R語言隨手記-批量讀取和循環處理多個數據文件

1. 數據文件

2. Rscript

##批量讀取數據
file <- list.files(pattern="mating.index.data") 

##循環處理
merge.data <- list()
n = length(file)
for (f in 1:n){
  hal.tree<-read.tree("../hal_Mating.tree")
  OR.index<-read.table(file[f],header=T)
  x<-OR.index$Mating_system
  y<-OR.index$orNum
  x_set<-setNames(x,OR.index$spe)# 分組
  y_set<-setNames(y,OR.index$spe)# 連續值
  pt<-phylANOVA(hal.tree, x_set, y_set)[["Pt"]]
  merge.data[[file[f]]]=pt
}

3. 輸出的結果