We fitted the psychometric functions based on raw responses of an individual subject in a given experimental session as follows. (No conditions were shared between sessions.) We denote by
tin the duration of the test stimulus (the fifth stimulus in Experiments 1.1, 2, and 3, and 5 and the first stimulus in Experiment 4) on the
nth trial in the
ith experimental condition. We assume that the probability of the participant's response
rn on that trial is
where
λ is the probability that the participant guesses (lapse rate; can depend on session; Wichmann & Hill,
2001),
μi is the test duration that a participant perceives as equally long as the other stimuli (point of subjective equality; PSE),
σi reflects the participants' sensitivity (just noticeable difference; JND), and Φ(
tin;
μi,
σi) is the cumulative Gaussian distribution function with mean
μi and standard deviation
σi. We use boldface
μ and
σ to denote the vectors of PSE and JNDs across all conditions tested in the session. To calculate the likelihood of the parameters,
L(
μ,
σ,
λ) =
p(data|
μ,
σ,
λ), we assume conditional independence between trials, allowing us to multiply across trials the factors
p(
rn|
tin,
μi,
σi,
λ) that correspond to the participant's responses:
where
C is the number of conditions in the session and
Ni is the number of trials in condition
i. For each session separately, the parameters
μ,
σ, and
λ were estimated simultaneously to maximize log
L(
μ,
σ,
λ), using fmincon in Matlab. For example, if an experiment had two conditions in a session, we simultaneously fitted five parameters for that session:
μ1,
μ2,
σ1,
σ2, and
λ. We quantify the relative duration distortion (RDD) of the test stimulus by the equation RDD
= (
tref-PSE)
/PSE, where
tref is the duration of the first four stimuli (the last four in Experiment 4). For example, an RDD of 0.05 means that the last stimulus was judged to be 5% longer than the other stimuli. All parameter estimates of each experiment are reported in
Supplementary Tables S1 through
S5.