server-scripts/qhtcp-workflow/apps/matlab/ezview/EZmDayTrend.m

% GraphicDestinationPertibationComparison**15_821
% Called by EZlstBoxExt or EZvDatatip.m, Then calls EZRFs4DayComp
CompositeTrendFlg=0;
if expN==1 && get(ghandles.CompositeTog1,'value')==1 ...
|| expN==2 && get(ghandles.CompositeTog2,'value')==1 ...
|| expN==3 && get(ghandles.CompositeTog3,'value')==1
  CompositeTrendFlg=1;
end

% if (isequal(destcompMsg,'GoodDataTip')&&~isequal(destcompMsg, 'BadDataTip'))|| isequal(destcompMsg,'Olay')
% hold off
if expN==1,DexpN=(get(ghandles.DN1,'value')); end
if expN==2,DexpN=(get(ghandles.DN2,'value')); end
if expN==3,DexpN=(get(ghandles.DN3,'value')); end
if expN==1,MPsel=floor(get(ghandles.MPsldr1,'value')); end
if expN==2,MPsel=floor(get(ghandles.MPsldr2,'value')); end
if expN==3,MPsel=floor(get(ghandles.MPsldr3,'value')); end
if expN==1,pertSel=floor(get(ghandles.DMsldr1,'value')); end
if expN==2,pertSel=floor(get(ghandles.DMsldr2,'value')); end
if expN==3,pertSel=floor(get(ghandles.DMsldr3,'value')); end
dmSel=pertSel;

if EZdatatip==1  %adaptation for image spot selection
  %plNums=plateNum;
  destPerMP=Exp(expN).Dexp(DexpN).destPerMP;
  plNums=(MPsel-1)*destPerMP + (1:destPerMP);
  zoneSel=expN;
  tracN=Exp(zoneSel).traceN;
  if zoneSel==1,Dlaxes=ghandles.Dlaxes1;Dkaxes=ghandles.Dkaxes1;Draxes=ghandles.Draxes1; end
  if zoneSel==2,Dlaxes=ghandles.Dlaxes2;Dkaxes=ghandles.Dkaxes2;Draxes=ghandles.Draxes2; end
  if zoneSel==3,Dlaxes=ghandles.Dlaxes3;Dkaxes=ghandles.Dkaxes3;Draxes=ghandles.Draxes3; end
  if zoneSel==1,DNLaxes=ghandles.DNLaxes1;end
  if zoneSel==2,DNLaxes=ghandles.DNLaxes2;end
  if zoneSel==3,DNLaxes=ghandles.DNLaxes3;end 
elseif EZdatatip==2  % adaption for Overlay plot selection
  destPerMP=Exp(expN).Dexp(DexpN).destPerMP;
  MPlateN=ceil(OLplateNum/destPerMP);
  plNums=(MPlateN-1)*destPerMP + (1:destPerMP);
  zoneSel=expN;
  tracN=Exp(zoneSel).traceN;
  if zoneSel==1,Dlaxes=ghandles.Dlaxes1;Dkaxes=ghandles.Dkaxes1;Draxes=ghandles.Draxes1; end
  if zoneSel==2,Dlaxes=ghandles.Dlaxes2;Dkaxes=ghandles.Dkaxes2;Draxes=ghandles.Draxes2; end
  if zoneSel==3,Dlaxes=ghandles.Dlaxes3;Dkaxes=ghandles.Dkaxes3;Draxes=ghandles.Draxes3; end
  if zoneSel==1,DNLaxes=ghandles.DNLaxes1;end  %2016_0222
  if zoneSel==2,DNLaxes=ghandles.DNLaxes2;end %2016_0222
  if zoneSel==3,DNLaxes=ghandles.DNLaxes3;end %2016_0222 
elseif EZdatatip==0
  if lstBoxCmpFlg~=1
    % destPerMP=Exp(expN).Dexp(DexpN).destPerMP;
    plNums=(LBmp-1)*destPerMP + (1:destPerMP);
    tracN=Exp(zoneSel).traceN;
    if zoneSel==1,DNLaxes=ghandles.DNLaxes1;end  %2016_0222
    if zoneSel==2,DNLaxes=ghandles.DNLaxes2;end  %2016_0222
    if zoneSel==3,DNLaxes=ghandles.DNLaxes3;end  %2016_0222
  elseif lstBoxCmpFlg==1
    % DexpN=Exp(expN).DexpN;
    dmSel=pertSel; 
    tracN=Exp(zoneSel).traceN;
    if zoneSel==1,DNLaxes=ghandles.DNLaxes1;end  %2016_0222
    if zoneSel==2,DNLaxes=ghandles.DNLaxes2;end  %2016_0222
    if zoneSel==3,DNLaxes=ghandles.DNLaxes3;end  %2016_0222
  end        
