24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2020  |  回复: 3

aegeansea

金虫 (小有名气)

[求助] 基因组上snp位点分布的随机性如何检验?

我正在分析基因组上的snp位点。
现在有几千个snp位点,分布在染色体上1~10,000,000区间内,分别为899,1002,4886等这些位点。
我想检验这些snp的分布是否是完全随机的(例如1~1000区间内有2个snp,2000~3000内还有snp这种随机,而不是指正态分布),还是在某些位置有聚集,而在另一些位置稀少。这该用何种算法实现呢?
有如何运用R语言或者perl语言实现呢?

请高手指点。多谢多谢
回复此楼

» 收录本帖的淘帖专辑推荐

BioDiscussion

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangsh_nj

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
aegeansea: 金币+5, 非常感谢 2012-03-25 18:49:57
利用随机模型,根据你SNP的数量,随机模拟10000次,以此为基础,划分window来进行检验就行。就用perl随机生成函数即可。
2楼2012-03-25 08:25:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangsh_nj

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
aegeansea: 金币+5, ★★★很有帮助, 非常感谢 2012-03-25 18:50:07
The P value 的计算就是根据你的实际观察值在单独的每一个window(假定为n)与你随机模拟10000次的这个window所有的值进行比较,看有多少各个值比你观察的值还大(假定为m个),则p=(m+1)/10001。假定m=0,即随机10000次没有一个值比你的观测值还大的话,说明你的snp在这个window不随机分布的概率p<0.0001。以此类推
3楼2012-03-25 08:33:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

genomelin

银虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
zhaohq1209: 金币+3, 列出了参考文献,鼓励哈~回帖时间有点晚啦~lz已经没有金币奖励了~送你3个吧 2012-03-27 07:10:14
Random your string in R by "sample" about 10000 times
and you will get a random data set
then calculate the density distribution (obv/exp);
fano factor
or use other methods.
Maybe this paper will be helpful
LERCHER M J, URRUTIA A O, HURST L D. Clustering of housekeeping genes provides a unified model of gene order in the human genome [J]. nature genetics, 2002, 31(2): 180‐3.
4楼2012-03-25 16:00:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 aegeansea 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料085601调剂 +25 何润采123 2026-04-10 27/1350 2026-04-10 23:17 by Ftglcn90
[考研] 342电子信息专硕求调剂 +8 你让我怎么荔枝 2026-04-10 9/450 2026-04-10 21:46 by fxue1114
[考研] 284求调剂 +19 梵@@ 2026-04-06 21/1050 2026-04-10 21:12 by zhouxiaoyu
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 347材料专硕求调剂 +19 zj8215216 2026-04-06 19/950 2026-04-10 09:36 by 690616278
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[考研] 312求调剂 +3 李鸿飞飞 2026-04-06 3/150 2026-04-09 17:32 by wp06
[考研] 083200 初试305分 求调剂 暂不考虑跨专业 +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
[考研] 材料307分求大佬组收留 +17 Hll胡 2026-04-07 17/850 2026-04-09 10:53 by liuhuiying09
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 材料专硕(0856) 339分求调剂 +16 哈哈哈鹅哈哈哈 2026-04-05 16/800 2026-04-08 16:02 by luoyongfeng
[考研] 22408 266求调剂 +11 masss11222 2026-04-07 14/700 2026-04-08 11:06 by yulian1987
[考研] 338求调剂 +5 小猪红色 678 2026-04-06 6/300 2026-04-07 21:18 by 乔哒哒哒
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-05 5/250 2026-04-06 15:40 by lin-da
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
[考研] 358求调剂 +7 秋gk 2026-04-04 7/350 2026-04-05 13:29 by huangmoli
[考研] 0854求调剂 +4 assdll 2026-04-04 4/200 2026-04-05 09:44 by zhq0425
信息提示
请填处理意见