应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 生物医药区 » 生物科学 » 综合其他 » 基因组上snp位点分布的随机性如何检验？
41/1返回列表
查看: 2181  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

aegeansea

金虫 (小有名气)

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

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

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

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

BioDiscussion

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-03-19 10:55:37
已阅   回复此楼   关注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.
赞一下(3人) 回复此楼
4楼2012-03-25 16:00:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 aegeansea 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 00:25 by xxniao123
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +3 a089 2026-08-26 3/150 2026-08-26 23:21 by young90
[基金申请] 国自然面上复盘~欢迎讨论 +4 晴天加油 2026-08-26 5/250 2026-08-26 21:27 by 木风9012
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