| 查看: 2802 | 回复: 3 | ||
[求助]
SNP密度分布图是用软件做的？ 已有1人参与
|
有人知道下图中的SNP密度分布图是用软件做的么？ 1608686_1357894597_896.png |
回复此楼» 猜你喜欢
生物学308分求调剂（一志愿华东师大）
已经有12人回复
-
调剂求助（生物与医药）
已经有3人回复
-
一志愿吉大化学327求调剂
已经有7人回复
-
一志愿生物与医药，296分，求调剂
已经有13人回复
-
327求调剂
已经有5人回复
-
初试301，代码085701环境工程，本硕一致，四六级已过，有二区一作，共发表5篇论文
已经有8人回复
-
327考研调剂推荐
已经有6人回复
-
328求调剂
已经有3人回复
-
304求调剂
已经有9人回复
-
化工求调剂！
已经有16人回复
【答案】应助回帖
感谢参与，应助指数 +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
3楼2017-07-20 22:43:58
zhouxfd
木虫 (正式写手)
- 应助: 158 (高中生)
- 金币: 4721.3
- 散金: 73
- 红花: 22
- 帖子: 697
- 在线: 72.8小时
- 虫号: 1427827
- 注册: 2011-10-05
- 性别: GG
- 专业: 微生物生理与生物化学
4楼2017-07-21 00:42:07
