Semi-empirical Likelihood Confidence Intervals for the Differences of Quantiles with Missing Data

Detecting population （group） differences is useful in many applications, such as medical research. In this paper, we explore the probabilistic theory for identifying the quantile differences .between two populations. Suppose that there are two populations x and y with missing data on both of them, where x is nonparametric and y is parametric. We are interested in constructing confidence intervals on the quantile differences of x and y. Random hot deck imputation is used to fill in missing data. Semi-empirical likelihood confidence intervals on the differences are constructed.

Empirical likelihood for spatial cross-sectional data models with matrix exponential spatial specification

In this paper,we study spatial cross-sectional data models in the form of matrix exponential spatial specification(MESS),where MESS appears in both dependent and error terms.The empirical likelihood(EL)ratio statistics are established for the parameters of the MESS model.It is shown that the limiting distributions of EL ratio statistics follow chi-square distributions,which are used to construct the confidence regions of model parameters.Simulation experiments are conducted to compare the performances of confidence regions based on EL method and normal approximation method.

Beamspace maximum likelihood algorithm based on sum and difference beams for elevation estimation

Beamspace super-resolution methods for elevation estimation in multipath environment has attracted significant attention, especially the beamspace maximum likelihood(BML)algorithm. However, the difference beam is rarely used in superresolution methods, especially in low elevation estimation. The target airspace information in the difference beam is different from the target airspace information in the sum beam. And the use of difference beams does not significantly increase the complexity of the system and algorithms. Thus, this paper applies the difference beam to the beamformer to improve the elevation estimation performance of BML algorithm. And the direction and number of beams can be adjusted according to the actual needs. The theoretical target elevation angle root means square error(RMSE) and the computational complexity of the proposed algorithms are analyzed. Finally, computer simulations and real data processing results demonstrate the effectiveness of the proposed algorithms.

Vulnerable brain regions in adolescent major depressive disorder:A resting-state functional magnetic resonance imaging activation likelihood estimation meta-analysis

BACKGROUND Adolescent major depressive disorder(MDD)is a significant mental health concern that often leads to recurrent depression in adulthood.Resting-state functional magnetic resonance imaging(rs-fMRI)offers unique insights into the neural mechanisms underlying this condition.However,despite previous research,the specific vulnerable brain regions affected in adolescent MDD patients have not been fully elucidated.AIM To identify consistent vulnerable brain regions in adolescent MDD patients using rs-fMRI and activation likelihood estimation(ALE)meta-analysis.METHODS We performed a comprehensive literature search through July 12,2023,for studies investigating brain functional changes in adolescent MDD patients.We utilized regional homogeneity(ReHo),amplitude of low-frequency fluctuations(ALFF)and fractional ALFF(fALFF)analyses.We compared the regions of aberrant spontaneous neural activity in adolescents with MDD vs healthy controls(HCs)using ALE.RESULTS Ten studies(369 adolescent MDD patients and 313 HCs)were included.Combining the ReHo and ALFF/fALFF data,the results revealed that the activity in the right cuneus and left precuneus was lower in the adolescent MDD patients than in the HCs(voxel size:648 mm3,P<0.05),and no brain region exhibited increased activity.Based on the ALFF data,we found decreased activity in the right cuneus and left precuneus in adolescent MDD patients(voxel size:736 mm3,P<0.05),with no regions exhibiting increased activity.CONCLUSION Through ALE meta-analysis,we consistently identified the right cuneus and left precuneus as vulnerable brain regions in adolescent MDD patients,increasing our understanding of the neuropathology of affected adolescents.

Optimization of Generator Based on Gaussian Process Regression Model with Conditional Likelihood Lower Bound Search

