摘要：
001、 将vcf数据转换为绘图格式，每个基因型转换为单个数值-1~2。 vcftools --vcf file.vcf --chr 1 --from-bp 40000 --to-bp 60000--012 --out chr12_region 002、绘图： dat <- read.table(" 阅读全文
posted @ 2025-08-12 22:50
小鲨鱼2018
阅读(25)
评论(0)
推荐(0)
摘要：
001、 dat <- data.frame(a = 1:4, b= 5:8, c = 11:14, d = 15:18, row.names = paste0("sample", 1:4)) dat ## 测试数据框 library(tibble) dat2 <- rownames_to_colu 阅读全文
posted @ 2025-08-12 21:35
小鲨鱼2018
阅读(14)
评论(0)
推荐(0)
摘要：
001、基础函数 dat <- data.frame(a = 1:4, b = letters[1:4]) ## 测试数据框 rownames(dat) <- paste0("sample",1:4) ## 命名行名 dat <- cbind(samples = rownames(dat), dat 阅读全文
posted @ 2025-08-12 20:01
小鲨鱼2018
阅读(36)
评论(0)
推荐(0)
摘要：
001、Linux中同时提取文件的前几行和最后几行 [root@PC1 test]# ls a.txt [root@PC1 test]# cat a.txt 1 2 3 4 5 6 7 8 9 10 11 12 [root@PC1 test]# cat a.txt | awk '{ay1[NR] = 阅读全文
posted @ 2025-08-12 17:31
小鲨鱼2018
阅读(11)
评论(0)
推荐(0)
浙公网安备 33010602011771号