Supp6WTCNO.m 587 B

12345678910111213141516171819202122232425262728293031323334
  1. load('CtrlDreaddCumValsNew.mat', 'Ctrl', 'CNO')
  2. %%
  3. Expert = Ctrl(1:3);
  4. WT = setdiff(Ctrl, [CNO;Expert]);
  5. %%
  6. figure
  7. PlotAllDataPoints3(WT, CNO, Expert)
  8. set(gca, 'XTick',1:3, 'XTickLabel', {'WT', 'CNO','Expert'})
  9. ylabel('Normalized learning score')
  10. xlim([0.5 3.5])
  11. %%
  12. x = [WT; CNO; Expert];
  13. group = [ones(length(WT),1); ones(length(CNO),1)+1; ones(length(Expert),1)+2]
  14. p = kruskalwallis(x,group)
  15. %%
  16. [p1 h] = ranksum(WT, CNO)
  17. text(1.2, 0.7, ['p = ', num2str(p1)])
  18. [p2 h] = ranksum(WT, Expert)
  19. text(2.2, 0.7, ['p = ', num2str(p2)])
  20. ylim([-0.5 1])
  21. set(gca,'YTick', [-0.5:0.5:1])