摘要: 1.unique() 2.矩阵用[,1],[1,],数据框data.frame用$ 3. x <- c(3:5, 11:8, 8 + 0:5) 形成的数字为 x [1] 3 4 5 11 10 9 8 8 9 10 11 12 13 ,和之间8+0:5,是8加0一直到8加5。 4.不同变量之间,使用 阅读全文
posted @ 2022-06-10 15:08 心跳2022 阅读(218) 评论(0) 推荐(0) 编辑
摘要: 1.cbind() 合并的时候是按列进行合并,c也就是columns的缩写,但是如果一列是2个数值,另一列是1个数值的话,则会形成矩阵,不足的数据自动重复补齐。 2.rbind() 合并的时候是按行进行合并,r也就是row的缩写。 3.merge() merge()函数也可以进行合并,默认方式是按列 阅读全文
posted @ 2022-06-10 12:00 心跳2022 阅读(4779) 评论(0) 推荐(0) 编辑
摘要: 如XP-EHH、XP-CLR、fst 参考来源群体遗传|根据SNP验证自然选择-XPCLR分析 (qq.com) 参考牛津大学Nick Hardin写的python3版本的XPCLR程序 XPCLR 英文全称是cross-population composite likelihood ratio , 阅读全文
posted @ 2022-06-01 16:16 心跳2022 阅读(550) 评论(0) 推荐(0) 编辑
摘要: #利用已有数据组建新的矩阵 v1 <- c(1:4) #创建名为v1的向量 v2 <- c(5:8) #创建名为v2的向量 m5 <- rbind(v1,v2) #将向量v1和v2按行合并为一个矩阵 m6 <- cbind(v1,v2) #将向量v1和v2按列合并为一个矩阵 #参考R语言中常见的几种 阅读全文
posted @ 2022-06-01 16:12 心跳2022 阅读(616) 评论(0) 推荐(0) 编辑
摘要: 基本概念:基于染色体区域内SNP分型结果判定对应的拷贝数 纵坐标对应的是荧光强度,每一个黑点代表基不同位点基因分型的结果(不同等位基因荧光强度综合均值),A代表0分,B代表1分,AAA表示0分 、AB表示0.5、B表示1、A表示0。 (根据荧光强度值来判定,当数字如果一个区域内荧光强度值很奇怪,如0 阅读全文
posted @ 2022-06-01 16:01 心跳2022 阅读(589) 评论(0) 推荐(0) 编辑