Sensitivity Analysis of Electromagnetic Scattering from Dielectric Targets with Polynomial Chaos Expansion and Method of Moments

In this paper,an adaptive polynomial chaos expansion method(PCE)based on the method of moments(MoM)is proposed to construct surrogate models for electromagnetic scattering and further sensitivity analysis.The MoM is applied to accurately solve the electric field integral equation(EFIE)of electromagnetic scattering from homogeneous dielectric targets.Within the bistatic radar cross section(RCS)as the research object,the adaptive PCE algorithm is devoted to selecting the appropriate order to construct the multivariate surrogate model.The corresponding sensitivity results are given by the further derivative operation,which is compared with those of the finite difference method(FDM).Several examples are provided to demonstrate the effectiveness of the proposed algorithm for sensitivity analysis of electromagnetic scattering from homogeneous dielectric targets.

Moments of inertia of triaxial nuclei in covariant density functional theory

The covariant density functional theory(CDFT)and five-dimensional collective Hamiltonian(5DCH)are used to analyze the experimental deformation parameters and moments of inertia(MoIs)of 12 triaxial nuclei as extracted by Allmond and Wood[J.M.Allmond and J.L.Wood,Phys.Lett.B 767,226(2017)].We find that the CDFT MoIs are generally smaller than the experimental values but exhibit qualitative consistency with the irrotational flow and experimental data for the relative MoIs,indicating that the intermediate axis exhibites the largest MoI.Additionally,it is found that the pairing interaction collapse could result in nuclei behaving as a rigid-body flow,as exhibited in the^(186-192)Os case.Furthermore,by incorporating enhanced CDFT MoIs(factor of f≈1.55)into the 5DCH,the experimental low-lying energy spectra and deformation parameters are reproduced successfully.Compared with both CDFT and the triaxial rotor model,the 5DCH demonstrates superior agreement with the experimental deformation parameters and low-lying energy spectra,respectively,emphasizing the importance of considering shape fluctuations.

The deterministic condition for the ground reaction force acting point on the combined knee valgus and tibial internal rotation moments in early phase of cutting maneuvers in female athletes

Background:Combined knee valgus and tibial internal rotation(VL+IR)moments have been shown to stress the anterior cruciate ligament(ACL)in several in vitro cadaveric studies.To utilize this knowledge for non-contact ACL injury prevention in sports,it is necessary to elucidate how the ground reaction force(GRF)acting point(center of pressure(CoP))in the stance foot produces combined knee VL+IR moments in risky maneuvers,such as cuttings.However,the effects of the GRF acting point on the development of the combined knee VL+IR moment in cutting are still unknown.Methods:We first established the deterministic mechanical condition that the CoP position relative to the tibial rotational axis differentiates the GRF vector’s directional probability for developing the combined knee VL+IR moment,and theoretically predicted that when the CoP is posterior to the tibial rotational axis,the GRF vector is more likely to produce the combined knee VL+IR moment than when the CoP is anterior to the tibial rotational axis.Then,we tested a stochastic aspect of our theory in a lab-controlled in vivo experiment.Fourteen females performed 60˚cutting under forefoot/rearfoot strike conditions(10 trials each).The positions of lower limb markers and GRF data were measured,and the knee moment due to GRF vector was calculated.The trials were divided into anterior-and posterior-CoP groups depending on the CoP position relative to the tibial rotational axis at each 10 ms interval from 0 to 100 ms after foot strike,and the occurrence rate of the combined knee VL+IR moment was compared between trial groups.Results:The posterior-CoP group showed significantly higher occurrence rates of the combined knee VL+IR moment(maximum of 82.8%)at every time point than those of the anterior-CoP trials,as theoretically predicted by the deterministic mechanical condition.Conclusion:The rearfoot strikes inducing the posterior CoP should be avoided to reduce the risk of non-contact ACL injury associated with the combined knee VL+IR stress.

