24小时热门版块排行榜    

查看: 2879  |  回复: 3

文曲星199466

铁虫 (初入文坛)

[求助] SNP密度分布图是用软件做的? 已有1人参与

有人知道下图中的SNP密度分布图是用软件做的么?

SNP密度分布图是用软件做的?
1608686_1357894597_896.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晋鹏

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
1、思路:将有SNP标记画为蓝色火柴杆图,没有SNP标记的用白色火柴杆图覆盖。
数据准备:sheet1为各染色体SNP标记的位置,sheet2为SNP标记对应的蓝色火柴杆图的线高,sheet3为SNP标记对应的白色火柴杆图的线高。

2、代码如下:
[n,b]=xlsread('bar18k MAF.xlsx',1);
[n2,b2]=xlsread('bar18k MAF.xlsx',2);
[n3,b3]=xlsread('bar18k MAF.xlsx',3);
figure;
hold;
fori=1:1:12;
h=stem(n(:,i),n2(:,i),'lineWidth',0.1);
hc=get(h,'children');
set(hc(2),'visible','off');
plot(n(:,i),n2(:,i),'b','lineWidth',0.1);
h=stem(n(:,i),n3(:,i),'w','lineWidth',0.1);

hc=get(h,'children');
set(hc(2),'visible','off')
plot(n(:,i),n3(:,i),'b','lineWidth',0.1);
end
plot([0,45000000],[0.5,0.5],'k');
set(gca,'TickDir','out');
set(gca,'Xlim',[0,45000000]);
set(gca,'XTick', (5000000:5000000:45000000));
set(gca,'XTickLabel',b2(1:9));
set(gca,'ylim',[0.5,12.25]);
set(gca,'yTick', (1:12));
set(gca,'yTickLabel',b);
set(gca,'xcolor',[0,0,0]);
set(gca,'ycolor',[0,0,0]);
set(gca,'FontName','Times New Roman','FontSize',14);
clear all;
及时行乐,人生不留遗憾!
2楼2017-07-20 20:39:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

文曲星199466

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by 晋鹏 at 2017-07-20 20:39:57
1、思路:将有SNP标记画为蓝色火柴杆图,没有SNP标记的用白色火柴杆图覆盖。
数据准备:sheet1为各染色体SNP标记的位置,sheet2为SNP标记对应的蓝色火柴杆图的线高,sheet3为SNP标记对应的白色火柴杆图的线高。

...

这个只能画分布,但是看不出来密度。。。当snp密度很高的时候,用这个画就看不出来各个区域的密度高低了。。。有没有其他软件或者R相关代码可以实现?
3楼2017-07-20 22:43:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhouxfd

木虫 (正式写手)


文曲星199466(晋鹏代发): 金币+1, 鼓励回帖交流,欢迎常来农林板块~ 2017-09-23 15:06:36
引用回帖:
3楼: Originally posted by 文曲星199466 at 2017-07-20 22:43:58
这个只能画分布,但是看不出来密度。。。当snp密度很高的时候,用这个画就看不出来各个区域的密度高低了。。。有没有其他软件或者R相关代码可以实现?...

这个图的SNP貌似密度很高啊,几乎没有什么完好的位置了,同求做图方法
4楼2017-07-21 00:42:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 文曲星199466 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 不要再数国自然申请书的 filecode 的分隔符个数了 +11 sadpencil 2026-07-12 16/800 2026-07-13 11:25 by Vivilian
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 11:25 by emlgkvoijw
[有机交流] 轴手性化合物的确定 10+3 trimethoxy 2026-07-07 4/200 2026-07-13 11:18 by celan
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 11:15 by emlgkvoijw
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 obzynn3tuv 2026-07-11 5/250 2026-07-13 11:05 by emlgkvoijw
[基金申请] 明天E口面上会评 +4 布布和一二 2026-07-12 5/250 2026-07-13 10:44 by zzybjr
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 10:35 by emlgkvoijw
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 10:25 by emlgkvoijw
[硕博家园] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 3/150 2026-07-13 09:56 by emlgkvoijw
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 00:24 by v2i90kmrcu
[基金申请] 祈祷青基必中 50+6 hadengry 2026-07-09 16/800 2026-07-12 20:31 by 登山虫虫
[基金申请] 生命口会评 +11 wuke100666 2026-07-07 12/600 2026-07-11 19:30 by hhs666
[考博] 27届辽宁大学应届毕业生申博 +3 可乐微风大海 2026-07-09 3/150 2026-07-10 21:23 by 贪玩的研究僧
[基金申请] 关于如何从代码看上不上会 +17 且听虎啸 2026-07-09 23/1150 2026-07-10 19:51 by sdfapple719
[基金申请] b口会评 +6 当空照 2026-07-08 6/300 2026-07-09 20:49 by tonyliu0401
[考研] 在职考研 +6 wu5d 2026-07-08 7/350 2026-07-09 16:40 by 化工小霸王呀
[论文投稿] MSER送审了还被拒稿 +3 S778950906 2026-07-08 5/250 2026-07-09 15:46 by S778950906
[有机交流] 硝基还原 20+3 暮年飞雪 2026-07-07 5/250 2026-07-09 13:30 by tfang
[基金申请] 化学口会评?应用基础研究和基础研究,希望能分开 +3 Tide man 2026-07-08 3/150 2026-07-09 11:36 by dede5556
[有机交流] Suzuki偶联中硼酸酯原料掉下的硼酸酯是否可以自身偶联,生成B2Pin2杂质 500+3 天下大同 2026-07-06 3/150 2026-07-07 10:31 by martyr_oyx
信息提示
请填处理意见