Inference of Median Subjective Threshold in Psychophysical Experiments

We consider the response of a test subject upon a skin area being heated with an electromagnetic wave or a contact surface. When the specifications of the electromagnetic beam are fixed, the stimulus is solely described by the heating duration. The binary response of a subject, escape or no escape, is determined by the stimulus and a subjective threshold that varies among test realizations. We study four methods for inferring the median subjective threshold in psychophysical experiments: 1) sample median, 2) maximum likelihood estimation (MLE) with 2 variables, 3) MLE with 1 variable, and 4) adaptive Bayesian method. While methods 1 - 3 require samples of time to escape measured in the method of limits, method 4 utilizes binary outcomes observed in the method of constant stimuli. We find that a) the adaptive Bayesian method converges and is as efficient as the sample median even when the assumed model distribution is incorrect;b) this robust convergence is lost if we infer the mean instead of the median;c) for the optimal performance in an uncertain situation, it is best to use a wide model distribution;d) the predicted error from the posterior standard deviation is unreliable, dominated by the assumed model distribution.

Extracting Human Reaction Time from Observations in the Method of Constant Stimuli

We consider the psychophysical experiments in which the test subject’s binary reaction is determined by the prescribed exposure duration to a stimulus and a random variable subjective threshold. For example, when a subject is exposed to a millimeter wave beam for a prescribed duration, the occurrence of flight action is binary (yes or no). In experiments, in addition to the binary outcome, the actuation time of flight action is also recorded if it occurs;the delay from the initiation time to the actuation time of flight action is the human reaction time, which is not measurable. In this study, we model the random subjective threshold as a Weibull distribution and formulate an inference method for estimating the human reaction time, from data of prescribed exposure durations, binary outcomes and actuation times of flight action collected in a sequence of tests. Numerical simulations demonstrate that the inference of human reaction time based on the Weibull distribution converges to the correct value even when the underlying true model deviates from the inference model. This robustness of the inference method makes it applicable to real experimental data where the underlying true model is unknown.

Minimum Energy Requirement for Inducing Withdrawal Reflex in Millimeter Wave Exposures

We consider the problem of inducing withdrawal reflex on a test subject via exposure to a millimeter wave beam. In our physical model, there are 10 physical parameters affecting the occurrence of withdrawal reflex. Our goal is to pinpoint the roles of these physical parameters in inducing withdrawal reflex. We first carry out non-dimensionalization to reduce the model to a non-dimensional system of only 3 composite parameters: non-dimensional beam power density, non-dimensional beam radius, and non-dimensional exposure time. If the beam power is kept on and steady, withdrawal reflex occurs eventually;the shortest exposure time for inducing withdrawal reflex corresponds to the smallest energy consumption at the given power density and beam radius. In the 2D space of power density and beam radius, the overall minimum energy occurs at the corner of very large power density and very small beam radius, which also produces a very large value of maximum skin temperature and a long time to withdrawal reflex. To reduce the burn injury risk, we introduce a cap on the maximum skin temperature. At each given total beam power, we carry out optimizations with respect to the beam radius, constrained by the prescribed temperature cap. The energy consumption varies negatively with the prescribed temperature cap: a lower temperature cap can be accommodated only with a higher energy consumption via the venue of a larger beam radius. The energy consumption is relatively flat with respect to the total beam power and attains a minimum at a moderately large total beam power. The time to withdrawal reflex is approximately inversely proportional to the total beam power. Our analysis demonstrates that a moderately large total beam power is a good compromise to achieve both low energy consumption and short time to withdrawal reflex.