The noise that comes from finite element simulation often causes the model to fall into the local optimal solution and over fitting during optimization of generator.Thus,this paper proposes a Gaussian Process Regression(GPR)model based on Conditional Likelihood Lower Bound Search(CLLBS)to optimize the design of the generator,which can filter the noise in the data and search for global optimization by combining the Conditional Likelihood Lower Bound Search method.Taking the efficiency optimization of 15 kW Permanent Magnet Synchronous Motor as an example.Firstly,this method uses the elementary effect analysis to choose the sensitive variables,combining the evolutionary algorithm to design the super Latin cube sampling plan;Then the generator-converter system is simulated by establishing a co-simulation platform to obtain data.A Gaussian process regression model combing the method of the conditional likelihood lower bound search is established,which combined the chi-square test to optimize the accuracy of the model globally.Secondly,after the model reaches the accuracy,the Pareto frontier is obtained through the NSGA-II algorithm by considering the maximum output torque as a constraint.Last,the constrained optimization is transformed into an unconstrained optimizing problem by introducing maximum constrained improvement expectation(CEI)optimization method based on the re-interpolation model,which cross-validated the optimization results of the Gaussian process regression model.The above method increase the efficiency of generator by 0.76%and 0.5%respectively;And this method can be used for rapid modeling and multi-objective optimization of generator systems.

QSAR Study of Nitrobenzenes’ Toxicity to Tetrahymena Pyriformis Using Semi-empirical Quantum Chemical Methods

The molecular geometries and electronic structures of 30 nitrobenzenes have been calculated by using semi-empirical MO AM1 and PM3 methods. EHOMO, ELUMO, ENHOMO, ENLUMO, AE, QNO2, Qc and V were selected as the structural descriptors. The acute toxicity （-log/C50） of nitrobenzenes to tetrahymena pyriformis along with the above eight structural parameters was used to establish the quantitative structure-activity relationships （QSARs）. The results indicate that the established model based on AM I method is superior to that on PM3 method not only for the stability but also for the predictive powers of the model. Based on AM1 parameters, a further classifying discussion was presented for the study of nitrobenzene toxic mechanism. The results show that the substituents, nitro group and halogen substituents on the aromatic ring are crucial to the chemicals＇ toxicity. For nitrobenzenes without halogen or other substituent, the reduction of nitro group is the main route. However, for those with halogen substituents, their next lowest unoccupied molecular orbital may take part in the toxic action betweeen the chemicals and macromolecules, and ENLUMO has the most important effect on these chemicals＇ toxicity.

Relationships between moment magnitude and fault parameters:theoretical and semi-empirical relationships

Fault parameters are important in earthquake hazard analysis.In this paper,theoretical relationships between moment magnitude and fault parameters including subsurface rupture length,downdip rupture width,rupture area,and average slip over the fault surface are deduced based on seismological theory.These theoretical relationships are further simplified by applying similarity conditions and an unique form is established.Then,combining the simplified theoretical relationships between moment magnitude and fault parameters with seismic source data selected in this study,a practical semi-empirical relationship is established.The seismic source data selected is also to used to derive empirical relationships between moment magnitude and fault parameters by the ordinary least square regression method.Comparisons between semi-empirical relationships and empirical relationships show that the former depict distribution trends of data better than the latter.It is also observed that downdip rupture widths of strike slip faults are saturated when moment magnitude is more than 7.0,but downdip rupture widths of dip slip faults are not saturated in the molnent magnitude rangcs of this study.

Towards a semi-empirical model of the sea ice thickness based on hyperspectral remote sensing in the Bohai Sea

Sea ice thickness is one of the most important input parameters for the prevention and mitigation of sea ice disasters and the prediction of local sea environments and climates. Estimating the sea ice thickness is currently the most important issue in the study of sea ice remote sensing. With the Bohai Sea as the study area, a semiempirical model of the sea ice thickness（SEMSIT） that can be used to estimate the thickness of first-year ice based on existing water depth estimation models and hyperspectral remote sensing data according to an optical radiative transfer process in sea ice is proposed. In the model, the absorption and scattering properties of sea ice in different bands（spectral dimension information） are utilized. An integrated attenuation coefficient at the pixel level is estimated using the height of the reflectance peak at 1 088 nm. In addition, the surface reflectance of sea ice at the pixel level is estimated using the 1 550–1 750 nm band reflectance. The model is used to estimate the sea ice thickness with Hyperion images. The first validation results suggest that the proposed model and parameterization scheme can effectively reduce the estimation error associated with the sea ice thickness that is caused by temporal and spatial heterogeneities in the integrated attenuation coefficient and sea ice surface. A practical semi-empirical model and parameterization scheme that may be feasible for the sea ice thickness estimation using hyperspectral remote sensing data are potentially provided.