Dynamic connectedness and network in the high moments of cryptocurrency,stock,and commodity markets

This study examines the connectedness in high-order moments between cryptocurrency,major stock(U.S.,U.K.,Eurozone,and Japan),and commodity(gold and oil)markets.Using intraday data from 2020 to 2022 and the time and frequency connectedness models of Diebold and Yilmaz(Int J Forecast 28(1):57–66,2012)and Barunik and Křehlik(J Financ Econom 16(2):271–296,2018),we investigate spillovers among the markets in realized volatility,the jump component of realized volatility,realized skewness,and realized kurtosis.These higher-order moments allow us to identify the unique characteristics of financial returns,such as asymmetry and fat tails,thereby capturing various market risks such as downside risk and tail risk.Our results show that the cryptocurrency,stock,and commodity markets are highly connected in terms of volatility and in the jump component of volatility,while their connectedness in skewness and kurtosis is smaller.Moreover,jump and volatility connectedness are more persistent than that of skewness and kurtosis connectedness.Our rolling-window analysis of the connectedness models shows that connectedness varies over time across all moments,and tends to increase during periods of high uncertainty.Finally,we show the potential of gold and oil as hedging and safe-haven investments for other markets given that they are the least connected to other markets across all moments and investment horizons.Our findings provide useful information for designing effective portfolio management and cryptocurrency regulations.

Unusual magnetic relaxation in a single-molecule magnet with toroidal magnetic moments

We report the synthesis and characterization of a single-molecule magnet composed of triangular clusters of dysprosium ions.The structural study shows that the symmetry changes from one polar point group(mm2)at room temperature to another polar point group(m)at low temperature.Magnetic studies and theory calculations illustrate that the vortex distribution of magnetic dipoles in the triangular dysprosium clusters forms a toroidal magnetic moment.Interestingly,the analysis of AC magnetic susceptibility reveals the coexistence of three distinct magnetic relaxation processes,corresponding to the Raman,Orbach,and QTM relaxation pathways,respectively.The sum of three modified Debye functions is successfully used to describe the multiple relaxation behavior.

Explaining Pomeranchuk Effect by Parity of Magnetic Moments of Leptons and Hadrons for Superconductivity in 3He and Graphene*

The mystery of superconductivity has intrigued scientists for 110 years now. The author in 2014 specifically predicted the superconductivity in carbon, sulfur and hydrogen compounds and generally predicted carbonaceous, hydrogeneous and sulfurous compounds in 2005 with reference to scattering to asymmetric orbital motions and associated spin and orbital exchanges between nuclei and electrons. The emphasis was in 2005 upon stronger electron and nuclear interactions and electron-phonon effects. But here the author develops more the un-gerade parity of the p and f orbitals and their contributions to the superconductivity at lower pressures and higher temperatures. On the bases of such, the role of parity from the origin and inflation of the Universe is noted and dark and bright energies and matters in the mature Universe are reasoned. Moreover, the superconductors are all reasoned by positive and negative nuclear magnetic moments (NMMs) with availability of un-gerade parities of p and f subshells and their orbitals. In addition to superconductivity, such positive and negative NMMs by Little Effect is presented for explaining Pomeranchuk Effect and thereby further explaining superconductivity and superfluidity of 3He. On the bases of successes of Little Effect via positive and negative NMMs, in particular negative NMMs of 3He, the superconductivity in twisted graphene is explained and also its recently discovered Pomeranchuk Effect.

Moments of Inertia, Magnetic Dipole Moments, and Electric Quadrupole Moments of the Lithium Isotopes

