Comparison of castleCSF with chromatic Pyramid of visibility as functions of spatial and temporal frequency, luminance, and colour directions with lower spatial and temporal frequencies removed

Fitting error

Dataset Fitting error Sensitivity adjustment
castleCSF Combined Chromatic Pyramid of Visibility CSF castleCSF Combined Chromatic Pyramid of Visibility CSF
Average training 3.45 ± 0.08 [dB] 4.73 ± 0.08 [dB] N/A N/A
Average testing 3.55 ± 0.17 [dB] 4.79 ± 0.36 [dB] N/A N/A
modelfest 2.51 [dB] 2.08 [dB] 0.876 1.414
hdrvdp_csf 2.74 [dB] 5.76 [dB] 1.127 1.861
rovamo1993 2.79 [dB] 7.33 [dB] 1.742 2.696
laird2006 5.22 [dB] 5.63 [dB] 0.904 1.526
snowden1995 4.11 [dB] 5.86 [dB] 1.233 2.527
robson1966 2.17 [dB] 2.98 [dB] 1.123 3.188
virsu1979 5.93 [dB] 10.83 [dB] 1.546 2.860
virsu1982 1.60 [dB] 2.07 [dB] 1.003 1.891
wright1983 2.74 [dB] 2.01 [dB] 0.963 1.450
colorfest 2.59 [dB] 5.30 [dB] 0.824 2.134
hdr_csf 2.28 [dB] 10.82 [dB] 1.000 1.000
kim2013 4.11 [dB] 4.38 [dB] 0.887 2.336
five_centres 3.53 [dB] 3.36 [dB] 0.929 4.914
lucassen2018 2.33 [dB] 2.33 [dB] 0.488 1.689
hdr_csf_disc 2.77 [dB] 5.24 [dB] 1.584 1.807
kong2018 3.97 [dB] 2.45 [dB] 1.073 4.726
vanderHorst1969_b 3.18 [dB] 3.71 [dB] 0.626 2.455

Model comparison statistics

Model No. of free parameters Sum of Square Errors (SS) Degrees of freedom (df) F-test AIC
F-statistic p-value
castleCSF (Reference Model) 53 42.644 1315 N/A N/A -4638.52
Combined Chromatic Pyramid of Visibility CSF 15 77.166 1353 28.0142 0.0000 ✓ -3903.2

We use AIC and F-test to test whether the difference in fitting error is statistically significant at alpha=0.05 level. Both statistical metrics take the number of optimized parameters into account.

F-test: For F-test, we compare the fitting results from castleCSF with those of other models. The F-static is calculated using the residual sum of squares and degrees of freedom (number of data points - number of optimized parameters) from both models. The corresponding p-value indicates whether or not the null hypothesis is rejected, where H0: the castleCSF does not provide significant better fit than the other model. The p-values less than 0.05 indicates that castleCSF provides a better fit to the data at the significance level of 0.05 (marked with ✓). We performed the F-test for all individual datasets as well as for all datasets combined. For smaller datasets, where the number of data points are comparable to the number of model parameters, F-test can not provide any results since it indicates there is more variance within the models' fits than between.

AIC: Akaike information criterion is a statistical estimator of prediction error and relative quality of the models, which accounts for the number of parameters of each model. The model with the lower AIC score is considered to be better and with a good balance of error value and the number of parameters.

The sensitivity adjustment column contains a multiplier that is used to adjust the sensitivity of each datasets. It corresponds to sd in the paper.

Model parameters

