[E,cues,cueLabels] = baumgartner2021( target,template )

target |
binaural impulse response(s) referring to the directional transfer function(s) (DFTs) of the target sound(s). Option 1: given in SOFA format -> sagittal plane DTFs will be extracted internally. Option 2: binaural impulse responses of all available listener-specific DTFs of the sagittal plane formatted according to the following matrix dimensions: time x direction x channel/ear |

template |
binaural impulse responses of all available
listener-specific DTFs of the sagittal plane referring to
the perceived lateral angle of the target sound.
Options 1 & 2 equivalent to target. |

E |
predicted degree of externalization |

cues |
outcomes of individual cues |

cueLabels |
cue labels; cell array with 1st col. denoting acronyms and 2nd col. for descriptions |

`baumgartner2021(...)` is a model framework for auditory
externalization perception. It enables to probe the contribution of
cue-specific expectation errors and to contrast dynamic versus static
strategies for combining those errors within static listening environments.

`baumgartner2021` accepts the following optional parameters:

'cueWeights',cW |
Set the weights of individual cues to determine the final externalization score. |

Cue-specific weights (entered as a vector) are ordered as follows:

- monaural spectral similarity (MSS)
- interaural spectral similarity of ILDs (ISS)
- spectral standard deviation of monaural gradients (MSSD)
- spectral standard deviation of ILDs (ISSD)
- interaural broadband time-intensity coherence (ITIT)
- interaural coherence (IC)
- monaural intensity difference (MI)
- temporal standard deviation of ILDs (ITSD).
Default weights are 0.6 for MSS, 0.4 for ISS, and 0 for all others.

'S',S |
Set the cue-specific sensitivity parameter to S.
1/S represents the slope of sigmoidal mapping
function.
Vector order equivalent to cueWeights.
Default values are determined by the weighted average
sensitivities determined in Baumgartner and Majdak
(2020) - run exp_baumgartner2021('tab2') to show them. |

'lat',lat |
Set the apparent lateral angle of the target sound to
lat. Default value is 0 degree (median SP). |

'range',c1 |
Set the range factor of the externalization scores to c1.
Default value is 3.78 from Hassager et al. (2016). |

'offset',c2 |
Set the offset of the externalization score to c2.
Default value is 1 from Hassager et al. (2016). |

'ILD_JND',L |
Set the just noticeable ILD difference to L from the
internal template. Default value is 1 (dB). |

'ITD_JND',T |
Set the just noticeable ITD difference to T from the
internal template. Default value is 20e-6 (s). |

`baumgartner2021` accepts the following flags:

'LTA' |
Looser-takes-all strategy: Model selects minimal predicted externalization scores across cues with weights larger than zero. |

'MTA' |
Median-takes-all strategy: Model selects median predicted externalization scores across cues with weights larger than zero. |

'WTA' |
Winner-takes-all strategy: Model selects maximal predicted externalization scores across cues with weights larger than zero. |

- SOFA API from http://sourceforge.net/projects/sofacoustics for Matlab (in e.g. thirdparty/SOFA)
- Data in hrtf/baumgartner2017
- Circular Statistics Toolbox from http://www.mathworks.com/matlabcentral/fileexchange/10676-circular-statistics-toolbox--directional-statistics-