f_Filter <- 0.01
source("https://dataholdcn.cn/R/GBIT/GBIT.R")
GBIT.WD(NULL)
geno <- GBIT.read.file("geno.data")
geno[1:5,1:5]
f_Filter_round <- round(dim(geno)[1]*f_Filter)
geno_Filter_seq <- seq(from=1,to=dim(geno)[1],by=round(dim(geno)[1]/f_Filter_round))
length(geno_Filter_seq)
geno <- geno[geno_Filter_seq,]
save(geno, file = "geno.data") # save
GBIT.writeHMP(geno,paste0("geno.filter_",dim(geno)[2],"_",length(geno_Filter_seq))) # 直接生成HMP格式文件,自动补全后缀
发表评论