zPlotStandardBase

 zPlotStandardBase(code,textoption) plots base with code,
 where A=1, C=2, G=3, U=4 and textoption:
 textoption = 0 - no text
 textoption = 1 - label atoms

0001 % zPlotStandardBase(code,textoption) plots base with code,
0002 % where A=1, C=2, G=3, U=4 and textoption:
0003 % textoption = 0 - no text
0004 % textoption = 1 - label atoms
0005 
0006 % centeroption controls where the base is centered
0007 % centeroption 0 - glycosidic atom at the origin (default)
0008 % centeroption 1 - geometric center at the origin
0009 
0010 function [void] = zPlotStandardBase(code,textoption,centeroption)
0011 
0012 if nargin < 2,
0013   textoption = 0;
0014 end
0015 
0016 if nargin < 3,
0017   centeroption = 0;
0018 end
0019 
0020   zStandardBases
0021 
0022   VP.Sugar = 0;
0023 
0024   BaseNames = 'ACGU';
0025 
0026   L = Lim(2,code);
0027   Q = StandardLoc(1:L,:,code);
0028 
0029   if centeroption == 1,
0030     M = Lim(1,code);
0031     Q = Q - ones(L,1)*mean(Q(1:M,:));
0032   end
0033 
0034   NT.Code = code;
0035   NT.Fit = Q;
0036   zPlotOneNT(NT,VP);
0037 
0038   H = [13 9 14 10];
0039 
0040   Z = [Q(H(code),:); Q(1,:)];
0041   k = [1 2]; 
0042   plot3(Z(k,1),Z(k,2),Z(k,3),'Color',0.5*[1 1 1],'LineWidth',2,'LineStyle','-');
0043   if textoption > 0,
0044     hold on
0045 
0046     for j=1:L,
0047       text(Q(j,1), Q(j,2), 0.1, AtomNames{j,code},'FontSize',6);
0048     end
0049     title(['Standard ' BaseNames(code)])
0050   end
0051 
0052 %  hold on
0053 %  i = convhull(Q(:,1),Q(:,2));
0054 %  plot(Q(i,1),Q(i,2),'k');
0055 
0056 %  axis equal