Comparative Study among Different Semi-Empirical Models for Soil Salinity Prediction in an Arid Environment Using OLI Landsat-8 Data

Salt-affected soils, caused by natural or human activities, are a common environmental hazard in semi-arid and arid landscapes. Excess salts in soils affect plant growth and production, soil and water quality and, therefore, increase soil erosion and land degradation. This research investigates the performance of five different semi-empirical predictive models for soil salinity spatial distribution mapping in arid environment using OLI sensor image data. This is the first attempt to test remote sensing based semi-empirical salinity predictive models in this area: the Kingdom of Bahrain. To achieve our objectives, OLI data were standardized from the atmosphere interferences, the sensor radiometric drift, and the topographic and geometric distortions. Then, the five semi-empirical predictive models based on the Normalized Difference Salinity Index (NDSI), the Salinity Index-ASTER (SI-ASTER), the Salinity Index-1 (SI-1), the Soil Salinity and Sodicity Index-1 and Index-2 (SSSI-1 and SSSI-2), developed for slight and moderate salinity in agricultural land, were implemented and applied to OLI image data. For validation purposes, a fieldwork was organized and different important spots-locations representing different salinity levels were visited, photographed, and localized using an accurate GPS (σ ≤ ±30 cm). Based on this a priori knowledge of the soil salinity, six validation sites were selected to reflect non-saline, low, moderate, high and extreme salinity classes, descriptive statistics extracted from polygons and/or transects over these sites were used. The obtained results showed that the models based on NDSI, SI-1 and SI-ASTER all failed to detect salinity bounds for both extreme salinity (Sabkhah) and non-saline conditions. In Fact, NDSI and SI-ASTER gave respectively only 35% dS/m and 25% dS/m in extreme salinity validation site, while SI-1 and SI-ASTER indicated 38% dS/m and 39% dS/m in non-saline validation site. Therefore, these three models were deemed inadequate for the study site. However, both SSSI-1 and SSSI-2 allowed a detection of the previous salinity bounds and furthermore described similarly and correctly the urban-vegetation areas and the open-land areas. Their predicted EC is around 10% dS/m for non-saline urban soil, about 25% dS/m for low salinity urban-vegetation soil, approximately 30% to 75% dS/m, respectively, for moderate to high salinity soils. SSSI-2 based semi-empirical salinity models was able to differentiate the high salinity versus extreme salinity in areas where both exist and was very accurate to highlight the pure salt where SSSI-1 has reach saturation for both salinity classes. In conclusion, reliable salinity map was produced using the model based on SSSI-2 and OLI sensor data that allows a better characterization of the soil salinity problem in an Arid Environment.

Semi-empiricial Likelihood Confidence Intervals for the Differences of Two Populations Based on Fractional Imputation

Suppose that there are two populations x and y with missing data on both of them, where x has a distribution function F（·） which is unknown and y has a distribution function Gθ（·） with a probability density function gθ（·） with known form depending on some unknown parameter θ. Fractional imputation is used to fill in missing data. The asymptotic distributions of the semi-empirical likelihood ration statistic are obtained under some mild conditions. Then, empirical likelihood confidence intervals on the differences of x and y are constructed.

A semi-empirical formula for evaluating residual strain of soils under earthquake loading