castleCSF
M_lms2acc = 1.0000 1.0000 0 1.0000 -2.3112 0 -1.0000 -1.0000 50.9875 p.rg.sigma_sust = 14.0413; p.rg.beta_sust = 1.11202; p.rg.ch_sust.S_max = [ 238.852 47.4726 0.507981 ]; p.rg.ch_sust.f_max = 0.341427; p.rg.ch_sust.bw = 1.261; p.rg.A_0 = 1717.92; p.rg.f_0 = 0.100677; p.rg.ecc_drop = 0.0591431; p.rg.ecc_drop_nasal = 2.89648e-05; p.rg.ecc_drop_f = 2.04986e-69; p.rg.ecc_drop_f_nasal = 0.180118; p.yv.sigma_sust = 7.06754; p.yv.beta_sust = 0.952515; p.yv.ch_sust.S_max = [ 259.921 61.0117 0.435826 ]; p.yv.ch_sust.f_max = 0.00146478; p.yv.ch_sust.bw = 2.93857; p.yv.A_0 = 2.77056e+07; p.yv.f_0 = 0.000611985; p.yv.ecc_drop = 0.00357397; p.yv.ecc_drop_nasal = 5.85804e-141; p.yv.ecc_drop_f = 0.0080878; p.yv.ecc_drop_f_nasal = 0.0147658; p.ach.ach_sust.S_max = [ 48.3078 6.37828 0.025689 6.34061e-07 1.13679e+10 ]; p.ach.ach_sust.f_max = [ 1.69023 38.5445 0.29831 ]; p.ach.ach_sust.bw = 0.000215444; p.ach.ach_sust.a = 0.054946; p.ach.ach_sust.A_0 = 157.103; p.ach.ach_sust.f_0 = 0.702338; p.ach.ach_trans.S_max = [ 0.175725 1620.11 ]; p.ach.ach_trans.f_max = 0.00300867; p.ach.ach_trans.bw = 2.20582; p.ach.ach_trans.a = 0.000241177; p.ach.ach_trans.A_0 = 1.23754; p.ach.ach_trans.f_0 = 4.1154; p.ach.sigma_trans = 0.104299; p.ach.sigma_sust = 6.14951; p.ach.omega_trans_sl = 1.84581; p.ach.omega_trans_c = 3.42145; p.ach.ecc_drop = 0.0259781; p.ach.ecc_drop_nasal = 0.0452708; p.ach.ecc_drop_f = 0.0217926; p.ach.ecc_drop_f_nasal = 0.0068348; Parameters for Ach component: p.ach_sust.S_max = [ 48.3078 6.37828 0.025689 6.34061e-07 1.13679e+10 ]; p.ach_sust.f_max = [ 1.69023 38.5445 0.29831 ]; p.ach_sust.bw = 0.000215444; p.ach_sust.a = 0.054946; p.ach_trans.S_max = [ 0.175725 1620.11 ]; p.ach_trans.f_max = 0.00300867; p.ach_trans.bw = 2.20582; p.ach_trans.a = 0.000241177; p.ach_trans.A_0 = 1.23754; p.ach_trans.f_0 = 4.1154; p.sigma_trans = 0.104299; p.sigma_sust = 6.14951; p.omega_trans_sl = 1.84581; p.omega_trans_c = 3.42145; p.ecc_drop = 0.0259781; p.ecc_drop_nasal = 0.0452708; p.ecc_drop_f = 0.0217926; p.ecc_drop_f_nasal = 0.0068348; Parameters for RG component: p.ch_sust.S_max = [ 238.852 47.4726 0.507981 ]; p.ch_sust.f_max = 0.341427; p.ch_sust.bw = 1.261; p.A_0 = 1717.92; p.f_0 = 0.100677; p.sigma_sust = 14.0413; p.beta_sust = 1.11202; p.ecc_drop = 0.0591431; p.ecc_drop_nasal = 2.89648e-05; p.ecc_drop_f = 2.04986e-69; p.ecc_drop_f_nasal = 0.180118; Parameters for YV component: p.ch_sust.S_max = [ 259.921 61.0117 0.435826 ]; p.ch_sust.f_max = 0.00146478; p.ch_sust.bw = 2.93857; p.A_0 = 2.77056e+07; p.f_0 = 0.000611985; p.sigma_sust = 7.06754; p.beta_sust = 0.952515; p.ecc_drop = 0.00357397; p.ecc_drop_nasal = 5.85804e-141; p.ecc_drop_f = 0.0080878; p.ecc_drop_f_nasal = 0.0147658;
Combined Chromatic Pyramid of Visibility CSF
p.P_ach = 0.681333; p.c_w_neg_ach = 0.0396845; p.c_f_neg_ach = 0.0538477; p.c_i_ach = 0.403671; p.c_0_ach = 0.552272; p.P_rg = 0.895066; p.c_w_neg_rg = 0.0436539; p.c_f_neg_rg = 0.0544095; p.c_i_rg = 0.374631; p.c_0_rg = 0.402751; p.P_yv = 0.0198153; p.c_w_neg_yv = 0.0472675; p.c_f_neg_yv = 0.122009; p.c_i_yv = 0.32123; p.c_0_yv = 0.0620767; Parameters for Ach subclass Parameters for RG subclass Parameters for YV subclass

CSF model: castleCSF

L+M
L-M
S-(L+M)

CSF model: Combined Chromatic Pyramid of Visibility CSF

L+M
L-M
S-(L+M)

Legend

To keep the plots legible, only up to 3 models are plotted.

Dataset: [modelfest] ModelFest

Achromatic CSF as a function of frequency

Dataset: [hdrvdp_csf] HDR-VDP CSF

Achromatic CSF as a function of frequency

Achromatic CSF as a function of size

Dataset: [rovamo1993] Rovamo et al. 1993

CSF as the funcation of stimulus area

CSF as the function of spatial frequency

Dataset: [laird2006] Laird et al. 2006

Achromatic CSF as a function of temporal frequency for different spatial frequencies

Dataset: [snowden1995] Snowden et al. 1995

Temporal contrast sensitivity at different spatial frequencies and luminance levels

Dataset: [robson1966] Robson 1966

Spatial CSF for different temporal frequencies

Temporal CSF for different spatial frequencies

Dataset: [virsu1979] Virsu & Rovamo 1979

Contrast sensitivity of central and peripheral vision as a function of spatial frequency and eccentricity

Dataset: [virsu1982] Virsu et al. 1982

Contrast sensitivity as the function of frequency

Dataset: [wright1983] Wright and Johnson 1983

CSF as function of eccentricity

Dataset: [colorfest] ColorFest

Chromatic CSF as a function of frequency

Dataset: [hdr_csf] High Dynamic Range CSF

CSF as the function of frequency at different luminance levels (fixed number of cycles)

Dataset: [kim2013] Chromatic CSF [Kim et al. 2013]

Chromatic CSF as a function of frequency

Dataset: [five_centres] Five centres [Xu et al. 2020]

Chromatic CSF as ellipses

Chromatic CSF as a function of frequency

Dataset: [lucassen2018] Lucassen et al. 2018

Lucassen2018 CSF sensitivity vs spatial frequency

Dataset: [hdr_csf_disc] High Dynamic Range Disc CSF

CSF as the function of size at different luminance levels

Dataset: [kong2018] Kong et al. 2018

Chromatic CSF as a function of temporal frequency (9 background colors and 4 color directions in uv space

Dataset: [vanderHorst1969_b] Van der Horst & Bouman 1969

Static chromatic contrast thresholds as a function of spatial frequency

Static chromatic contrast thresholds as function of luminance

Chromatic contrast thresholds for travelling sine waves as function of spatial frequency

Chromatic contrast thresholds for travelling sine waves as function of temporal frequency