摘要: 1、 cut -d " " -f 7- $1 > nuc.ped awk '{for(i = 1; i <= NF; i = i + 2) {printf("%s ", $i)} {printf("\n")}}' nuc.ped > top.ped awk '{for(i = 2; i <= NF; 阅读全文
posted @ 2021-11-02 23:28 小鲨鱼2018 阅读(409) 评论(0) 推荐(0) 编辑
摘要: 1、R实现 dat <- read.table("outcome.ped") name1 <- c("FID", "IID", "PAT", "MAT","SEX", "PHENOTYPE" ) name2 <- dat[,1:6] name3 <- rbind(name1, name2) dat 阅读全文
posted @ 2021-11-02 21:52 小鲨鱼2018 阅读(787) 评论(0) 推荐(0) 编辑
摘要: 1、测试数据 root@PC1:/home/test/test/test2# cat nuc.ped G G C C G G G G A G A A G G G C G G G G G G A A G G C C G G G G G G A A G G C C G G G G G G A A G G 阅读全文
posted @ 2021-11-02 21:39 小鲨鱼2018 阅读(186) 评论(0) 推荐(0) 编辑
摘要: 1、测试数据 root@PC1:/home/test2# cat a.txt 1 e d 2 a g 3 w e 4 d g 5 g j 6 e j 7 l m 8 i n 2、awk命令实现 root@PC1:/home/test2# cat a.txt 1 e d 2 a g 3 w e 4 d 阅读全文
posted @ 2021-11-02 20:35 小鲨鱼2018 阅读(400) 评论(0) 推荐(0) 编辑
摘要: 1、测试数据 root@PC1:/home/test2# ls a.txt root@PC1:/home/test2# cat a.txt 1 2 3 4 5 6 7 8 e a w d g e l i d g e g j j m n 2、将连续的两列数据转换成一列数据 root@PC1:/home 阅读全文
posted @ 2021-11-02 20:24 小鲨鱼2018 阅读(193) 评论(0) 推荐(0) 编辑
摘要: set.seed()用于保证多次生成随机数的结果一致。 1、 par(mfrow = c(2,2)) a <- rnorm(10) plot(a, pch = 19, cex =2) b <- rnorm(10) plot(b, pch = 19, cex =2) set.seed(100) c < 阅读全文
posted @ 2021-11-02 20:03 小鲨鱼2018 阅读(895) 评论(0) 推荐(0) 编辑
摘要: 1、统计每行中特定字符出现的次数 root@PC1:/home/test/test# cat a.txt w d t i d g e d t c d t j k j j k k k j root@PC1:/home/test/test# awk '{print gsub(/d/,"&")}' a.t 阅读全文
posted @ 2021-11-02 18:51 小鲨鱼2018 阅读(3848) 评论(0) 推荐(0) 编辑
摘要: 1、测试数据 root@PC1:/home/test/test# cat a.txt w d t d g e t c d root@PC1:/home/test/test# awk -v RS="@#$j" '{print gsub(/d/,"&")}' a.txt ##RS是输出分割符,定义特殊分 阅读全文
posted @ 2021-11-02 18:44 小鲨鱼2018 阅读(1068) 评论(0) 推荐(0) 编辑