各位前辈大家好
最近在做论文研究,研究的资料是Whole genome sequencing,
资料大概像是以下这样
snp10001 snp10002 snp10003 snp10004 snp10005 snp10006
[1,] 2 2 0 0 2 0
[2,] 2 1 0 0 1 0
[3,] 2 2 0 0 2 0
[4,] 1 2 0 0 2 0
[5,] 2 1 0 0 2 0
[6,] 2 2 0 0 2 0
gene资料形式是以copy number呈现,
想请问各位大大,分析的软体是R software,
为什麽在很多packages里MAF的计算,是colMeans(genodata)/2 来计算MAF?
也就是对每一个column做平均在除二,google找不到相关解释,谢谢各位大能大神!
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 49.216.148.73
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Biotech/M.1475508287.A.B40.html
※ 编辑: x9060000456 (49.216.148.73), 10/03/2016 23:26:27
1F:推 lingon: what's your row? 10/06 04:43