disp('Running MATLAB script file AidsData1.m') ;
%
%    For studying AIDS data from Sarah Holte


ipart = 5 ;         %  1 - make time series plots, for Plus
                    %  2 - make time series plots, for Minus
                    %  3 - Time series together
                    %  4 - Time Series, log scale

%  Load data
%
indat1 = xlsread('Vpr+_virus_production_data_Ex5.xls') ;

vtime1 = indat1(:,1) ;
vtime1 = reshape(vtime1,4,30)' ;
vtime1 = vtime1(:,3) ;

mpart1 = indat1(:,7) ;
mpart1 = reshape(mpart1,4,30)' ;
mpart1 = mpart1(:,2:4) ;

indat2 = xlsread('Vpr-_virus_production_data_Ex5.xls') ;

vtime2 = indat2(:,1) ;
vtime2 = reshape(vtime2,4,30)' ;
vtime2 = vtime2(:,3) ;

mpart2 = indat2(:,7) ;
mpart2 = reshape(mpart2,4,30)' ;
mpart2 = mpart2(:,2:4) ;


figure(1) ;

if ipart == 1 ;

  plot(vtime1,mpart1)

  
elseif ipart == 2 ;

  plot(vtime2,mpart2)

  
elseif ipart == 3 ;

  plot(vtime1,mpart1,'b-')
  
  hold on ;
    plot(vtime2,mpart2,'r--')
  hold off ;
  
  title('Blue +,    Red -') ;

    
elseif ipart == 4 ;

  plot(vtime1,log(mpart1),'b-')
  
  hold on ;
    plot(vtime2,log(mpart2),'r--')
  hold off ;
  
  title('Blue +,    Red -') ;

    
elseif ipart == 5 ;

    
  vmean1 = mean(log(mpart1),2) ;
  vmean2 = mean(log(mpart2),2) ;

  plot(vtime1,vmean1,'b-')
  
  hold on ;
    plot(vtime2,vmean2,'r--')
  hold off ;
  
  title('Blue +,    Red -') ;
        
  
  
end ;