A new semi-empirical formula for evaluating the residual strain of soils under earthquake loading is presented in this paper based on the incremental method and the increment model proposed by the authors.When the incident loading is uniform,the results calculated by the new formula are nearly the same as those by the existing formula.For excitation of the random earthquake loading,the results calculated by the new formula are compared to the results obtained by dynamic triaxial tests.The dynamic triaxial tests had been performed considering different seismic waves,confining stresses, consolidation ratios,and types of cohesive soils.The comparison between the calculated and tested results indicate that the presented formula can efficiently and practically describe the time-dependent process of the soil residual strains under actual seismic loads.

An improved semi-empirical friction model for gas-liquid two-phase flow in horizontal and near horizontal pipes

Pressure drop and liquid hold-up are two very important fluid flow parameters in design and control of multiphase flow pipelines.Friction factors play an important role in the accurate calculation of pressure drop.Various empirical and semi-empirical closure relations exist in the literature to calculate the liquid-wall,gas-wall and interfacial friction in two-phase pipe flow.However most of them are empirical correlations found under special experimental conditions.In this paper by modification of a friction model available in the literature,an improved semiempirical model is proposed.The proposed model is incorporated in the two-fluid correlations under equilibrium conditions and solved.Pressure gradient and velocity profiles are validated against experimental data.Using the improved model,the pressure gradient deviation from experiments diminishes by about 3%;the no-slip condition at the interface is satisfied and the velocity profile is predicted in better agreement with the experimental data.

Adaptive semi-empirical model for non-contact atomic force microscopy

Non-contact atomic force microscope is a powerful tool to investigate the surface topography with atomic resolution.Here we propose a new approach to estimate the interaction between its tips and samples,which combines a semi-empirical model with density functional theory(DFT)calculations.The generated frequency shift images are consistent with the experiment for mapping organic molecules using CuCO,Cu,CuCl,and CuO_(x)tips.This approach achieves accuracy close to DFT calculation with much lower computational cost.

Simulation of a DI Diesel Engine Performance Fuelled on Biodiesel Using a Semi-Empirical 0D Model

Diesel engines have proven over the years important in terms of efficiency and fuel consumption to power generation ratio. Many research works show the potential of biodiesel as a substitute for conventional gasoil. Mainly, previous and recent researches have focused on experimental investigation of diesel engine performance fuelled by biodiesel. Researches on the mathematical description of diesel engine process running on biodiesel are scarce, and mostly about chemical and thermodynamic description of the combustion process of biodiesel rather than performance studies. This work describes a numerical investigation on the performance analysis of a diesel engine fuelled by palm oil biodiesel. The numerical investigation was made using a semi empirical 0D model based on Wiebe’s and Watson’s model which was implemented via the open access numerical calculation software Scilab. The model was validated first by comparing with experimental pressure and performance data of a one cylinder engine at rated speed and secondly by comparing with a six cylinders engine performance data at various crankshaft rotational speeds. Simulations were then made to analyze the engine performance when running on biodiesel. The calculations were made at constant combustion duration and constant coefficient of excess air. Results showed that the model matches the overall experimental data, such as the power output and peak cylinder pressure. The ignition delay was somehow underestimated by the model for the first experiment, which caused a slight gap on in cylinder pressure curve, whereas it predicted the average ignition delay fairly well for the second set of validation. The simulations of engine performance when running on biodiesel confirmed results obtained in previous experimental researches on biodiesel. The model will be further investigated for engine control when shifting to biodiesel fuel.

The study of a neutron spectrum unfolding method based on particle swarm optimization combined with maximum likelihood expectation maximization

