rmeddis@28: function plotProfile(longTone,shortTone,gaps,BFs,TMC,offBFs,IFMCs)
rmeddis@28: 
rmeddis@28: %% plot profile
rmeddis@28: if nargin<1
rmeddis@28:     load profile
rmeddis@28: end
rmeddis@28: 
rmeddis@29: % comparison data (e.g. participants)
rmeddis@29: % rows are BFs
rmeddis@28: 
rmeddis@29: % -------------------------------------------JSan
rmeddis@29: compareBFs=[250	500	1000	2000	3000
rmeddis@29: ];
rmeddis@29: compareLongTone=  [67.46485	56.95655	65.01985	61.46655	73.33265
rmeddis@29: ];
rmeddis@29: compareShortTone=[	72.3185	63.2818	69.0373	65.2853	76
rmeddis@29: ];
rmeddis@28: 
rmeddis@29: compareGaps=[0.02 0.05 0.08];
rmeddis@29: compareTMC=	[
rmeddis@29:     84	69	77	75	93
rmeddis@29: 88	73	81	79	97
rmeddis@29: 95	79	85	83	98
rmeddis@28: ];
rmeddis@29: 
rmeddis@29: compareMaskerFreqs=[0.7  0.9 1 1.1  1.3 ];
rmeddis@29: compareIFMCs=[83.1698	77.3165	79.8474	82.9074	82.3294
rmeddis@29: 80.9667	73.6653	80.9446	80.7005	79.0022
rmeddis@29: 82.0135	71.2284	78.7345	74.3342	84.126
rmeddis@29: 79.3348	70.3347	78.5769	79.7274	91.9849
rmeddis@29: 79.2308	76.3164	83.6881	86.2105	NaN
rmeddis@29:   ];
rmeddis@29: 
rmeddis@29: % % -------------------------------------------JE
rmeddis@29: % BFs=[250 500 1000 2000 4000 8000];
rmeddis@29: % compareLongTone=  [32 30	31	40	54 NaN];
rmeddis@29: % compareShortTone=[	49	50	47	56	63	NaN];
rmeddis@29: % 
rmeddis@29: % compareGaps=0.01:0.01:0.09;
rmeddis@29: % compareTMC=	[
rmeddis@29: %     69	83	82	NaN	NaN	NaN	NaN	NaN	NaN
rmeddis@29: %     61	68	79	88	93	NaN	NaN	NaN	NaN
rmeddis@29: %     63	69	79	84	92	NaN	NaN	NaN	NaN
rmeddis@29: %     67	71	75	80	82	84	88	93	NaN
rmeddis@29: %     82	82	86	86	NaN	83	88	90	75
rmeddis@29: %     NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
rmeddis@29: %     ];
rmeddis@29: % compareTMC=compareTMC';
rmeddis@29: % 
rmeddis@29: % compareMaskerFreqs=[0.5	0.7	0.9	1	1.1	1.3	1.6];
rmeddis@29: % compareIFMCs=[
rmeddis@29: %     64	60	58	58	57	60	64
rmeddis@29: %     65	63	58	55	54	59	69
rmeddis@29: %     68	64	60	59	62	73	79
rmeddis@29: %     76	75	71	67	68	71	77
rmeddis@29: %     79	71	68	69	73	75	77
rmeddis@29: %     76	73	75	75	76	80	NaN
rmeddis@29: %     ];
rmeddis@29: % compareIFMCs=compareIFMCs';
rmeddis@29: 
rmeddis@29: % -------------------------------------------CMR
rmeddis@29: % CMR 
rmeddis@29: % BFs=[250 500 1000 2000 4000 8000];
rmeddis@29: % compareLongTone=[	11.4	1.55	-13.5	-6.35	-6.4	7.45];
rmeddis@29: % compareShortTone=[	23.85	18.9	9.85	10.6	9.55	21.9];
rmeddis@29: % 
rmeddis@29: % compareGaps=0.01:0.01:0.09;
rmeddis@29: % compareTMC=	[
rmeddis@29: %     28.5	35.0	49.3	70.1	80.5	85.5    NaN      NaN    NaN;
rmeddis@29: %     31.1	44.3	48.4	59.5	56.4	76.7	70.2	82.4	76.3;
rmeddis@29: %     33.4	38.4	48.8	55.8	64.5	78.7	84.2	88.3	90.3;
rmeddis@29: %     25.4	37.0	49.2	49.7	58.2	69.6	87.7	95.8	93.0;
rmeddis@29: %     18.2	23.5	27.4	41.5	64.3	82.1	86.7	91.2	NaN;
rmeddis@29: %     32.5	35.8	43.5	52.1	69.1	78.6	86.6	86.0	NaN;
rmeddis@29: %     ];
rmeddis@29: % compareTMC=compareTMC';
rmeddis@29: % 
rmeddis@29: % compareMaskerFreqs=[0.5	0.7	0.9	1	1.1	1.3	1.6];
rmeddis@29: % compareIFMCs=[
rmeddis@29: %     50	42	34	35	34	33	37;
rmeddis@29: %     58	51	38	33	28	41	49;
rmeddis@29: %     57	41	27	20	28	37	66;
rmeddis@29: %     61	49	27	20	34	68	79;
rmeddis@29: %     67	45	27	22	46	74	87;
rmeddis@29: %     62	62	43	22	47	56	83;
rmeddis@29: %     ];
rmeddis@29: % compareIFMCs=compareIFMCs';
rmeddis@28: 
rmeddis@28: % absolute thresholds
rmeddis@28: figure(90), clf
rmeddis@28: subplot(2,1,2)
rmeddis@28: semilogx(BFs,longTone,'ko-','lineWidth',2); hold on
rmeddis@28: semilogx(BFs,shortTone,'bo-','lineWidth',2); hold on
rmeddis@29: semilogx(compareBFs,compareLongTone,'ko:'); hold on
rmeddis@29: semilogx(compareBFs,compareShortTone,'bo:'); hold on
rmeddis@28: ylim([0 100])
rmeddis@28: 
rmeddis@28: % TMC
rmeddis@28: for BFno=1:length(BFs)
rmeddis@28:     subplot(2,6,BFno)
rmeddis@28:     plot(gaps,TMC(:,BFno)-longTone(BFno),'r','lineWidth',3), hold on
rmeddis@29:     plot(gaps,TMC(:,BFno),'b','lineWidth',3), hold on
rmeddis@29:     ylim([-10 110])
rmeddis@28:     xlim([0.01 0.1])
rmeddis@29:     grid on
rmeddis@28:     if BFno==1
rmeddis@28:         ylabel('masker dB SL')
rmeddis@28:         xlabel('gap')
rmeddis@29: %         text(0.02,80,' TMC','backgroundColor','w')
rmeddis@28:     end
rmeddis@29:     title([num2str(BFs(BFno)) ' Hz'])
rmeddis@29:     set(gca,'XTick',[ 0.1],'xTickLabel', { '0.1'})
rmeddis@28: end
rmeddis@28: 
rmeddis@28: % IFMCs
rmeddis@28: for BFno=1:length(BFs)
rmeddis@28:     freq=offBFs'*BFs(BFno);
rmeddis@28:     subplot(2,1,2)
rmeddis@28:     semilogx(freq,IFMCs(:,BFno),'r','lineWidth',3), hold on
rmeddis@28:     ylim([0 100])
rmeddis@28:     xlim([100 12000])
rmeddis@28:     grid on
rmeddis@28: end
rmeddis@28: xlabel('frequency (Hz)')
rmeddis@28: ylabel('masker dB / probe dB')
rmeddis@28: set(gca,'XTick',BFs)
rmeddis@29: 
rmeddis@29: for BFno=1:length(compareBFs)
rmeddis@29:     subplot(2,6,BFno)
rmeddis@29:     plot(compareGaps,compareTMC(:,BFno)-longTone(BFno),'k:')
rmeddis@29:     plot(compareGaps,compareTMC(:,BFno),'k:')
rmeddis@29:     ylim([-10 110])
rmeddis@29:     xlim([0.01 0.1])
rmeddis@29:     grid on
rmeddis@29:     if BFno==1
rmeddis@29:         ylabel('masker dB SL')
rmeddis@29:         xlabel('gap')
rmeddis@29: %         text(0.02,80,' TMC','backgroundColor','w')
rmeddis@29:     end
rmeddis@29:     title([num2str(BFs(BFno)) ' Hz'])
rmeddis@29:     set(gca,'XTick',[ 0.1],'xTickLabel', { '0.1'})
rmeddis@29: end
rmeddis@29: 
rmeddis@29: % IFMCs
rmeddis@29: for BFno=1:length(compareBFs)
rmeddis@29:     compareFreq=compareMaskerFreqs'*BFs(BFno);
rmeddis@29:     subplot(2,1,2)
rmeddis@29:     semilogx(compareFreq,compareIFMCs(:,BFno),'k:')
rmeddis@29:     ylim([0 100])
rmeddis@29:     xlim([100 12000])
rmeddis@29:     grid on
rmeddis@29: end