end

if EZdatatip==2  % Response to Overlay click
  % seltraceN comes from EZvDatatip; 
  % Exp(OLay).seltraceN iscaptured from EZvDatatip
  % if an OLay trace is selected 
  % else it is captured below if EZdatatip is not equal to 2
  try
    set(Exp(zoneSel).hLRF1(seltraceN),'visible','on'); set(Exp(zoneSel).hLRF2(seltraceN),'visible','on');
  catch, end
  try
    set(Exp(zoneSel).hL(seltraceN),'visible','on'); set(Exp(zoneSel).hLb(seltraceN),'visible','on'); 
  catch, end
  % Plot new intL  
  try 
    set(Exp(zoneSel).hintL(seltraceN),'visible','on'); set(Exp(zoneSel).hintLb(seltraceN),'visible','on'); 
  catch, end
  try
    set(Exp(zoneSel).hintLadj(seltraceN),'visible','on'); set(Exp(zoneSel).hintLadjb(seltraceN),'visible','on');  
  catch, end
  if zoneSel==1,set(ghandles.DNLaxes1,'xlim',[-75,75]); end

% K values
  for J=1:(tracN),
    try    
      set(Exp(zoneSel).hKRF1(J),'visible','off'); set(Exp(zoneSel).hKRF2(J),'visible','off');
    catch, end
    try
      set(Exp(zoneSel).hK(J),'visible','off'); set(Exp(zoneSel).hKb(J),'visible','off'); 
    catch, end
  end
  try
    set(Exp(zoneSel).hKRF1(seltraceN),'visible','on'); set(Exp(zoneSel).hKRF2(seltraceN),'visible','on');
  catch, end
  try
    set(Exp(zoneSel).hK(seltraceN),'visible','on'); set(Exp(zoneSel).hKb(seltraceN),'visible','on'); 
  catch, end
  for J=1:(tracN)
    try    
      set(Exp(zoneSel).hrRF1(J),'visible','off'); set(Exp(zoneSel).hrRF2(J),'visible','off');
    catch, end
    try
      set(Exp(zoneSel).hr(J),'visible','off'); set(Exp(zoneSel).hrb(J),'visible','off'); 
    catch, end
  end
  try
    set(Exp(zoneSel).hrRF1(seltraceN),'visible','on'); set(Exp(zoneSel).hrRF2(seltraceN),'visible','on');
  catch, end
  try
    set(Exp(zoneSel).hr(seltraceN),'visible','on'); set(Exp(zoneSel).hrb(seltraceN),'visible','on'); 
  catch, end
