-
Extreme-value Statistics: Rudiments and applications
Authors:
Evangelos Matsinos
Abstract:
This study provides a summary of the theory which enables the analysis of extreme values, i.e., of measurements acquired from the observation of extraordinary/rare physical phenomena. The formalism is developed in a transparent way, tailored to the application to real-world problems. Three examples of the application of the theory are detailed: first, an old problem, relating to the return period…
▽ More
This study provides a summary of the theory which enables the analysis of extreme values, i.e., of measurements acquired from the observation of extraordinary/rare physical phenomena. The formalism is developed in a transparent way, tailored to the application to real-world problems. Three examples of the application of the theory are detailed: first, an old problem, relating to the return period of floods caused by the Rh{ô}ne river, is revisited; second, the frequency of occurrence and the severity of meteorite impacts on the Earth are examined; the data, which are analysed in the third example, relate to the longevity of supercentenarians.
△ Less
Submitted 30 June, 2024;
originally announced July 2024.
-
Modelling of the monthly global average temperature anomalies after the year 1900
Authors:
Evangelos Matsinos
Abstract:
Two ways of parameterising the temporal dependence of the global average temperature anomalies after 1900 are put forth in this technical note. The models are fitted to the data of the Berkeley Earth Group (up to the end of 2023), after the application of corrections which translate the input so that its average value over the last five decades of the $19$-th century vanishes; as the corrections a…
▽ More
Two ways of parameterising the temporal dependence of the global average temperature anomalies after 1900 are put forth in this technical note. The models are fitted to the data of the Berkeley Earth Group (up to the end of 2023), after the application of corrections which translate the input so that its average value over the last five decades of the $19$-th century vanishes; as the corrections are evaluated separately for each month of the year, the seasonal effects are removed from the data. The results are given in a convenient form for use.
△ Less
Submitted 27 February, 2024; v1 submitted 23 January, 2024;
originally announced February 2024.
-
Analysis of four sets of global average temperature anomalies
Authors:
Evangelos Matsinos
Abstract:
Four sets of global average temperature anomalies, altered so that they refer to pre-industrial temperature levels (baseline), as recommended by the Intergovernmental Panel on Climate Change, are analysed in this study. Expectation values of the time, at which the crossing of the $1.5^\circ$C threshold above the baseline will occur, are extracted from the two datasets which are accompanied by mean…
▽ More
Four sets of global average temperature anomalies, altered so that they refer to pre-industrial temperature levels (baseline), as recommended by the Intergovernmental Panel on Climate Change, are analysed in this study. Expectation values of the time, at which the crossing of the $1.5^\circ$C threshold above the baseline will occur, are extracted from the two datasets which are accompanied by meaningful uncertainties; the two estimates are in the 2030s and differ by about four years: mid-2032 and mid-2036.
△ Less
Submitted 22 October, 2023;
originally announced October 2023.
-
A discussion on two celebrated examples of application of linear regression
Authors:
Evangelos Matsinos
Abstract:
This work aims at providing some (further) perspective into the analysis of two well-known datasets in terms of the methods of linear regression. The first set has been taken from Hubble's 1929 investigation of the relation between the distances of galaxies and their recessional velocities. The second set relates to Galton's family data on human stature, collected for the purposes of his 1889 book…
▽ More
This work aims at providing some (further) perspective into the analysis of two well-known datasets in terms of the methods of linear regression. The first set has been taken from Hubble's 1929 investigation of the relation between the distances of galaxies and their recessional velocities. The second set relates to Galton's family data on human stature, collected for the purposes of his 1889 book on natural inheritance.
△ Less
Submitted 2 October, 2023;
originally announced October 2023.
-
COVID-19: The extraction of the effective reproduction number from the time series of new cases
Authors:
Evangelos Matsinos
Abstract:
Addressed in this work is the performance of five popular algorithms, which aim at assessing the dissemination dynamics of the COVID-19 disease on the basis of the time series of new confirmed cases. The tests are based on simulated data, generated by means of a deterministic compartmental epidemiological model \cite{Matsinos2020a}, adapted herein to also include the possibility of the loss of imm…
▽ More
Addressed in this work is the performance of five popular algorithms, which aim at assessing the dissemination dynamics of the COVID-19 disease on the basis of the time series of new confirmed cases. The tests are based on simulated data, generated by means of a deterministic compartmental epidemiological model \cite{Matsinos2020a}, adapted herein to also include the possibility of the loss of immunity by the group of the recovered (or vaccinated) subjects. Assuming a simple temporal dependence of the effective reproduction number (the exact details are of no relevance as far as the conclusions of this work are concerned), time series of new cases were generated in a time domain of nearly one year for the five top-ranking countries in the cumulative number of infections by January 1, 2021. These countries are (in descending order of infections): the United States of America, India, Brazil, Russia, and the United Kingdom. The processing of each simulated time series led to the establishment of relations between the input (actual) and the reconstructed values of the effective reproduction number for each country and algorithm, separately; this work argues that all five algorithms underestimate the effective reproduction number when the latter exceeds the critical value of $1$. The five algorithms were subsequently applied to the real-life time series of new cases for the aforementioned five countries, which also span a temporal interval of nearly one year; corrected values of the effective reproduction number are obtained for these countries in 2020.
△ Less
Submitted 7 May, 2021;
originally announced May 2021.
-
COVID-19: On the quarantine duration after short visits to high-risk regions
Authors:
Evangelos Matsinos
Abstract:
A simple Monte-Carlo method will be put forward herein, to enable the extraction of an estimate for the quarantine duration, applicable to visitors to high-risk regions. Results will be obtained on the basis of an analysis of the upper tail of the cumulative distribution function of the time span between the departure of the travellers from the place where the infection occurs and the time instant…
▽ More
A simple Monte-Carlo method will be put forward herein, to enable the extraction of an estimate for the quarantine duration, applicable to visitors to high-risk regions. Results will be obtained on the basis of an analysis of the upper tail of the cumulative distribution function of the time span between the departure of the travellers from the place where the infection occurs and the time instant when COVID-19 infections may currently be detected. As expected, the quarantine duration is a decreasing function of the fraction of the infected travellers, which one is prepared to identify as `acceptable risk'. The analysis suggests that a maximal $5~\%$ risk (of new infections originating from subjects who become infective after their quarantine is lifted) may be associated with a minimal quarantine duration of about eight days, $1~\%$ with about twelve, and $0.1~\%$ with about sixteen. Unless the distribution of the duration of short (typically, up to three weeks) travels departs significantly from the shape assumed in this study, the results of the present analysis do not provide support for the plans to shorten the quarantine duration of about ten days to two weeks, which currently applies to travellers entering most European countries from regions with a high risk of infection.
△ Less
Submitted 6 October, 2020;
originally announced October 2020.
-
Analysis of luminosity measurements of the pre-white dwarf PG 1159-035
Authors:
Evangelos Matsinos
Abstract:
The study of the luminosity measurements of the pre-white dwarf PG 1159-035 has established the properties of the rich power spectrum of the detected radiation and, derived thereof, the physical properties of this celestial body. Those of the measurements which are available online are analysed in this work from a different perspective. After the measurements were band-passed, they were split into…
▽ More
The study of the luminosity measurements of the pre-white dwarf PG 1159-035 has established the properties of the rich power spectrum of the detected radiation and, derived thereof, the physical properties of this celestial body. Those of the measurements which are available online are analysed in this work from a different perspective. After the measurements were band-passed, they were split into two parts (of comparable sizes), one yielding the training (learning) set (i.e., the database of embedding vectors and associated predictions), the other the test set. The optimal embedding dimension $m_0=10$ was obtained using Cao's method; this result was confirmed by an analysis of the correlation dimension. Subsequently, the extraction of the maximal Lyapunov exponent $λ$ was pursued for embedding dimensions $m$ between $3$ and $12$; results were obtained after removing the prominent undulations of the out-of-sample prediction-error arrays $S (k)$ by fitting a monotonic function to the data. The grand mean of the values, obtained for sufficient embedding dimensions ($10 \leq m \leq 12$), was: $λ= (9.2 \pm 1.0 ({\rm stat.}) \pm 2.7 ({\rm syst.})) \cdot 10^{-2}~Δτ^{-1}$, where $Δτ=10$ s is the sampling interval in the measurements. On the basis of this significantly non-zero result, it may be concluded that the physical processes, underlying the variation of the luminosity of PG 1159-035, are non-linear. The aforementioned result for $λ$ was obtained using the $L^\infty$-norm distance; a larger, yet not incompatible, result was extracted with the Euclidean ($L^2$-norm) distance.
△ Less
Submitted 26 August, 2018; v1 submitted 12 August, 2018;
originally announced August 2018.
-
Processing of acoustical signals via a wavelet-based analysis
Authors:
Evangelos Matsinos
Abstract:
In the present paper, details are given on the implementation of a wavelet-based analysis tailored to the processing of acoustical signals. The family of the suitable wavelets (`Reimann wavelets') are obtained in the time domain from a Fourier transform, extracted in Ref.~\cite{r1} after invoking theoretical principles and time-frequency localisation constraints. A scheme is set forth to determine…
▽ More
In the present paper, details are given on the implementation of a wavelet-based analysis tailored to the processing of acoustical signals. The family of the suitable wavelets (`Reimann wavelets') are obtained in the time domain from a Fourier transform, extracted in Ref.~\cite{r1} after invoking theoretical principles and time-frequency localisation constraints. A scheme is set forth to determine the optimal values of the parameters of this type of wavelet on the basis of the goodness of the reproduction of a $30$-s audio file containing harmonic signals corresponding to six successive $A$ notes of the chromatic musical scale, from $A_2$ to $A_7$. The quality of the reproduction over about six and a half octaves is investigated. Finally, details are given on the incorporation of the re-assignment method in the analysis framework, as the means a) to determine the important contributions of the wavelet transforms and b) to suppress noise present in the signal.
△ Less
Submitted 30 September, 2015;
originally announced September 2015.
-
On using the Microsoft Kinect$^{\rm TM}$ sensors to determine the lengths of the arm and leg bones of a human subject in motion
Authors:
M. J. Malinowski,
E. Matsinos
Abstract:
The present study is part of a broader programme, exploring the possibility of involving the Microsoft Kinect$^{\rm TM}$ sensor in the analysis of human motion. We examine the output obtained from the two available versions of this sensor in relation to the variability of the estimates of the lengths of eight bones belonging to the subject's extremities: of the humerus (upper arm), ulna (lower arm…
▽ More
The present study is part of a broader programme, exploring the possibility of involving the Microsoft Kinect$^{\rm TM}$ sensor in the analysis of human motion. We examine the output obtained from the two available versions of this sensor in relation to the variability of the estimates of the lengths of eight bones belonging to the subject's extremities: of the humerus (upper arm), ulna (lower arm, forearm), femur (upper leg), and tibia (lower leg, shank). Large systematic effects in the output of the two sensors have been observed.
△ Less
Submitted 2 May, 2015;
originally announced May 2015.
-
Comparative study of the two versions of the Microsoft Kinect$^{\rm TM}$ sensor in regard to the analysis of human motion
Authors:
M. J. Malinowski,
E. Matsinos
Abstract:
The present paper is part of a broader programme, exploring the possibility of involving the Microsoft Kinect$^{\rm TM}$ sensor in the analysis of human motion. In this study, the output obtained from the two available versions of this sensor is critically examined. We demonstrate that the two outputs differ in regard to the variation of the physical quantities involved in the modelling of the hum…
▽ More
The present paper is part of a broader programme, exploring the possibility of involving the Microsoft Kinect$^{\rm TM}$ sensor in the analysis of human motion. In this study, the output obtained from the two available versions of this sensor is critically examined. We demonstrate that the two outputs differ in regard to the variation of the physical quantities involved in the modelling of the human motion. As the original sensor has been found unsuitable for applications requiring high precision, the observed differences in the output of the two sensors call for the validation of the upgraded sensor on the basis of a marker-based system.
△ Less
Submitted 1 April, 2015;
originally announced April 2015.
-
On using the Microsoft Kinect$^{\rm TM}$ sensors in the analysis of human motion
Authors:
M. J. Malinowski,
E. Matsinos,
S. Roth
Abstract:
The present paper aims at providing the theoretical background required for investigating the use of the Microsoft Kinect$^{\rm TM}$ (`Kinect', for short) sensors (original and upgraded) in the analysis of human motion. Our methodology is developed in such a way that its application be easily adaptable to comparative studies of other systems used in capturing human-motion data. Our future plans in…
▽ More
The present paper aims at providing the theoretical background required for investigating the use of the Microsoft Kinect$^{\rm TM}$ (`Kinect', for short) sensors (original and upgraded) in the analysis of human motion. Our methodology is developed in such a way that its application be easily adaptable to comparative studies of other systems used in capturing human-motion data. Our future plans include the application of this methodology to two situations: first, in a comparative study of the performance of the two Kinect sensors; second, in pursuing their validation on the basis of comparisons with a marker-based system (MBS). One important feature in our approach is the transformation of the MBS output into Kinect-output format, thus enabling the analysis of the measurements, obtained from different systems, with the same software application, i.e., the one we use in the analysis of Kinect-captured data; one example of such a transformation, for one popular marker-placement scheme (`Plug-in Gait'), is detailed. We propose that the similarity of the output, obtained from the different systems, be assessed on the basis of the comparison of a number of waveforms, representing the variation within the gait cycle of quantities which are commonly used in the modelling of the human motion. The data acquisition may involve commercially-available treadmills and a number of velocity settings: for instance, walking-motion data may be acquired at $5$ km/h, running-motion data at $8$ and $11$ km/h. We recommend that particular attention be called to systematic effects associated with the subject's knee and lower leg, as well as to the ability of the Kinect sensors in reliably capturing the details in the asymmetry of the motion for the left and right parts of the human body. The previous versions of the study have been withdrawn due to the use of a non-representative database.
△ Less
Submitted 27 May, 2015; v1 submitted 4 December, 2014;
originally announced December 2014.
-
Modelling of the motion of a Mecanum-wheeled vehicle
Authors:
Evangelos Matsinos
Abstract:
The present document aims at developing the formalism needed in order to describe the two-dimensional motion of a vehicle with Mecanum wheels, including the effects of resistive friction. The description of recently-acquired experimental data, on the basis of this model, is found satisfactory.
The present document aims at developing the formalism needed in order to describe the two-dimensional motion of a vehicle with Mecanum wheels, including the effects of resistive friction. The description of recently-acquired experimental data, on the basis of this model, is found satisfactory.
△ Less
Submitted 10 November, 2012;
originally announced November 2012.
-
A method for the determination of nuclear cross sections of proton beams by the collective model and extended nuclear-shell theory
Authors:
W. Ulmer,
E. Matsinos
Abstract:
A recapitulatory analysis of total nuclear cross sections of various nuclei is presented, which yields detailed knowledge on the different physical processes such as potential/resonance scatter and nuclear reactions. The physical base for potential/resonance scatter and the threshold energy resulting from Coulomb repulsion of nuclei are collective/oscillator models. The part pertaining to the nucl…
▽ More
A recapitulatory analysis of total nuclear cross sections of various nuclei is presented, which yields detailed knowledge on the different physical processes such as potential/resonance scatter and nuclear reactions. The physical base for potential/resonance scatter and the threshold energy resulting from Coulomb repulsion of nuclei are collective/oscillator models. The part pertaining to the nuclear reactions can only be determined by the microscopic theory (Schrödinger equation and strong interactions). The physical impact is the fluence decrease of proton beams in different media and the scatter behavior of secondary particles.
△ Less
Submitted 4 September, 2010;
originally announced September 2010.
-
Theoretical methods for the calculation of Bragg curves and 3D distributions of proton beams
Authors:
Waldemar Ulmer,
Evangelos Matsinos
Abstract:
The well-known Bragg-Kleeman rule RCSDA = A dot E0p has become a pioneer work in radiation physics of charged particles and is still a useful tool to estimate the range RCSDA of approximately monoenergetic protons with initial energy E0 in a homogeneous medium. The rule is based on the continuous-slowing-down-approximation (CSDA). It results from a generalized (nonrelativistic) Langevin equation a…
▽ More
The well-known Bragg-Kleeman rule RCSDA = A dot E0p has become a pioneer work in radiation physics of charged particles and is still a useful tool to estimate the range RCSDA of approximately monoenergetic protons with initial energy E0 in a homogeneous medium. The rule is based on the continuous-slowing-down-approximation (CSDA). It results from a generalized (nonrelativistic) Langevin equation and a modification of the phenomenological friction term. The complete integration of this equation provides information about the residual energy E(z) and dE(z)/dz at each position z (0 <= z <= RCSDA). A relativistic extension of the generalized Langevin equation yields the formula RCSDA = A dot (E0 +E02/2M dot c2)p. The initial energy of therapeutic protons satisfies E0 << 2M dot c2 (M dot c2 = 938.276 MeV), which enables us to consider the relativistic contributions as correction terms. Besides this phenomenological starting-point, a complete integration of the Bethe-Bloch equation (BBE) is developed, which also provides the determination of RCSDA, E(z) and dE(z)/dz and uses only those parameters given by the BBE itself (i.e., without further empirical parameters like modification of friction). The results obtained in the context of the aforementioned methods are compared with Monte-Carlo calculations (GEANT4); this Monte-Carlo code is also used with regard to further topics such as lateral scatter, nuclear interactions, and buildup effects. In the framework of the CSDA, the energy transfer from protons to environmental atomic electrons does not account for local fluctuations.
△ Less
Submitted 21 August, 2010;
originally announced August 2010.
-
Collimator effects in proton planning
Authors:
Evangelos Matsinos
Abstract:
The present paper pertains to corrections which are due to the presence of beam-limiting and beam-shaping devices in proton planning. Two types of corrections are considered: those which are due to the nonzero thickness of such devices (geometrical effects) and those relating to the scattering of beam particles off their material. The application of these two types of corrections is greatly faci…
▽ More
The present paper pertains to corrections which are due to the presence of beam-limiting and beam-shaping devices in proton planning. Two types of corrections are considered: those which are due to the nonzero thickness of such devices (geometrical effects) and those relating to the scattering of beam particles off their material. The application of these two types of corrections is greatly facilitated by decomposing the physical effects (i.e., the contribution to the fluence) of two-dimensional objects (i.e., of the apertures of the devices) into one-dimensional, easily-calculable contributions. To minimise the time requirements in the derivation of the scattering corrections, a two-step process is introduced. The first step occurs at beam-configuration phase and comprises the analysis of half-block fluence measurements and the extraction of the one parameter of the model which is used in the description of the beamline characteristics; subsequently, a number of Monte-Carlo runs lead to the determination of the parameters of a convenient parameterisation of the relevant fluence contributions. The second step involves (at planning time) the reconstruction of the parameters (which are used in the description of the scattering contributions) via simple interpolations, performed on the results obtained during the beam-configuration phase. It is shown that the inclusion of the scattering effects leads to substantial improvement in the reproduction of the experimental data. The contributions from the block-thickness and block-scattering effects have been presented separately in the case of a simple water phantom. In this example, the maximal contribution of the block-relating effects amounts to a few percent of the prescribed dose.
△ Less
Submitted 5 December, 2008; v1 submitted 7 November, 2008;
originally announced November 2008.
-
The dual-dose imaging technique: a way to enhance the dynamic range of X-ray detectors
Authors:
Evangelos Matsinos,
Wolfgang Kaissl
Abstract:
We describe a method aiming at increasing the dynamic range of X-ray detectors. Two X-ray exposures of an object are acquired at different dose levels and constitute the only input data. The values of the parameters which are needed to process these images are determined from information contained in the images themselves; the values of two parameters are extracted from the input data. The two i…
▽ More
We describe a method aiming at increasing the dynamic range of X-ray detectors. Two X-ray exposures of an object are acquired at different dose levels and constitute the only input data. The values of the parameters which are needed to process these images are determined from information contained in the images themselves; the values of two parameters are extracted from the input data. The two input images are finally merged in such a way as to create one image containing useful information in all its entirety. This selective use of parts of each image allows both the contour of the irradiated object to be visible and the high-attenuation areas to retain their image quality corresponding to the information contained in the high-dose image. The benefits of the method are demonstrated with an example involving a head phantom.
△ Less
Submitted 4 July, 2006;
originally announced July 2006.
-
The dual-gain mode: a way to enhance the dynamic range of X-ray detectors
Authors:
Evangelos Matsinos,
Wolfgang Kaissl
Abstract:
Varian Medical Systems has manufactured and recently put into operation a clinically-applicable solution for image-guided radiation therapy. Cone-beam CT imaging, one of the operation modes of the imaging unit of this device, aims at high-quality volumetric reconstruction. To boost the image quality, the dual-gain mode, a successful means for enhancing the dynamic range of the flat-panel detecto…
▽ More
Varian Medical Systems has manufactured and recently put into operation a clinically-applicable solution for image-guided radiation therapy. Cone-beam CT imaging, one of the operation modes of the imaging unit of this device, aims at high-quality volumetric reconstruction. To boost the image quality, the dual-gain mode, a successful means for enhancing the dynamic range of the flat-panel detectors and obtaining better results in the contrast of the reconstructed image, was developed and successfully tested during the last few years. The important steps in the calibration of this mode involve a correction to the pulse widths associated with the X-ray production, the assessment of the detector signal above which nonlinear effects become significant and the determination of some properties of the detector pixels, namely, of dark fields, flatness corrections, etc. Finally, the defect-pixel map is obtained, containing dead and flickering pixels, as well as pixels with properties which are sufficiently `out of range'. An effect observed in the offset correction of raw images seems to originate in the mismatch between the measured and extrapolated values of the dark field; a correction scheme is proposed to take account of this effect. In the data processing, which is achieved on the basis of the files produced in the calibration of the mode, the high-gain signal is used whenever meaningful; otherwise, it is substituted by the low-gain signal, properly scaled. Defect pixels are interpolated from their good neighbours; a procedure achieving a correction in case of small defect-pixel clusters is outlined.
△ Less
Submitted 4 July, 2006;
originally announced July 2006.
-
The geometric calibration of cone-beam imaging and delivery systems in radiation therapy
Authors:
Evangelos Matsinos,
Wolfgang Kaissl
Abstract:
We propose a method to achieve the geometric calibration of cone-beam imaging and delivery systems in radiation therapy; our approach applies to devices where an X-ray source and a flat-panel detector, facing each other, move in circular orbits around the irradiated object. In order to extract the parameters of the geometry from the data, we use a light needle phantom which is easy to manufactur…
▽ More
We propose a method to achieve the geometric calibration of cone-beam imaging and delivery systems in radiation therapy; our approach applies to devices where an X-ray source and a flat-panel detector, facing each other, move in circular orbits around the irradiated object. In order to extract the parameters of the geometry from the data, we use a light needle phantom which is easy to manufacture. A model with ten free parameters (spatial lengths and distortion angles) has been put forth to describe the geometry and the mechanical imperfections of the units being calibrated; a few additional parameters are introduced to account for residual effects (small effects which lie beyond our model). The values of the model parameters are determined from one complete scan of the needle phantom via a robust optimisation scheme. The application of this method to two sets of five counterclockwise (ccw) and five clockwise (cw) scans yielded consistent and reproducible results. A number of differences have been observed between the ccw and cw scans, suggesting a dissimilar behaviour of the devices calibrated in the two rotation modes. The description of the geometry of the devices was achieved at the sub-pixel level.
△ Less
Submitted 4 July, 2006;
originally announced July 2006.