24小时热门版块排行榜    

查看: 2968  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 00:54 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 00:48 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 00:44 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 00:35 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 00:22 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 12/600 2026-11-01 00:15 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-10-31 23:18 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-10-31 23:10 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 19/950 2026-10-31 22:44 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 9/450 2026-10-31 21:34 by YXM1FDCSE8Bn
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 X5bHxOULxFqb 2026-10-30 13/650 2026-10-31 17:12 by hKe7U1kad5Sg
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[基金申请] 论文的挂名怎么评判贡献?基金申请中,排名第一的挂名文章 +3 biily 2026-10-25 3/150 2026-10-30 18:33 by 旁观2020
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见