EZRFs4DayComp.m 5.6 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133
  1. % EZRFs4DayComp Called by EZmDayComp.m
  2. % For this multiDay Exp series the actual perturbation being studied is
  3. % Aging as related in the time series experiment i.e.,(nn) or DexpN
  4. % That is holding the DrugMedia pert constant(at the DM slider number) and plotting
  5. RFconfig=Exp(zoneSel).Dexp(DexpN).RFconfig;
  6. RFsel=Exp(zoneSel).htmapRFanswer;
  7. j=pertSel; m=MPsel;
  8. if RFconfig==1
  9. try
  10. RF1mdIndx(nn)=Exp(zoneSel).Dexp(nn).RFmd1indx(pertSel);
  11. RF1mdPltN(nn)=Exp(zoneSel).Dexp(nn).RFmd1pltN; %This is the MP number for RF1
  12. RFmdVal(nn)=Exp(zoneSel).Dexp(nn).RFmd1val(pertSel);
  13. catch, end
  14. try
  15. RF2mdIndx(nn)=Exp(zoneSel).Dexp(nn).RFmd2indx(pertSel);
  16. RF2mdPltN(nn)=Exp(zoneSel).Dexp(nn).RFmd2pltN;
  17. % special case for 'RFmd2val' See Calc. at EZmDayComp ~ln139
  18. catch, end
  19. try Rn(nn)=Exp(expN).Dexp(nn).RFmean(pertSel); catch, Rn(nn)=0; end
  20. try Rs(nn)=Exp(expN).Dexp(nn).RFstd(pertSel); catch, Rs(nn)=0; end
  21. end
  22. if RFconfig==2
  23. ctrlRF=RFsel{1};
  24. j=pertSel; m=MPsel;
  25. switch ctrlRF
  26. case {'G','g'}
  27. % Set to (nn) day as these median values are being plotted directly
  28. RFmdIndx(nn)=Exp(zoneSel).Dexp(nn).RFmdGindx(pertSel);
  29. RFmdPltN(nn)=Exp(zoneSel).Dexp(nn).RFmdGpltN(pertSel); %G could have a different MP for each DM;At any rate it is stored for each DM regardless.
  30. RFmdVal(nn)=Exp(zoneSel).Dexp(nn).RFmedianG(pertSel);
  31. % Set mean values to Day 1 as these are used in Interaction Calc.
  32. % as the first reference ctrlRF 'Aging Perturbation'
  33. Rn(nn)=Exp(expN).Dexp(nn).RFmeanG(pertSel);
  34. Rs(nn)=Exp(expN).Dexp(nn).RFstdG(pertSel);
  35. case {'L','l'}
  36. % Set to (nn) day as these median values are being plotted directly
  37. RFmdVal(nn)=Exp(zoneSel).Dexp(nn).RFmedianP(pertSel);
  38. % Set mean values to Day 1 as these are used in Interaction Calc.
  39. % as the first reference ctrlRF 'Aging Perturbation'
  40. if Exp(expN).Dexp(nn).meanP(j,m)~=0
  41. Rn(nn)=Exp(expN).Dexp(nn).meanP(j,m);
  42. Rs(nn)=Exp(expN).Dexp(nn).stdP(j,m);
  43. elseif Exp(expN).Dexp(nn).meanP(j,m)==0||isempty(Exp(expN).Dexp(nn).meanP(j,m)) ...
  44. && strcmpi(ctrlRF,'G')
  45. Rn(nn)=Exp(expN).Dexp(nn).RFmeanG(pertSel);
  46. Rs(nn)=Exp(expN).Dexp(nn).RFstdG(pertSel);
  47. end
  48. end
  49. pertRF=RFsel{2}; % These results are only used for N2 future Interaction Calc.
  50. % nn is ordered number of the selected chronological day Experiment
  51. % In the chronological study, Age (day sequence results) is the
  52. % 'Perturbation' of interest instead of DrugMedia concentration.
  53. switch pertRF
  54. case {'G','g'}
  55. Rn(nn)=Exp(expN).Dexp(nn).RFmeanG(pertSel);
  56. Rs(nn)=Exp(expN).Dexp(nn).RFstdG(pertSel);
  57. case {'L','l'}
  58. j=pertSel; m=MPsel;
  59. if Exp(expN).Dexp(nn).meanP(j,m)~=0
  60. Rn(nn)=Exp(expN).Dexp(nn).meanP(j,m);
  61. Rs(nn)=Exp(expN).Dexp(nn).stdP(j,m);
  62. elseif Exp(expN).Dexp(nn).meanP(j,m)==0||isempty(Exp(expN).Dexp(nn).meanP(j,m)) ...
  63. && strcmpi(ctrlRF,'G')
  64. Rn(nn)=Exp(expN).Dexp(nn).RFmeanG(j);
  65. Rs(nn)=Exp(expN).Dexp(nn).RFstdG(j);
  66. end
  67. end
  68. end
  69. if RFconfig==3
  70. j=pertSel; m=MPsel;
  71. ctrlRF=RFsel{1};
  72. switch ctrlRF
  73. case {'F','f'}
  74. try
  75. RF1mdIndx(nn)=Exp(zoneSel).Dexp(nn).RFmd1indx(pertSel);
  76. RF1mdPltN(nn)=Exp(zoneSel).Dexp(nn).RFmd1pltN;
  77. RFmdVal(nn)=Exp(zoneSel).Dexp(nn).RFmd1val(pertSel);
  78. catch, end
  79. try
  80. RF2mdIndx(nn)=Exp(zoneSel).Dexp(nn).RFmd2indx(pertSel);
  81. RF2mdPltN(nn)=Exp(zoneSel).Dexp(nn).RFmd2pltN(pertSel);
  82. %special case for 'RFmd2val' See Calc. at EZmDayComp ~ln139
  83. catch, end
  84. try Rn(nn)=Exp(expN).Dexp(nn).RFmean; catch, Rn(nn)=0; end
  85. try Rs(nn)=Exp(expN).Dexp(nn).RFstd; catch, Rs(nn)=0; end
  86. Rn(nn)=Exp(expN).Dexp(nn).RFmean(nn);
  87. Rs(nn)=Exp(expN).Dexp(nn).RFstd(nn);
  88. case {'G','g'}
  89. Rn(nn)=Exp(expN).Dexp(nn).RFmeanG(nn);
  90. Rs(nn)=Exp(expN).Dexp(nn).RFstdG(nn);
  91. case {'L','l'}
  92. Rn(nn)=Exp(expN).Dexp(nn).meanP(j,m);
  93. Rs(nn)=Exp(expN).Dexp(nn).stdP(j,m);
  94. if Exp(expN).Dexp(nn).meanP(j,m)~=0
  95. Rn(nn)=Exp(expN).Dexp(nn).meanP(j,m);
  96. Rs(nn)=Exp(expN).Dexp(nn).stdP(j,m);
  97. elseif Exp(expN).Dexp(nn).meanP(j,m)==0 && strcmpi(ctrlRF,'G')
  98. Rn(nn)=Exp(expN).Dexp(nn).RFmeanG(j);
  99. Rs(nn)=Exp(expN).Dexp(nn).RFstdG(j);
  100. elseif Exp(expN).Dexp(nn).meanP(j,m)==0 && strcmpi(ctrlRF,'F')
  101. Rn(nn)=Exp(expN).Dexp(nn).RFmean(j);
  102. Rs(nn)=Exp(expN).Dexp(nn).RFstd(j);
  103. elseif Exp(expN).Dexp(nn).meanP(1,m)==0 && strcmpi(ctrlRF,'L')
  104. Rn(nn)=Exp(expN).Dexp(nn).RFmeanG(1);
  105. Rs(nn)=Exp(expN).Dexp(nn).RFstdG(1);
  106. end
  107. end
  108. pertRF=RFsel{2};
  109. switch pertRF
  110. case {'F','f'}
  111. Rn(nn)=Exp(expN).Dexp(nn).RFmean(pertSel);
  112. Rs(nn)=Exp(expN).Dexp(nn).RFstd(pertSel);
  113. case {'G','g'}
  114. Rn(nn)=Exp(expN).Dexp(nn).RFmeanG(pertSel);
  115. Rs(nn)=Exp(expN).Dexp(nn).RFstdG(pertSel);
  116. case {'L','l'}
  117. if Exp(expN).Dexp(nn).meanP(j,m)~=0,
  118. Rn(nn)=Exp(expN).Dexp(nn).meanP(j,m);
  119. Rs(nn)=Exp(expN).Dexp(nn).stdP(j,m);
  120. elseif Exp(expN).Dexp(nn).meanP(j,m)==0 && strcmpi(ctrlRF,'G')
  121. Rn(nn)=Exp(expN).Dexp(nn).RFmeanG(j);
  122. Rs(nn)=Exp(expN).Dexp(nn).RFstdG(j);
  123. elseif Exp(expN).Dexp(nn).meanP(j,m)==0 && strcmpi(ctrlRF,'F')
  124. Rn(nn)=Exp(expN).Dexp(nn).RFmean(j);
  125. Rs(nn)=Exp(expN).Dexp(nn).RFstd(j);
  126. elseif Exp(expN).Dexp(nn).meanP(j,m)==0 && strcmpi(ctrlRF,'L')
  127. Rn(nn)=Exp(expN).Dexp(nn).RFmeanG(j);
  128. Rs(nn)=Exp(expN).Dexp(nn).RFstdG(j);
  129. end
  130. end
  131. end