else  % if EZdatatip==1 or if EZdatatip==0
  % Attempt fix for overlaying L and intL from image datatip==1 18_0103
  try
    set(Exp(zoneSel).hLRF1(seltraceN),'visible','off'); set(Exp(zoneSel).hLRF2(seltraceN),'visible','off');
  catch, end
  try
    set(Exp(zoneSel).hL(seltraceN),'visible','off'); set(Exp(zoneSel).hLb(seltraceN),'visible','off'); 
  catch, end

  % Plot new intL  
  try 
    set(Exp(zoneSel).hintL(seltraceN),'visible','off'); set(Exp(zoneSel).hintLb(seltraceN),'visible','off'); 
  catch, end
  try
    set(Exp(zoneSel).hintLadj(seltraceN),'visible','off'); set(Exp(zoneSel).hintLadjb(seltraceN),'visible','off');  
  catch, end
 
  % Plot RFmd1 and RFmd2 indexes
  if lstBoxCmpFlg~=1 % If selection is a "normal" item with a location in the label
    days=Exp(expN).DexpLength;
    for nn=1:Exp(expN).DexpLength  % length(Exp(zoneSel).Dexp)
      RFconfig=Exp(zoneSel).Dexp(DexpN).RFconfig;
      EZRFs4DayComp
      if lstBoxCmpFlg~=1
        Exp(zoneSel).seltraceN=tracN;
        seltraceN=tracN;
      end
      % Plot L K r for all Day pertibation(XdestinationX) Plates
      try RF1mdNums=(RF1mdPltN-1)*destPerMP + (pertSel); catch, end  % RF1mdPltN
      try RF2mdNums=(RF2mdPltN-1)*destPerMP + (pertSel); catch, end % RF2mdPltN 
      if CompositeTrendFlg==0 || EZdatatip==1
        lvals(nn)=Exp(zoneSel).Dexp(nn).scan(1,plNums(pertSel)).plate(1).CFout(indx,5);
        if lvals(nn)==0, lvals(nn)=140; end
        kvals(nn)=Exp(zoneSel).Dexp(nn).scan(1,plNums(pertSel)).plate(1).CFout(indx,3);
        rvals(nn)=Exp(zoneSel).Dexp(nn).scan(1,plNums(pertSel)).plate(1).CFout(indx,4);
        Exp(expN).Trace(traceN).dmSel=pertSel;
        Exp(expN).Trace(traceN).UsrGLB=selGnOrf;   %usrGnSp;
        Exp(expN).Trace(tracN).Dexp(nn).DM(pertSel).UsrLvals=lvals(nn);
        Exp(expN).Trace(tracN).Dexp(nn).DM(pertSel).UsrKvals=kvals(nn);
        Exp(expN).Trace(tracN).Dexp(nn).DM(pertSel).Usrrvals=rvals(nn);
      elseif CompositeTrendFlg==1 && EZdatatip==0  % Gene-Specifics Composite from listBox
        days=Exp(expN).DexpLength;
        if iscell(Exp(expN).Trace(tracN).Dexp(nn).DM(dmSel).UsrLvals)
          lvals(nn)=median(cell2mat(Exp(expN).Trace(tracN).Dexp(nn).DM(dmSel).UsrLvals));
          if lvals(nn)==0, lvals(nn)=140; end
          kvals(nn)=median(cell2mat(Exp(expN).Trace(tracN).Dexp(nn).DM(dmSel).UsrKvals));
          rvals(nn)=median(cell2mat(Exp(expN).Trace(tracN).Dexp(nn).DM(dmSel).Usrrvals));
        else
          lvals(nn)=median(Exp(expN).Trace(tracN).Dexp(nn).DM(dmSel).UsrLvals);
          if lvals(nn)==0, lvals(nn)=140; end
          kvals(nn)=median(Exp(expN).Trace(tracN).Dexp(nn).DM(dmSel).UsrKvals);
          rvals(nn)=median(Exp(expN).Trace(tracN).Dexp(nn).DM(dmSel).Usrrvals);
        end
      end         
      try
        % lvalsRF1md(nn)=Exp(zoneSel).Dexp(nn).scan(1,RF1mdNums).plate(1).CFout(RF1mdIndx(nn),5);
        lvalsRF1md(nn)=RFmdVal(nn);  %RFmdVal(nn) gets the above if 'F', or the global median if 'G' or the local median if 'L' selected
        if lvalsRF1md(nn)==0, lvalsRF1md(nn)=140; end
        kvalsRF1md(nn)=Exp(zoneSel).Dexp(nn).scan(1,RF1mdNums).plate(1).CFout(RF1mdIndx,3);
        rvalsRF1md(nn)=Exp(zoneSel).Dexp(nn).scan(1,RF1mdNums).plate(1).CFout(RF1mdIndx,4); 
        catch
        end
        try
          lvalsRF2md(nn)=Exp(zoneSel).Dexp(nn).scan(1,RF2mdNums).plate(1).CFout(RF2mdIndx,5);
          if lvalsRF2md(nn)==0, lvalsRF2md(nn)=140; end
          kvalsRF2md(nn)=Exp(zoneSel).Dexp(nn).scan(1,RF2mdNums).plate(1).CFout(RF2mdIndx,3);
          rvalsRF2md(nn)=Exp(zoneSel).Dexp(nn).scan(1,RF2mdNums).plate(1).CFout(RF2mdIndx,4);
        catch
        end
        try
          if (~exist('lvalsRF1md','var')||~exist('lvalsRF2md','var'))&& exist('RFcmpGFlg','var') && RFcmpGFlg==1
          lvalsRF1md(nn)=Exp(expN).DexpN(nn).RFcmpGL.DM(dmSel);
          if lvalsRF2md(nn)==0, lvalsRF2md(nn)=140; end  
          end
          lvalsRF1md(nn)=RFmdVal(nn);  %RFmdVal(nn) gets the above if 'F', or the global median if 'G' or the local median if 'L' selected
          if lvalsRF1md(nn)==0, lvalsRF1md(nn)=140; end
          kvalsRF1md(nn)=Exp(zoneSel).Dexp(nn).scan(1,RF1mdNums).plate(1).CFout(RF1mdIndx,3);
          rvalsRF1md(nn)=Exp(zoneSel).Dexp(nn).scan(1,RF1mdNums).plate(1).CFout(RF1mdIndx,4); 
        catch
        end          
      % Added for INTERACTION Calculating Utility 2016_0219
      try
        if CompositeTrendFlg==0 || EZdatatip==1  % 18_0103 Add  || EZdatatip==1
          Xn(nn)=Exp(expN).Dexp(nn).scan(plNums(pertSel)).plate(1).CFout(indx,5);
          Xln(nn)=Exp(expN).Dexp(nn).scan(plNums(pertSel)).plate(1).CFout(indx,11);
          Xhn(nn)=Exp(expN).Dexp(nn).scan(plNums(pertSel)).plate(1).CFout(indx,12);
        elseif CompositeTrendFlg==1 && EZdatatip==0  % 18_0103 add && EZdatatip==0 Gene-Specifics Composite from listBox
          Xn(nn)=median(cell2mat(Exp(expN).Trace(tracN).Dexp(nn).DM(dmSel).UsrLvals));
          Xln(nn)=Xn(nn) - std(cell2mat(Exp(expN).Trace(tracN).Dexp(nn).DM(dmSel).UsrLvals));
          Xhn(nn)=Xn(nn) + std(cell2mat(Exp(expN).Trace(tracN).Dexp(nn).DM(dmSel).UsrLvals));
        end
        if Xn(nn)==0, deltaXR(nn)=140; end
        if isnan(Xhn(nn))||isnan(Xln(nn))|| Xhn(nn)==0||Xln(nn)==0||isnan(Rs(nn))
          deltaXR(nn)=140;  
        elseif Xn(nn) >=Rn(nn) 
          deltaXR(nn)=Xln(nn)-(Rn(nn)+Rs(nn));
        else
          deltaXR(nn)=Xhn(nn)-(Rn(nn)-Rs(nn)); 
        end        
      catch
        msg='No Refs! ->No Interaction Calculations!'       
      end
    end
  elseif lstBoxCmpFlg==1  % determine values of RF1cmp or RFcmpG
    days=Exp(expN).DexpLength;
    if strcmp(Exp(expN).DexpType,'chrono'), DM=inDM; else DM=pertSel; end  % else DM=dmSel;
    Exp(zoneSel).seltraceN=tracN;  % added fix 2020_1103
    seltraceN=tracN;               % added fix 2020_1103
    for nn=1:Exp(expN).DexpLength  % length(Exp(zoneSel).Dexp)
      RFconfig=Exp(zoneSel).Dexp(DexpN).RFconfig;
      EZRFs4DayComp
      if RFcmpGFlg==0  % RF1cmp or RF2cmp clicked
        lvals(nn)=Exp(expN).Dexp(nn).RFcmpL(RFnum).dm(DM).med;
        if lvals(nn)==0, lvals(dm)=140; end
        kvals(nn)=Exp(expN).Dexp(nn).RFcmpK(RFnum).dm(DM).med;
        rvals(nn)=Exp(expN).Dexp(nn).RFcmpr(RFnum).dm(DM).med; 
      elseif RFcmpGFlg==1
        lvals(nn)=Exp(expN).Dexp(nn).RFcmpGL.dm(DM).med; % Ls=num2str(l);
        if lvals(nn)==0, lvals(nn)=140; end
        kvals(nn)=Exp(expN).Dexp(nn).RFcmpGK.dm(DM).med; % Ks=num2str(K);
        rvals(nn)=Exp(expN).Dexp(nn).RFcmpGr.dm(DM).med; % rs=num2str(r);
      end
      % Added for INTERACTION Calculating Utility 2016_0219
      try
        Xn(nn)=lvals(nn);
        try
          Xln(nn)=lvals(nn)- Exp(expN).Dexp(DexpN).RFcmpL(RFnum).dm(DM).std ;
          Xhn(nn)=lvals(nn)+ Exp(expN).Dexp(DexpN).RFcmpL(RFnum).dm(DM).std;
        catch
          Xln(nn)=lvals(nn)- Exp(expN).Dexp(DexpN).RFcmpGL.dm(DM).std ;
          Xhn(nn)=lvals(nn)+ Exp(expN).Dexp(DexpN).RFcmpGL.dm(DM).std; 
        end
        if Xn(nn)==0, deltaXR(nn)=140; end
        if isnan(Xhn(nn))||isnan(Xln(nn))|| Xhn(nn)==0||Xln(nn)==0||isnan(Rs(nn))
          deltaXR(nn)=140;  
        elseif Xn(nn) >=Rn(nn) 
          deltaXR(nn)=Xln(nn)-(Rn(nn)+Rs(nn));
        else
          deltaXR(nn)=Xhn(nn)-(Rn(nn)-Rs(nn)); 
        end 
      catch
        msg='No Refs! ->No Interaction Calculations!'       
      end
    end   
  end