The single-particle Schrödinger fluid model is designed mainly to calculate the moments of inertia of the axially symmetric deformed nuclei by assuming that each nucleon in the nucleus is moving in a single-particle potential which is deformed with time t, through its parametric dependence on a classical shape variable α(t). Also, the Nilsson model is designed for the calculations of the single-particle energy levels, the magnetic dipole moments, and the electric quadrupole moments of axially symmetric deformed nuclei by assuming that all the nucleons are moving in the field of an anisotropic oscillator potential. On the other hand, the nuclear superfluidity model is designed for the calculations of the nuclear moments of inertia and the electric quadrupole moments of deformed nuclei which have no axes of symmetry by assuming that the nucleons are moving in a quadruple deformed potential. Furthermore, the cranked Nilsson model is designed for the calculations of the total nuclear energy and the quadrupole moments of deformed nuclei which have no axes of symmetry by modifying the Nilsson potential to include second and fourth order oscillations. Accordingly, to investigate whether the six p-shell isotopes 6Li, 7Li, 8Li, 9Li, 10Li, and 11Li have axes of symmetry or not, we applied the four mentioned models to each nucleus by calculating their moments of inertia, their magnetic dipole moments, and their electric quadrupole moments by varying the deformation parameter β and the non-axiality parameter γ in wide ranges of values for this reason. Hence for the assumption that these isotopes are deformed and have axes of symmetry, we applied the single-particle Schrödinger fluid model and the Nilsson model. On the other hand, for the assumption that these isotopes are deformed and have no axes of symmetry, we applied the cranked Nilsson model and the nuclear super fluidity model. As a result of our calculations, we can conclude that the nucleus 6Li may be assumed to be deformed and has an axis of symmetry.

The MW5.5 earthquake on August 6,2023,in Pingyuan,Shandong,China:A rupture on a buried fault

On August 6,2023,a magnitude MW5.5 earthquake struck Pingyuan County,Dezhou City,Shandong Province,China.This event was significant as no large earthquakes had been recorded in the region for over a century,and no active fault had been previously identified.This study collects 1309 P-wave arrival times and 866 S-wave arrival times from 74 seismic stations less than 200 km to the epicenter to constrain the spatial distribution of the mainshock and its 125 early aftershocks by the double difference earthquake relocation method,and selects 864 P-waveforms from 288 stations located within 800 km of the epicenter to constrain the focal mechanism solution of the mainshock through centroid moment tensor inversion.The relocation and the inversion indicate,the Pingyuan MW5.5 earthquake was caused by a rupture on a buried fault,likely an extensive segment of the Gaotang fault.This buried fault exhibited a dip of approximately 75°to the northwest,with a strike of 222°,similar to the Gaotang fault.The rupture initiated at the depth of 18.6 km and propagated upward and northeastward.However,the ground surface was not broken.The total duration of the rupture was~6.0 s,releasing the scalar moment of 2.5895×1017 N·m,equivalent to MW5.54.The moment rate reached the maximum only 1.4 seconds after the rupture initiation,and the 90%scalar moment was released in the first 4.6 s.In the first 1.4 seconds of the rupture process,the rupture velocity was estimated to be 2.6 km/s,slower than the local S-wave velocity.As the rupture neared its end,the rupture velocity decreased significantly.This study provides valuable insights into the seismic characteristics of the Pingyuan MW5.5 earthquake,shedding light on the previously unidentified buried fault responsible for the seismic activity in the region.Understanding the behavior of such faults is crucial for assessing seismic hazards and enhancing earthquake preparedness in the future.

Speed and surface steepness affect internal tibial loading during running