The neutron spectrum unfolding by Bonner sphere spectrometer(BSS) is considered a complex multidimensional model,which requires complex mathematical methods to solve the first kind of Fredholm integral equation. In order to solve the problem of the maximum likelihood expectation maximization(MLEM) algorithm which is easy to suffer the pitfalls of local optima and the particle swarm optimization(PSO) algorithm which is easy to get unreasonable flight direction and step length of particles, which leads to the invalid iteration and affect efficiency and accuracy, an improved PSO-MLEM algorithm, combined of PSO and MLEM algorithm, is proposed for neutron spectrum unfolding. The dynamic acceleration factor is used to balance the ability of global and local search, and improves the convergence speed and accuracy of the algorithm. Firstly, the Monte Carlo method was used to simulated the BSS to obtain the response function and count rates of BSS. In the simulation of count rate, four reference spectra from the IAEA Technical Report Series No. 403 were used as input parameters of the Monte Carlo method. The PSO-MLEM algorithm was used to unfold the neutron spectrum of the simulated data and was verified by the difference of the unfolded spectrum to the reference spectrum. Finally, the 252Cf neutron source was measured by BSS, and the PSO-MLEM algorithm was used to unfold the experimental neutron spectrum.Compared with maximum entropy deconvolution(MAXED), PSO and MLEM algorithm, the PSO-MLEM algorithm has fewer parameters and automatically adjusts the dynamic acceleration factor to solve the problem of local optima. The convergence speed of the PSO-MLEM algorithm is 1.4 times and 3.1 times that of the MLEM and PSO algorithms. Compared with PSO, MLEM and MAXED, the correlation coefficients of PSO-MLEM algorithm are increased by 33.1%, 33.5% and 1.9%, and the relative mean errors are decreased by 98.2%, 97.8% and 67.4%.

Application of GC/EIMS in Combination with Semi-Empirical Calculations for Identification and Investigation of Some Volatile Components in Basil Essential Oil

Volatile components in the extracts of basil leaves (Ocimum basilicum L.) were identified by gas chromatography/mass spectrometry (GC/MS) with electron ionization (EI) mode. The major volatile components of basil under investigation are α-pinene, sabinene, β-pinene, d-limonene, eucalyptol, l-linalool and estragole. Electron ionization mass spectra of these compounds have been obtained and investigated. Furthermore, the semi-empirical MNDO [Modified Neglect of Diatomic Overlap] method was used to calculate the thermochemical data for the structural properties of these compounds.

A synthetic semi-empirical physical model of secondary electron yield of metals under E-beam irradiation

Calculations of secondary electron yield（SEY） by physical formula can hardly accord with experimental results precisely. Simplified descriptions of internal electron movements in the calculation and complex surface contamination states of real sample result in notable difference between simulations and experiments. In this paper, in order to calculate SEY of metal under complicated surface state accurately, we propose a synthetic semi-empirical physical model. The processes of excitation of internal secondary electron（SE） and movement toward surface can be simulated using this model.This model also takes into account the influences of incident angle and backscattering electrons as well as the surface gas contamination. In order to describe internal electronic states accurately, the penetration coefficient of incident electron is described as a function of material atom number. Directions of internal electrons are set to be uniform in each angle. The distribution of internal SEs is proposed by considering both the integration convergence and the cascade scattering process.In addition, according to the experiment data, relationship among desorption gas quantities, sample ultimate temperature and SEY is established. Comparing with experiment results, this synthetic semi-empirical physical model can describe the SEY of metal better than former formulas, especially in the aspect of surface contaminated states. The proposed synthetic semi-empirical physical model and presented results in this paper can be helpful for further studying SE emission, and offer an available method for estimating and taking advantage of SE emission accurately.

Salt-Affected Soil Mapping in an Arid Environment Using Semi-Empirical Model and Landsat-OLI Data