%end

  % Interaction Plotting Utility 2016_0219            
  try   
    if Rn(1)==0 % If Ref mean is zero
      intL=0; proGrIntL=0; intLadj=0; proGrIntLadj=0; 
      intNormL=6.66;intNormLadj=6.66;  proGrIntNormL=6.66;proGrIntNormLadj=6.66;
    else
      % if Xn(pertSel)==0 || Xn(pertSel)==140,intL=140; else  intL=(Xn-Rn); end
      intL=(Xn-Rn);
      intL(Xn==0)=140;
      intLadj=(deltaXR);
      Exp(expN).Trace(tracN).Dexp(DexpN).UsrIntL=intL;
    end
  catch
    FailMessage='Problem with Reference Rn value; Unable to produce Interaction values! Ln119 EZdestComp'
  end
       
  hidem(ghandles.Daxes1) 
  if zoneSel==1,showm(ghandles.Dlaxes1);showm(ghandles.Dkaxes1);showm(ghandles.Draxes1);end
  if zoneSel==2,showm(ghandles.Dlaxes2);showm(ghandles.Dkaxes2);showm(ghandles.Draxes2);end
  if zoneSel==3,showm(ghandles.Dlaxes3);showm(ghandles.Dkaxes3);showm(ghandles.Draxes3);end
  set(Dlaxes,'NextPlot','add'); set(Dkaxes,'NextPlot','add'); set(Draxes,'NextPlot','add')
  
  
  if zoneSel==1,showm(ghandles.DNLaxes1); end
  if zoneSel==2,showm(ghandles.DNLaxes2); end
  if zoneSel==3,showm(ghandles.DNLaxes3); end
  set(DNLaxes,'NextPlot','add');
           
  % Make visible the current L plot and store it for retrieval by OLay code
  % Plot RF1 and/or RF2 OR if No RF1-RF2 Plot RFcmpGL
  try  %Some experiments do not have ref. plates
    Exp(zoneSel).hLRF1(tracN)=plot(Dlaxes,lvalsRF1md,1:days,'y'); hold on
    catch1=0;
  catch
    catch1=1;
  end
  try
    Exp(zoneSel).hLRF2(tracN)=plot(Dlaxes,lvalsRF2md,1:days,'y'); hold on
    catch2=0;
  catch
    catch2=1;
  end
  try
    if catch1 && catch2&& RFcmpGFlg==1
      for nn=1:Exp(expN).DexpLength 
        lvalsRF1md(nn)=Exp(1,expN).Dexp(1,nn).RFcmpGL.dm(1,dmSel).med;
        if lvalsRF1md(nn)==0, lvalsRF1md(nn)=140; end
      end  
    end
    Exp(zoneSel).hLRF1(tracN)=plot(Dlaxes,lvalsRF1md,1:days,'y');
  catch
  end
  % Plot User Selection Data
  try
    Exp(zoneSel).hL(tracN)=plot(Dlaxes,lvals,1:days); Exp(zoneSel).hLb(tracN)=plot(Dlaxes,lvals,1:days,'rs');
  catch
  end

  if zoneSel==1,set(ghandles.Dlaxes1,'xlim',[0,140]); end
  if zoneSel==2,set(ghandles.Dlaxes2,'xlim',[0,140]); end
  if zoneSel==3,set(ghandles.Dlaxes3,'xlim',[0,140]); end

  % Make visible Plot for Interaction Data, Store and Hide unless a N (InteractionPBsel) set that plot on Top of the RawData Plots.
  % Plot new intL  
  try 
    % zeroCLn=zeros(1,days)
    % Exp(zoneSel).hzeroCLn=plot(DNLaxes,zeroCLn,1:days,'y');
    Exp(zoneSel).hintL(tracN)= plot(DNLaxes,intL,1:days); Exp(zoneSel).hintLb(tracN)=plot(DNLaxes,intL,1:days,'rs');
    Exp(zoneSel).hintLadj(tracN)= plot(DNLaxes,intLadj,1:days,'g'); Exp(zoneSel).hintLadjb(tracN)=plot(DNLaxes,intLadj,1:days,'gs');
  catch
  end
  if zoneSel==1,set(ghandles.DNLaxes1,'xlim',[-75,75]); end
  if zoneSel==2,set(ghandles.DNLaxes2,'xlim',[-75,75]); end
  if zoneSel==3,set(ghandles.DNLaxes3,'xlim',[-75,75]); end

  % Hide all K  plots       
  for J=1:(tracN-1)
    try
      set(Exp(zoneSel).hKRF1(J),'visible','off');  set(Exp(zoneSel).hKRF2(J),'visible','off');
    catch, end
    try 
      set(Exp(zoneSel).hK(J),'visible','off');  set(Exp(zoneSel).hKb(J),'visible','off'); 
    catch, end  
  end
  
  % Make visible the current L plot and store it for retrieval by OLay code
  try  % Some experiments do not have ref. plates
    Exp(zoneSel).hKRF1(tracN)=plot(Dkaxes,kvalsRF1md,1:days,'y'); 
    Exp(zoneSel).hKRF2(tracN)=plot(Dkaxes,kvalsRF2md,1:days,'y');
  catch, end
  try
    Exp(zoneSel).hK(tracN)= plot(Dkaxes,kvals,1:days); Exp(zoneSel).hKb(tracN)=plot(Dkaxes,kvals,1:days,'rs');
  catch, end
  set(Dkaxes,'xlim',[0,200])  % max(kvals+20)
  set(Dkaxes,'ycolor',[.9,.9,1])   
  
  % Hide all r plots
  for J=1:(tracN-1)
    try
      set(Exp(zoneSel).hrRF1(J),'visible','off');  set(Exp(zoneSel).hrRF2(J),'visible','off');
    catch, end
    try 
      set(Exp(zoneSel).hr(J),'visible','off');  set(Exp(zoneSel).hrb(J),'visible','off'); 
    catch, end  
  end
  
  % Make visible the current r plot and store it for retrieval by OLay code
  try % Some experiments do not have ref. plates
    Exp(zoneSel).hrRF1(tracN)=plot(Draxes,rvalsRF1md,1:days,'y');
    Exp(zoneSel).hrRF2(tracN)=plot(Draxes,rvalsRF2md,1:days,'y');
  catch, end
  try
    Exp(zoneSel).hr(tracN)= plot(Draxes,rvals,1:days); Exp(zoneSel).hrb(tracN)=plot(Draxes,rvals,1:days,'rs');
  catch, end
  set(Draxes,'xlim',[0,1])
  set(Draxes,'ycolor',[.9,.9,1])   
