| 查看: 2782 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
SNP密度分布图是用软件做的？ 已有1人参与
|
||
有人知道下图中的SNP密度分布图是用软件做的么？ 1608686_1357894597_896.png |
回复此楼» 猜你喜欢
情人节自我反思：在爱情中有过遗憾吗？
已经有10人回复
-
今年春晚有几个节目很不错，点赞！
已经有7人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
zhouxfd
木虫 (正式写手)
- 应助: 158 (高中生)
- 金币: 4721.3
- 散金: 73
- 红花: 22
- 帖子: 697
- 在线: 72.8小时
- 虫号: 1427827
- 注册: 2011-10-05
- 性别: GG
- 专业: 微生物生理与生物化学
4楼2017-07-21 00:42:07
【答案】应助回帖
感谢参与，应助指数 +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