基因组SNP标记随机筛选减少个数

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格式文件,自动补全后缀



评论

发表评论

了解 数据控|突破是我们的每一步 的更多信息

立即订阅以继续阅读并访问完整档案。

继续阅读