Background:Internal tibial loading is influenced by modifiable factors with implications for the risk of stress injury.Runners encounter varied surface steepness(gradients)when running outdoors and may adapt their speed according to the gradient.This study aimed to quantify tibial bending moments and stress at the anterior and posterior peripheries when running at different speeds on surfaces of different gradients.Methods:Twenty recreational runners ran on a treadmill at 3 different speeds(2.5 m/s,3.0 m/s,and 3.5 m/s)and gradients(level:0%;uphill:+5%,+10%,and+15%;downhill:-5%,-10%,and-15%).Force and marker data were collected synchronously throughout.Bending moments were estimated at the distal third centroid of the tibia about the medial-lateral axis by ensuring static equilibrium at each 1%of stance.Stress was derived from bending moments at the anterior and posterior peripheries by modeling the tibia as a hollow ellipse.Two-way repeated-measures analysis of variance were conducted using both functional and discrete statistical analyses.Results:There were significant main effects for running speed and gradient on peak bending moments and peak anterior and posterior stress.Higher running speeds resulted in greater tibial loading.Running uphill at+10%and+15%resulted in greater tibial loading than level running.Running downhill at-10%and-15%resulted in reduced tibial loading compared to level running.There was no difference between+5%or-5%and level running.Conclusion:Running at faster speeds and uphill on gradients≥+10%increased internal tibial loading,whereas slower running and downhill running on gradients≥-10%reduced internal loading.Adapting running speed according to the gradient could be a protective mechanism,providing runners with a strategy to minimize the risk of tibial stress injuries.

Tensile strength and failure behavior of rock-mortar interfaces: Direct and indirect measurements

The tensile strength at the rock-concrete interface is one of the crucial factors controlling the failure mechanisms of structures,such as concrete gravity dams.Despite the critical importance of the failure mechanism and tensile strength of rock-concrete interfaces,understanding of these factors remains very limited.This study investigated the tensile strength and fracturing processes at rock-mortar interfaces subjected to direct and indirect tensile loadings.Digital image correlation(DIC)and acoustic emission(AE)techniques were used to monitor the failure mechanisms of specimens subjected to direct tension and indirect loading(Brazilian tests).The results indicated that the direct tensile strength of the rock-mortar specimens was lower than their indirect tensile strength,with a direct/indirect tensile strength ratio of 65%.DIC strain field data and moment tensor inversions(MTI)of AE events indicated that a significant number of shear microcracks occurred in the specimens subjected to the Brazilian test.The presence of these shear microcracks,which require more energy to break,resulted in a higher tensile strength during the Brazilian tests.In contrast,microcracks were predominantly tensile in specimens subjected to direct tension,leading to a lower tensile strength.Spatiotemporal monitoring of the cracking processes in the rock-mortar interfaces revealed that they show AE precursors before failure under the Brazilian test,whereas they show a minimal number of AE events before failure under direct tension.Due to different microcracking mechanisms,specimens tested under Brazilian tests showed lower roughness with flatter fracture surfaces than those tested under direct tension with jagged and rough fracture surfaces.The results of this study shed light on better understanding the micromechanics of damage in the rock-concrete interfaces for a safer design of engineering structures.

Two-Stream Approximation to the Radiative Transfer Equation:A New Improvement and Comparative Accuracy with Existing Methods

Mathematical modeling of the interaction between solar radiation and the Earth's atmosphere is formalized by the radiative transfer equation(RTE), whose resolution calls for two-stream approximations among other methods. This paper proposes a new two-stream approximation of the RTE with the development of the phase function and the intensity into a third-order series of Legendre polynomials. This new approach, which adds one more term in the expression of the intensity and the phase function, allows in the conditions of a plane parallel atmosphere a new mathematical formulation of γparameters. It is then compared to the Eddington, Hemispheric Constant, Quadrature, Combined Delta Function and Modified Eddington, and second-order approximation methods with reference to the Discrete Ordinate(Disort) method(δ –128 streams), considered as the most precise. This work also determines the conversion function of the proposed New Method using the fundamental definition of two-stream approximation(F-TSA) developed in a previous work. Notably,New Method has generally better precision compared to the second-order approximation and Hemispheric Constant methods. Compared to the Quadrature and Eddington methods, New Method shows very good precision for wide domains of the zenith angle μ 0, but tends to deviate from the Disort method with the zenith angle, especially for high values of optical thickness. In spite of this divergence in reflectance for high values of optical thickness, very strong correlation with the Disort method(R ≈ 1) was obtained for most cases of optical thickness in this study. An analysis of the Legendre polynomial series for simple functions shows that the high precision is due to the fact that the approximated functions ameliorate the accuracy when the order of approximation increases, although it has been proven that there is a limit order depending on the function from which the precision is lost. This observation indicates that increasing the order of approximation of the phase function of the RTE leads to a better precision in flux calculations. However, this approach may be limited to a certain order that has not been studied in this paper.

