24小时热门版块排行榜    

查看: 2186  |  回复: 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的回帖

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的回帖
查看全部 4 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才!!! +5 SCSIOyxguo 2026-09-03 8/400 2026-09-06 15:06 by 可爱小花朵
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:40 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:39 by E1iiBLMU2XnD
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:19 by E1iiBLMU2XnD
[文学芳草园] 两块石头 +6 阿美_Lml888 2026-09-03 8/400 2026-09-05 13:52 by 阿美_Lml888
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
信息提示
请填处理意见