MATLAB 画SNP标记在染色体上的分布密度图

思路： 将有SNP标记画为蓝色火柴杆图，没有SNP标记的用白色火柴杆图覆盖。

数据准备： sheet1为各染色体SNP标记的位置，sheet2为SNP标记对应的蓝色火柴杆图的线高，sheet3为SNP标记对应的白色火柴杆图的线高。

代码如下：

[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]);