Bayesian and Non-Bayesian Analysis for the Sine Generalized Linear Exponential Model under Progressively Censored Data

This article introduces a novel variant of the generalized linear exponential(GLE)distribution,known as the sine generalized linear exponential(SGLE)distribution.The SGLE distribution utilizes the sine transformation to enhance its capabilities.The updated distribution is very adaptable and may be efficiently used in the modeling of survival data and dependability issues.The suggested model incorporates a hazard rate function(HRF)that may display a rising,J-shaped,or bathtub form,depending on its unique characteristics.This model includes many well-known lifespan distributions as separate sub-models.The suggested model is accompanied with a range of statistical features.The model parameters are examined using the techniques of maximum likelihood and Bayesian estimation using progressively censored data.In order to evaluate the effectiveness of these techniques,we provide a set of simulated data for testing purposes.The relevance of the newly presented model is shown via two real-world dataset applications,highlighting its superiority over other respected similar models.

Construction of a complete set of similarity invariants using pseudo-Zernike moments

To resolve the completeness and independence of an invariant set derived by the traditional method, a systematic method for deriving a complete set of pseudo-Zernike moment similarity （translation, scale and rotation） invariants is described. First, the relationship between pseudo-Zernike moments of the original image and those of the image having the same shape but distinct orientation and scale is established. Based on this relationship, a complete set of similarity invariants can be expressed as a linear combination of the original pseudo-Zernike moments of the same order and lower order. The problem of image reconstruction from a finite set of the pseudo-Zernike moment invariants （PZMIs） is also investigated. Experimental results show that the proposed PZMIs have better performance than complex moment invariants.

Study on the Load Characteristics of Submerged Body Under Internal Solitary Waves on the Continental Shelf and Slope

Based on the high-quality observation data and the numerical simulation,the evolution characteristics of internal solitary waves(ISWs)and the load on the suspend submerged body are studied on the continental shelf and slope separately.The observed ISWs exhibit the first mode depression ISWs.The amplitudes of ISWs on the shelf and slope areas reach 50 m and 80 m,respectively.The upper layer velocity in the westward direction is about 0.8 m/s on the continental shelf and 0.9 m/s on the continental slope during the passing through of ISWs.The lower layer is dominated by the eastward compensating flow.In the vertical direction,the water in front of the wave flows downward,while the water behind the wave flows upward,and the maximum vertical velocity exceeds 0.2 m/s.Numerical simulation results show that the larger the amplitude of ISWs,the larger the load on the submerged body.The force on the submerged body by ISWs is dominated by the vertical force,and the corresponding maximum vertical forces on the continental shelf and slope are−25 kN and −27 kN.The submerged body is subjected to a large counterclockwise moment and the sudden change of the moment will also cause the submerged body to capsize.This paper not only gives a deeper understanding of the characteristics of ISWs from the deep continental slope to the shallow continental shelf,but also has a certain guiding value for the prediction of ISWs and for marine military activities.

Structural failure analysis with CMS-based ground motion selection using innovative cost function and weight factors