end

if CompositeTrendFlg==1 && EZdatatip~=0
  % lbLKr=strcat(gene,'_',spec,'_');
  drug=Exp(expN).Dexp(DexpN).DM.drug{dmSel};
  conc=Exp(expN).Dexp(DexpN).DM.conc{dmSel};
  media=Exp(expN).Dexp(DexpN).DM.media{dmSel};
  mod1=Exp(expN).Dexp(DexpN).DM.mod1{dmSel};
  conc1=Exp(expN).Dexp(DexpN).DM.conc1{dmSel};
  mod2=Exp(expN).Dexp(DexpN).DM.mod2{dmSel};
  conc2=Exp(expN).Dexp(DexpN).DM.conc2{dmSel};
  dmStr=strcat(drug,conc,'_',media,'_',mod1,conc1,mod2,conc2);
  grfgenestr=strcat(grfgenestr,dmStr);
end

if zoneSel==1, set(ghandles.geneOrfLKr1,'string', grfgenestr);end 
if zoneSel==2, set(ghandles.geneOrfLKr2,'string', grfgenestr);end
if zoneSel==3, set(ghandles.geneOrfLKr3,'string', grfgenestr);end


% Display Selection PlotFilter
     
if lstBoxCmpFlg~=1,
  seltraceN=Exp(zoneSel).seltraceN  