The aim of this research is to map the salt-affected soil in an arid environment using an advanced semi-empirical predictive model, Operational Land Imager (OLI) data, a digital elevation model (DEM), field soil sampling, and laboratory and statistical analyses. To achieve our objectives, the OLI data were atmospherically corrected, radiometric sensor drift was calibrated, and distortions of topography and geometry were corrected using a DEM. Then, the soil salinity map was derived using a semi-empirical predictive model based on the Soil Salinity and Sodicity Index-2 (SSSI-2). The vegetation cover map was extracted from the Transformed Difference Vegetation Index (TDVI). In addition, accurate DEM of 5-m pixels was used to derive topographic attributes (elevation and slope). Visual comparisons and statistical validation of the semi-empirical model using ground truth were undertaken in order to test its capability in an arid environment for moderate and strong salinity mapping. To accomplish this step, fieldwork was organized and 120 soil samples were collected with various degrees of salinity, including non-saline soil samples. Each one was automatically labeled using a digital camera and an accurate global positioning system (GPS) survey (σ ≤ ± 30 cm) connected in real time to the geographic information system (GIS) database. Subsequently, in the laboratory, the major exchangeable cations (Ca2+, Mg2+, Na+, K+, Cl- and SO42-), pH and the electrical conductivity (EC-Lab) were extracted from a saturated soil paste, as well as the sodium adsorption ratio (SAR) being calculated. The EC-Lab, which is generally accepted as the most effective method for soil salinity quantification was used for statistical analysis and validation purposes. The obtained results demonstrated a very good conformity between the derived soil salinity map from OLI data and the ground truth, highlighting six major salinity classes: Extreme, very high, high, moderate, low and non-saline. The laboratory chemical analyses corroborate these results. Furthermore, the semi-empirical predictive model provides good global results in comparison to the ground truth and laboratory analysis (EC-Lab), with correlation coefficient (R2) of 0.97, an index of agreement (D) of 0.84 (p < 0.05), and low overall root mean square error (RMSE) of 11%. Moreover, we found that topographic attributes have a substantial impact on the spatial distribution of salinity. The areas at a relatively high altitude and with hard bedrock are less susceptible to salinity, while areas at a low altitude and slope (≤2%) composed of Quaternary soil are prone to it. In these low areas, the water table is very close to the surface (≤1 m), and the absence of an adequate drainage network contributes significantly to waterlogging. Consequently, the intrusion and emergence of seawater at the surface, coupled with high temperature and high evaporation rates, contribute extensively to the soil salinity in the study area.

Improvement of Bare Soil Semi-Empirical Radar Backscattering Models (Oh and Dubois) with SAR Multi-Spectral Satellite Data (X-, C- and L-Bands)

The objective of this study is to improve the performance of semi-empirical radar backscatter models, which are mainly used in microwave remote sensing (Oh 1992, Oh 2004 and Dubois). The study is based on satellite and ground data collected on bare soil surfaces during the Multispectral Crop Monitoring experimental campaign of the CESBIO laboratory in 2010 over an agricultural region in southwestern France. The dataset covers a wide range of soil (viewing top soil moisture, surface roughness and texture) and satellite (at different frequencies: X-, C- and L-bands, and different incidence angles: 24.3° to 53.3°) configurations. The proposed methodology consists in identifying and correcting the residues of the models, depending on the surface properties (roughness, moisture, texture) and/or sensor characteristics (frequency, incidence angle). Finally, one model has been retained for each frequency domain. Results show that the enhancements of the models significantly increase the simulation performances. The coefficient of correlation increases of 23% in mean and the simulation errors (RMSE) are reduced to below 2 dB (at the X and C-bands) and to 1 dB at the L-band, compared to the initial models. At the X- and C-bands, the best performances of the modified models are provided by Dubois, whereas Oh 2004 is more suitable for the L-band (r is equal to 0.69, 0.65 and 0.85). Moreover, the modified models of Oh 1992 and 2004 and Dubois, developed in this study, offer a wider domain of validity than the initial formalism and increase the capabilities of retrieving the backscattering signal in view of applications of such approaches to stronglycontrasted agricultural surface states.

Weighted Maximum Likelihood Technique for Logistic Regression

In this paper, a weighted maximum likelihood technique (WMLT) for the logistic regression model is presented. This method depended on a weight function that is continuously adaptable using Mahalanobis distances for predictor variables. Under the model, the asymptotic consistency of the suggested estimator is demonstrated and properties of finite-sample are also investigated via simulation. In simulation studies and real data sets, it is observed that the newly proposed technique demonstrated the greatest performance among all estimators compared.