The selection and scaling of ground motion records is considered a primary and essential task in performing structural analysis and design.Conventional methods involve using ground motion models and a conditional spectrum to select ground motion records based on the target spectrum.This research demonstrates the influence of adopting different weighted factors for various period ranges during matching selected ground motions with the target hazard spectrum.The event data from the Next Generation Attenuation West 2(NGA-West 2)database is used as the basis for ground motion selection,and hazard de-aggregation is conducted to estimate the event parameters of interest,which are then used to construct the target intensity measure(IM).The target IMs are then used to select ground motion records with different weighted vector-valued objective functions.The weights are altered to account for the relative importance of IM in accordance with the structural analysis application of steel moment resisting frame(SMRF)buildings.Instead of an ordinary objective function for the matching spectrum,a novel model is introduced and compared with the conventional cost function.The results indicate that when applying the new cost function for ground motion selection,it places higher demands on structures compared to the conventional cost function.Moreover,submitting more weights to the first-mode period of structures increases engineering demand parameters.Findings demonstrate that weight factors allocated to different period ranges can successfully account for period elongation and higher mode effects.

Moment Redistribution Effect of the Continuous Glass Fiber Reinforced Polymer-Concrete Composite Slabs Based on Static Loading Experiment

This study aimed to investigate the moment redistribution in continuous glass fiber reinforced polymer(GFRP)-concrete composite slabs caused by concrete cracking and steel bar yielding in the negative bending moment zone.An experimental bending moment redistribution test was conducted on continuous GFRP-concrete composite slabs,and a calculation method based on the conjugate beam method was proposed.The composite slabs were formed by combining GFRP profiles with a concrete layer and supported on steel beams to create two-span continuous composite slab specimens.Two methods,epoxy resin bonding,and stud connection,were used to connect the composite slabs with the steel beams.The experimental findings showed that the specimen connected with epoxy resin exhibited two moments redistribution phenomena during the loading process:concrete cracking and steel bar yielding at the internal support.In contrast,the composite slab connected with steel beams by studs exhibited only one-moment redistribution phenomenon throughout the loading process.As the concrete at the internal support cracked,the bending moment decreased in the internal support section and increased in the midspan section.When the steel bars yielded,the bending moment further decreased in the internal support section and increased in the mid-span section.Since GFRP profiles do not experience cracking,there was no significant decrease in the bending moment of the mid-span section.All test specimens experienced compressive failure of concrete at the mid-span section.Calculation results showed good agreement between the calculated and experimental values of bending moments in the mid-span section and internal support section.The proposed model can effectively predict the moment redistribution behavior of continuous GFRP-concrete composite slabs.

Implementing and evaluating an automatic centroid moment tensor procedure for the Indonesia region and surrounding areas

The purpose of this research was to suggest an applicable procedure for computing the centroid moment tensor(CMT)automatically and in real time from earthquakes that occur in Indonesia and the surrounding areas.Gisola software was used to estimate the CMT solution by selecting the velocity model that best suited the local and regional geological conditions in Indonesia and the surrounding areas.The data used in this study were earthquakes with magnitudes of 5.4 to 8.0.High-quality,real-time broadband seismographic data were provided by the International Federation of Digital Seismograph Networks Web Services(FDSNWS)and the European Integrated Data Archive(EIDA)Federation in Indonesia and the surrounding areas.Furthermore,the inversion process and filter adjustment were carried out on the seismographic data to obtain good CMT solutions.The CMT solutions from Gisola provided good-quality solutions,in which all earthquake data had A-level quality(high quality,with good variant reduction).The Gisola CMT solution was justified with the Global CMT(GCMT)solution by using the Kagan angle value,with an average of approximately 11.2°.This result suggested that the CMT solution generated from Gisola was trustworthy and reliable.The Gisola CMT solution was typically available within approximately 15 minutes after an earthquake occurred.Once it met the quality requirement,it was automatically published on the internet.The catalog of local and regional earthquake records obtained through this technology holds great promise for improving the current understanding of regional seismic activity and ongoing tectonic processes.The accurate and real-time CMT solution generated by implementing the Gisola algorithm consisted of moment tensors and moment magnitudes,which provided invaluable insights into earthquakes occurring in Indonesia and the surrounding areas.

Seismic safety assessment with non-Gaussian random processes for train-bridge coupled systems