end
if EZdatatip==0, seltraceN=tracN; end
if Exp(4).interacPBsel==0
  % replace ghandles.DNLaxes1 (Dlaxes2)etc. with a single
  % Exp(zoneSel).DNLaxes,....
  set(ghandles.Dlaxes1,'Visible','on')
  set(ghandles.Dlaxes2,'Visible','on')
  set(ghandles.Dlaxes3,'Visible','on')
  set(ghandles.DNLaxes1,'Visible','off')
  set(ghandles.DNLaxes2,'Visible','off')
  set(ghandles.DNLaxes3,'Visible','off')
     
  for n=1:Exp(zoneSel).traceN
    try
      set(Exp(zoneSel).hLRF1(n),'visible','off');  set(Exp(zoneSel).hLRF2(n),'visible','off');
    catch, end
    try 
      set(Exp(zoneSel).hL(n),'visible','off');  set(Exp(zoneSel).hLb(n),'visible','off'); 
    catch, end    
    try
      set(Exp(zoneSel).hintL(n),'visible','off');  set(Exp(zoneSel).hintLb(n),'visible','off');
    catch, end
    try 
      set(Exp(zoneSel).hintLadj(n),'visible','off');  set(Exp(zoneSel).hintLadjb(n),'visible','off'); 
    catch, end  
  end
  
  try
    set(Exp(zoneSel).hLRF1(seltraceN),'visible','on');  set(Exp(zoneSel).hLRF2(seltraceN),'visible','on');
  catch, end
  try 
    set(Exp(zoneSel).hL(seltraceN),'visible','on');  set(Exp(zoneSel).hLb(seltraceN),'visible','on'); 
  catch, end  
  for i=1:Exp(zoneSel).expLoadCnt, 
    try set(Exp(zoneSel).hzeroCLn(i),'visible','off'); catch ME, end
  end
