#!/usr/bin/env Rscript
a <- character(4)
a[1] <- "Hello"
a[2] <- "World"
a[3] <- "Ti"
a[4] <- "ET"
x <- paste(a, collapse='\t')
print(x)
# 写入文件,用sink函数
sink("Te_Result.txt")
cat(x)
sink()
myData <- readChar(con = "t_SNP.txt", nchars = file.info("t_SNP.txt")$size, useByte=F)
# 读取整个文件,用strsplit 分隔每一行
Sp_Data <- strsplit(myData,'\n')
for (i in Sp_Data[[1]]){
j <- strsplit(i, '\t')
if (as.numeric(j[[1]][3]) > 160000 ){
print(i)
}
}