Extensive high-speed railway(HSR)network resembled the intricate vascular system of the human body,crisscrossing mainlands.Seismic events,known for their unpredictability,pose a significant threat to both trains and bridges,given the HSR’s extended operational duration.Therefore,ensuring the running safety of train-bridge coupled(TBC)system,primarily composed of simply supported beam bridges,is paramount.Traditional methods like the Monte Carlo method fall short in analyzing this intricate system efficiently.Instead,efficient algorithm like the new point estimate method combined with moment expansion approximation(NPEM-MEA)is applied to study random responses of numerical simulation TBC systems.Validation of the NPEM-MEA’s feasibility is conducted using the Monte Carlo method.Comparative analysis confirms the accuracy and efficiency of the method,with a recommended truncation order of four to six for the NPEM-MEA.Additionally,the influences of seismic magnitude and epicentral distance are discussed based on the random dynamic responses in the TBC system.This methodology not only facilitates seismic safety assessments for TBC systems but also contributes to standard-setting for these systems under earthquake conditions.

Towards fast focal mechanism inversion of shallow crustal earthquakes in the Chinese mainland

We have developed an automatic regional focal mechanism inversion system based on the Earthquake Rapid Report(ERR) system and the real-time three-component seismic waveform stream of 1 000 broadband seismic stations provided by the China Earthquake Networks Center(CENC). The system can rapidly provide a double couple solution and centroid depth within 5–15 min after receiving earthquake information from the ERR system.The data processing is triggered by earthquake information obtained from the ERR system. The system is capable of determining the focal mechanism of all shallow-depth earthquakes in the Chinese mainland with a magnitude of 5.5–6.5. It utilizes waveform data recorded by seismic stations located within 500 km from the epicenter,enabling the reporting of a focal mechanism solution within 5–15 min of an earthquake occurrence. Additionally,the system can assign a corresponding grade(A B C) to the focal mechanism solution. We processed a total of 301earthquakes that occurred from 2021 to June 2022, and after the quality control, 166 of them were selected.These selected solutions were manually checked, and 160 of them were compiled in a focal mechanism catalog.This catalog can be conveniently downloaded online via the Internet. The automatic focal mechanism solution of earthquakes in eastern China exhibits a good agreement with that provided by the Global Centroid Moment Tensor(GCMT), when available. The average Kagan angle between this catalog and GCMT is 22°, and the average difference in MWis 0.17. Furthermore, compared with GCMT, the minimum magnitude of our catalog has been reduced from approximately 5.0 to 4.0. The correlation between the centroid depth and crustal thickness in the Chinese mainland confirms the distribution of the centroid depth.

Parameters Estimation of a Bivariate Generalized Poisson Distribution with Applications to Metabolic Syndrome Data

Background: Bivariate count data are commonly encountered in medicine, biology, engineering, epidemiology and many other applications. The Poisson distribution has been the model of choice to analyze such data. In most cases mutual independence among the variables is assumed, however this fails to take into accounts the correlation between the outcomes of interests. A special bivariate form of the multivariate Lagrange family of distribution, names Generalized Bivariate Poisson Distribution, is considered in this paper. Objectives: We estimate the model parameters using the method of maximum likelihood and show that the model fits the count variables representing components of metabolic syndrome in spousal pairs. We use the likelihood local score to test the significance of the correlation between the counts. We also construct confidence interval on the ratio of the two correlated Poisson means. Methods: Based on a random sample of pairs of count data, we show that the score test of independence is locally most powerful. We also provide a formula for sample size estimation for given level of significance and given power. The confidence intervals on the ratio of correlated Poisson means are constructed using the delta method, the Fieller’s theorem, and the nonparametric bootstrap. We illustrate the methodologies on metabolic syndrome data collected from 4000 spousal pairs. Results: The bivariate Poisson model fitted the metabolic syndrome data quite satisfactorily. Moreover, the three methods of confidence interval estimation were almost identical, meaning that they have the same interval width.