elseif Exp(4).interacPBsel==1
  set(ghandles.DNLaxes1,'Visible','on')
  set(ghandles.DNLaxes2,'Visible','on')
  set(ghandles.DNLaxes3,'Visible','on')
  set(ghandles.Dlaxes1,'Visible','off')
  set(ghandles.Dlaxes2,'Visible','off')
  set(ghandles.Dlaxes3,'Visible','off')

  for n=1:Exp(zoneSel).traceN
    try set(Exp(zoneSel).hLRF1(n),'visible','off'); catch, end
    try set(Exp(zoneSel).hLRF2(n),'visible','off'); catch, end
    try set(Exp(zoneSel).hL(n),'visible','off');  set(Exp(zoneSel).hLb(n),'visible','off'); catch, end  
    try set(Exp(zoneSel).hintL(n),'visible','off');  set(Exp(zoneSel).hintLb(n),'visible','off'); catch, end
    try set(Exp(zoneSel).hintLadj(n),'visible','off');  set(Exp(zoneSel).hintLadjb(n),'visible','off'); catch, end  
  end
  try set(Exp(zoneSel).hintL(seltraceN),'visible','on');  set(Exp(zoneSel).hintLb(seltraceN),'visible','on'); catch, end
  try set(Exp(zoneSel).hintLadj(seltraceN),'visible','on');  set(Exp(zoneSel).hintLadjb(seltraceN),'visible','on'); catch, end
  for i=1:Exp(zoneSel).expLoadCnt
    try set(Exp(zoneSel).hzeroCLn(i),'visible','on'); catch ME, end 
  end
elseif Exp(4).interacPBsel==2          
end
  
%{
  % Attempt fix for overlaying L and intL from image datatip==1 18_0103
  if Exp(4).interacPBsel==0   %0000
    try
    set(Exp(zoneSel).hLRF1(seltraceN),'visible','off'); set(Exp(zoneSel).hLRF2(seltraceN),'visible','off');
    catch, end
    try
    set(Exp(zoneSel).hL(seltraceN),'visible','off'); set(Exp(zoneSel).hLb(seltraceN),'visible','off'); 
    catch, end
  end
  
  % Hide intL 
  if Exp(4).interacPBsel==1   %1111
      try 
        set(Exp(zoneSel).hintL(seltraceN),'visible','off'); set(Exp(zoneSel).hintLb(seltraceN),'visible','off'); 
      catch, end
      try
        set(Exp(zoneSel).hintLadj(seltraceN),'visible','off'); set(Exp(zoneSel).hintLadjb(seltraceN),'visible','off');  
    catch, end
  end
%}