minor update

functions/scriptsToMakeFigs/makeFigureS6.m

@@ -18,40 +18,3 @@ param_hi = PPIrisky_high;
 violinMOtwoGroups(param_lo,param_hi,'r')
 [~,p,~,stats] = ttest2(param_lo,param_hi);
 disp(['risky PPI (in dmPFC ROI) between groups: t(' num2str(stats.df) ') = ' num2str(stats.tstat) ', P = ' num2str(round(p,4,'significant'))])
-
-
-test=[OCUsafe_trans([1:4,6:31]) PPIsafe];
-[r,p] = corr(test(:,1),test(:,2));
-disp(['correlation between OCUsafe and PPI safe: r = ' num2str(r) ', P = ' num2str(p)])
-test=test(splitSubstNaive([1:4,6:31]),:);
-test2=test(~isnan(test(:,1)) & ~isnan(test(:,2)),:);
-stats=regstats(test2(:,2), test2(:,1));
-figure; hold on; plot(test2(:,1), test2(:,2),'.','MarkerSize',smallMarkerSize,'color',colorBlueDrugNo)
-plot([min(test2(:,1)) max(test2(:,1))], stats.tstat.beta(1) + stats.tstat.beta(2)*[min(test2(:,1)) max(test2(:,1))],'color',colorBlueDrugNo,'LineWidth',2)
-
-test=[OCUsafe_trans([1:4,6:31]) PPIsafe];
-test=test(splitSubstExposed([1:4,6:31]),:);
-test2=test(~isnan(test(:,1)) & ~isnan(test(:,2)),:);
-stats=regstats(test2(:,2), test2(:,1));
-plot(test2(:,1), test2(:,2),'.','MarkerSize',smallMarkerSize,'color',colorBlueDrugYes)
-plot([min(test2(:,1)) max(test2(:,1))], stats.tstat.beta(1) + stats.tstat.beta(2)*[min(test2(:,1)) max(test2(:,1))],'color',colorBlueDrugYes,'LineWidth',2)
-xlim([0 1])
-plot([0 1],[0 0],'k--')
-
-test=[OCUrisky_trans(1:31) PPIrisky];
-[r,p] = corr(test(:,1),test(:,2));
-disp(['correlation between OCUrisky and PPI risky: r = ' num2str(r) ', P = ' num2str(p)])
-test=test(splitSubstNaive(1:31),:);
-test2=test(~isnan(test(:,1)) & ~isnan(test(:,2)),:);
-stats=regstats(test2(:,2), test2(:,1));
-figure; hold on; plot(test2(:,1), test2(:,2),'.','MarkerSize',smallMarkerSize,'color',colorRedDrugNo)
-plot([min(test2(:,1)) max(test2(:,1))], stats.tstat.beta(1) + stats.tstat.beta(2)*[min(test2(:,1)) max(test2(:,1))],'color',colorRedDrugNo,'LineWidth',2)
-
-test=[OCUrisky_trans(1:31) PPIrisky];
-test=test(splitSubstExposed(1:31),:);
-test2=test(~isnan(test(:,1)) & ~isnan(test(:,2)),:);
-stats=regstats(test2(:,2), test2(:,1));
-plot(test2(:,1), test2(:,2),'.','MarkerSize',smallMarkerSize,'color',colorRedDrugYes)
-plot([min(test2(:,1)) max(test2(:,1))], stats.tstat.beta(1) + stats.tstat.beta(2)*[min(test2(:,1)) max(test2(:,1))],'color',colorRedDrugYes,'LineWidth',2)
-xlim([0 1])
-plot([0 1],[0 0],'k--')