In this paper, we consider the risk assessment problem under multi-levels and multiple mixture subpopulations. Our result is the generalization of the results of [1-5].1 Finite Mixture Normal ModelsIn dose-response s...In this paper, we consider the risk assessment problem under multi-levels and multiple mixture subpopulations. Our result is the generalization of the results of [1-5].1 Finite Mixture Normal ModelsIn dose-response studies, a class of phenomena that frequently occur are that experimental subjects (e.g., mice) may have different responses like ’none, mild, severe’ after a toxicant experiment, or ’getting worse, no change, getting better’ after a medical treatment, etc. These phenomena have attracted the attention of many researchers in recent years. Finite展开更多
We present a new nonparametric predictive inference(NPI)method using a power-normal model for accelerated life testing(ALT).Combined with the accelerating link function and imprecise probability theory,the proposed me...We present a new nonparametric predictive inference(NPI)method using a power-normal model for accelerated life testing(ALT).Combined with the accelerating link function and imprecise probability theory,the proposed method is a feasible way to predict the life of the product using ALT failure data.To validate the method,we run a series of simulations and conduct accelerated life tests with real products.The NPI lower and upper survival functions show the robustness of our method for life prediction.This is a continuous research,and some progresses have been made by updating the link function between different stress levels.We also explain how to renew and apply our model.Moreover,discussions have been made about the performance.展开更多
A problem that frequently occurs in biological experiments with laboratory animals is that some subjects are less susceptible to the treatment group than others. Finite mixture models have traditionally been used to d...A problem that frequently occurs in biological experiments with laboratory animals is that some subjects are less susceptible to the treatment group than others. Finite mixture models have traditionally been used to describe the distribution of responses in treated subjects for such studies. In this paper, we first study the mixture normal model with multi-levels and multiple mixture sub-populations under each level, with particular attention being given to the model in which the proportions of susceptibility are related to dose levels, then we use EM-algorithm to find the maximum likelihood estimators of model parameters. Our results are generalizations of the existing results. Finally, we illustrate realistic significance of the above extension based on a set of real dose-response data.展开更多
Extensive studies based on partition curve of gravity separation have been investigated. All created models are merely used to simulate density distribution at the same size fraction. However, they cannot be used to p...Extensive studies based on partition curve of gravity separation have been investigated. All created models are merely used to simulate density distribution at the same size fraction. However, they cannot be used to predictive distribution of materials depending on compound feature of density and size. According to this situation, an improved model of partition curve based on accumulation normal distribution, which was distinguished from conventional model of accumulation normal distribution for partition curve, was proposed in this paper. It could simulate density distribution at different size fractions by using the density-size compound index and conflating the partition curves at different size fractions as one partition curve. The feasibility of three compound indexes, including mass index, settlement index and transformation index, were investigated. Specific forms of the improved model were also proposed. It is found that transformation index leads to the best fitting results, while the fitting error is only 1.75 according to the fitting partition curve.展开更多
Joint location and scale models of the skew-normal distribution provide useful ex- tension for joint mean and variance models of the normal distribution when the data set under consideration involves asymmetric outcom...Joint location and scale models of the skew-normal distribution provide useful ex- tension for joint mean and variance models of the normal distribution when the data set under consideration involves asymmetric outcomes. This paper focuses on the maximum likelihood estimation of joint location and scale models of the skew-normal distribution. The proposed procedure can simultaneously estimate parameters in the location model and the scale model. Simulation studies and a real example are used to illustrate the proposed methodologies.展开更多
The plate flexure and normal faulting characteristics along the Tonga, Japan, Izu-Bonin and Mariana Trenches are investigated by combining observations and modeling of elastoplastic deformation of the subducting plate...The plate flexure and normal faulting characteristics along the Tonga, Japan, Izu-Bonin and Mariana Trenches are investigated by combining observations and modeling of elastoplastic deformation of the subducting plate. The observed average trench relief is found to be the smallest at the Japan Trench(3 km) and the largest at the Mariana Trench(4.9 km), and the average fault throw is the smallest at the Japan Trench(113 m) and the largest at the Tonga Trench(284 m). A subducting plate is modeled to bend and generate normal faults subjected to three types of tectonic loading at the trench axis: vertical loading, bending moment, and horizontal tensional force. It is inverted for the solutions of tectonic loading that best fit the observed plate flexure and normal faulting characteristics of the four trenches. The results reveal that a horizontal tensional force(HTF) for the Japan Trench is 33%, 50% and 60% smaller than those of the Mariana, Tonga and Izu-Bonin Trenches, respectively. The normal faults are modeled to penetrate to a maximum depth of 29, 23, 32 and 32 km below the sea floor for the Tonga,Japan, Izu-Bonin and Mariana Trenches, respectively, which is consistent with the depths of relocated normal faulting earthquakes in the Japan and Izu-Bonin Trenches. Moreover, it is argued that the calculated horizontal tensional force is generally positively correlated with the observed mean fault throw, while the integrated area of the reduction in the effective elastic thickness is correlated with the trench relief. These results imply that the HTF plays a key role in controlling the normal faulting pattern and that plate weakening can lead to significant increase in the trench relief.展开更多
We present an efficient three-dimensional coupled-mode model based on the Fourier synthesis technique. In principle, this model is a one-way model, and hence provides satisfactory accuracy for problems where the forwa...We present an efficient three-dimensional coupled-mode model based on the Fourier synthesis technique. In principle, this model is a one-way model, and hence provides satisfactory accuracy for problems where the forward scattering dominates. At the same time, this model provides an efficiency gain of an order of magnitude or more over two-way coupled-mode models. This model can be applied to three-dimensional range-dependent problems with a slowly varying bathymetry or internal waves. A numerical example of the latter is demonstrated in this work. Comparisons of both accuracy and efficiency between the present model and a benchmark model are also provided.展开更多
This paper based on the essay [1], studies in case that replicated observations are available in some experimental points., the parameters estimation of one dimensional linear errors-in-variables (EV) models. Asymptot...This paper based on the essay [1], studies in case that replicated observations are available in some experimental points., the parameters estimation of one dimensional linear errors-in-variables (EV) models. Asymptotic normality is established.展开更多
Conventional Akaike’s Information Criterion (AIC) for normal error models uses the maximum-likelihood estimator of error variance. Other estimators of error variance, however, can be employed for defining AIC for nor...Conventional Akaike’s Information Criterion (AIC) for normal error models uses the maximum-likelihood estimator of error variance. Other estimators of error variance, however, can be employed for defining AIC for normal error models. The maximization of the log-likelihood using an adjustable error variance in light of future data yields a revised version of AIC for normal error models. It also gives a new estimator of error variance, which will be called the “third variance”. If the model is described as a constant plus normal error, which is equivalent to fitting a normal distribution to one-dimensional data, the approximated value of the third variance is obtained by replacing (n-1) (n is the number of data) of the unbiased estimator of error variance with (n-4). The existence of the third variance is confirmed by a simple numerical simulation.展开更多
The distribution function is an important tool in the study of the stochastic variances. The normal distribution is very popular in the nature and our society. The idea of membership functions is the foundation of the...The distribution function is an important tool in the study of the stochastic variances. The normal distribution is very popular in the nature and our society. The idea of membership functions is the foundation of the fuzzy sets theory. While the fuzzy theory is widely used, the completely certain membership function that has no any fuzziness at all has been the bottleneck of the applications of this theory. Cloud models are the effective tools in transforming between qualitative concepts and their quantitative expressions. It can represent the fuzziness and randomness and their relations of uncertain concepts. Also cloud models can show the concept granularity in multi-scale spaces by the digital characteristic Entropy (En). The normal cloud model not only broadens the form conditions of the normal distribution but also makes the normal membership function be the expectation of the random membership degree. In this paper, the universality of the normal cloud model is proved, which is more superior and easier, and can fit the fuzziness and gentleness of human cognitive processing. It would be more applicable and universal in the representation of uncertain notions.展开更多
Seismic ground faulting is the greatest hazard for continuous buried pipelines.Over the years,researchers have attempted to understand pipeline behavior mostly via numerical modeling such as the finite element method....Seismic ground faulting is the greatest hazard for continuous buried pipelines.Over the years,researchers have attempted to understand pipeline behavior mostly via numerical modeling such as the finite element method.The lack of well-documented field case histories of pipeline failure from seismic ground faulting and the cost and complicated facilities needed for full-scale experimental simulation mean that a centrifuge-based method to determine the behavior of pipelines subjected to faulting is best to verify numerical approaches.This paper presents results from three centrifuge tests designed to investigate continuous buried steel pipeline behavior subjected to normal faulting.The experimental setup and procedure are described and the recorded axial and bending strains induced in a pipeline are presented and compared to those obtained via analytical methods.The influence of factors such as faulting offset,burial depth and pipe diameter on the axial and bending strains of pipes and on ground soil failure and pipeline deformation patterns are also investigated.Finally,the tensile rupture of a pipeline due to normal faulting is investigated.展开更多
Mesopelagic fish,the most important daily vertically migrating community in the oceans,are characterized by high lipid content which may obscure the interpretation of stable isotopes analysis.Demersal fish,which are i...Mesopelagic fish,the most important daily vertically migrating community in the oceans,are characterized by high lipid content which may obscure the interpretation of stable isotopes analysis.Demersal fish,which are important consumers in the food web dominated by mesopelagic fish,also have a high lipid content.Here we collected 127 fish samples from the South China Sea and evaluated the effect of lipid contents on△δ^(13)C of mesopelagic and demersal fish.In lipid-extracted mesopelagic fish,the C/N content ratio(<5.5)shows a clear correlation withΔδ^(13)C(the offset of bulk and lipid-extractedδ^(13)C values),especially in non-migratory and semi-migratory species;these values were less correlation in demersal fish.Based on our results,we suggest that mesopelagic and demersal fish in different regions of the South China Sea should be studied separately using appropriate correction models and less fit for the traditional model.Moreover,the C/N content ratio should be used cautiously for establishing the lipid normalization model,especially for the fish in migratory mesopelagic fish and demersal fish.Our results also reveal that mesopelagic fish across nearby regions could be analyzed together.The new models described here can be applied in future studies of mesopelagic and demersal fish in the South China Sea.展开更多
In this paper,several properties of one-way classification model with skew-normal random effects are obtained,such as moment generating function,density function and noncentral skew chi-square distribution,etc.Based o...In this paper,several properties of one-way classification model with skew-normal random effects are obtained,such as moment generating function,density function and noncentral skew chi-square distribution,etc.Based on the EM algorithm,we discuss the maximum likelihood(ML)estimation of unknown parameters.For testing problem of fixed effect,a parametric bootstrap(PB)approach is developed.Finally,some simulation results on the Type I error rates and powers of the PB approach are obtained,which show that the PB approach provides satisfactory performances on the Type I error rates and powers,even for small samples.For illustration,our main results are applied to a real data problem.展开更多
A tube model to simulate the normal black smoker system has been built, where the Darcy flow equation, the Ergun equation and the turbulent pipe flow equation are used respectively to describe the dynamic process of d...A tube model to simulate the normal black smoker system has been built, where the Darcy flow equation, the Ergun equation and the turbulent pipe flow equation are used respectively to describe the dynamic process of different key areas in the hydrothermal circulation system. At the same time, a convection - diffuse Equation for the temperature field is used for describe the exchange of thermal energy and the law of temperature variation. Combining the above facts and using efficient mathematical algorithms and programming with the MatLab programming language, the variation curves of temperature, pressure and mass flow rate are determined, while also the dynamic heat equilibrium and pressure equilibrium within the black smoker system are analyzed. On the basis of the model of the normal black smoker system, a megaplume formation model is further built. For instance, the hydrothermal venting plume on the Juan de Fuca Ridge has been simulated and the simulation results are fairly consistent with Baker＇ s imputed data in 1986. On the basis of the above productive simulation, a series of factors for megaplume formation and the effectiveness of the main parameters of the periodicity of the megaplume formation, temperature and the maximum mass flow rate are systematically discussed. Main conclusions are as follows： The normal black smoker system can evolve into a megaplum eruption. In fact, the passageway of the hydrothermal discharge is blocked by the hydrothermal sediments during the black smoker period, which leads to a hydrothermal fluid accumulation, rise of temperature and increase of buoyancy pressure under the seabed. After a period of 2 - 3 a, the megaplume hydrothermal eruption will occur when the increasing buoyancy pressure is high enough to crack the blockage （cap）. Meanwhile, the temperature of the heat source must exceed 500 ℃, while the highest temperature of the eruption fluid may be high up to 413 ℃, which is fairly consistent with the surveying data. If the temperature of the heat source were to be higher than 500 ℃, then the critical period for the megaplume formation could be obviously curtailed to be less than 1 a, while the critical temperature and the maximum mass flow rate are nearly invariable. As the permeability increases, the maximum mass flow rate increases gradually close to a steady value.展开更多
This letter presents a theoretical model of the normal （head-on） collisions between two soft spheres for predicting the experimental characteristic of the coefficient of restitution dependent on impact velocity. Aft...This letter presents a theoretical model of the normal （head-on） collisions between two soft spheres for predicting the experimental characteristic of the coefficient of restitution dependent on impact velocity. After the contact force law between the contacted spheres during a collision is phenomenologically formulated in terms of the compression or overlap displacement under considera- tion of an elastic-plastic loading and a plastic unloading subprocesses, the coefficient of restitution is gained by the dynamic equation of the contact process once an initial impact velocity is input. It is found that the theoretical predictions of the coefficient of restitution varying with the impact velocity are well in agreement with the existing experimental characteristics which are fitted by the explicit formula.展开更多
Generalized linear mixed models (GLMMs) are typically constructed by incorporating random effects into the linear predictor. The random effects are usually assumed to be normally distributed with mean zero and varianc...Generalized linear mixed models (GLMMs) are typically constructed by incorporating random effects into the linear predictor. The random effects are usually assumed to be normally distributed with mean zero and variance-covariance identity matrix. In this paper, we propose to release random effects to non-normal distributions and discuss how to model the mean and covariance structures in GLMMs simultaneously. Parameter estimation is solved by using Quasi-Monte Carlo (QMC) method through iterative Newton-Raphson (NR) algorithm very well in terms of accuracy and stabilization, which is demonstrated by real binary salamander mating data analysis and simulation studies.展开更多
In this paper, we propose a robust mixture regression model based on the skew scale mixtures of normal distributions (RMR-SSMN) which can accommodate asymmetric, heavy-tailed and contaminated data better. For the vari...In this paper, we propose a robust mixture regression model based on the skew scale mixtures of normal distributions (RMR-SSMN) which can accommodate asymmetric, heavy-tailed and contaminated data better. For the variable selection problem, the penalized likelihood approach with a new combined penalty function which balances the SCAD and l<sub>2</sub> penalty is proposed. The adjusted EM algorithm is presented to get parameter estimates of RMR-SSMN models at a faster convergence rate. As simulations show, our mixture models are more robust than general FMR models and the new combined penalty function outperforms SCAD for variable selection. Finally, the proposed methodology and algorithm are applied to a real data set and achieve reasonable results.展开更多
Normal mixture regression models are one of the most important statistical data analysis tools in a heterogeneous population. When the data set under consideration involves asymmetric outcomes, in the last two decades...Normal mixture regression models are one of the most important statistical data analysis tools in a heterogeneous population. When the data set under consideration involves asymmetric outcomes, in the last two decades, the skew normal distribution has been shown beneficial in dealing with asymmetric data in various theoretic and applied problems. In this paper, we propose and study a novel class of models： a skew-normal mixture of joint location, scale and skewness models to analyze the heteroscedastic skew-normal data coming from a heterogeneous population. The issues of maximum likelihood estimation are addressed. In particular, an Expectation-Maximization （EM） algorithm for estimating the model parameters is developed. Properties of the estimators of the regression coefficients are evaluated through Monte Carlo experiments. Results from the analysis of a real data set from the Body Mass Index （BMI） data are presented.展开更多
文摘In this paper, we consider the risk assessment problem under multi-levels and multiple mixture subpopulations. Our result is the generalization of the results of [1-5].1 Finite Mixture Normal ModelsIn dose-response studies, a class of phenomena that frequently occur are that experimental subjects (e.g., mice) may have different responses like ’none, mild, severe’ after a toxicant experiment, or ’getting worse, no change, getting better’ after a medical treatment, etc. These phenomena have attracted the attention of many researchers in recent years. Finite
基金the National Natural Science Foundation of China(No.11272082)the China Scholarship Council State Scholarship Fund(No.201506070017)
文摘We present a new nonparametric predictive inference(NPI)method using a power-normal model for accelerated life testing(ALT).Combined with the accelerating link function and imprecise probability theory,the proposed method is a feasible way to predict the life of the product using ALT failure data.To validate the method,we run a series of simulations and conduct accelerated life tests with real products.The NPI lower and upper survival functions show the robustness of our method for life prediction.This is a continuous research,and some progresses have been made by updating the link function between different stress levels.We also explain how to renew and apply our model.Moreover,discussions have been made about the performance.
基金Supported by the National Natural Science Foundation of China (No. 10571073)Specialized Research Fund for the Doctoral Program of Higher Education (No. 20070183023)Program for New Century Excellent Talents in University, Scientific Research Fund of Jilin University (No. 200810024)
文摘A problem that frequently occurs in biological experiments with laboratory animals is that some subjects are less susceptible to the treatment group than others. Finite mixture models have traditionally been used to describe the distribution of responses in treated subjects for such studies. In this paper, we first study the mixture normal model with multi-levels and multiple mixture sub-populations under each level, with particular attention being given to the model in which the proportions of susceptibility are related to dose levels, then we use EM-algorithm to find the maximum likelihood estimators of model parameters. Our results are generalizations of the existing results. Finally, we illustrate realistic significance of the above extension based on a set of real dose-response data.
基金the financial support from the National Natural Science Foundation of China (No. 51221462)
文摘Extensive studies based on partition curve of gravity separation have been investigated. All created models are merely used to simulate density distribution at the same size fraction. However, they cannot be used to predictive distribution of materials depending on compound feature of density and size. According to this situation, an improved model of partition curve based on accumulation normal distribution, which was distinguished from conventional model of accumulation normal distribution for partition curve, was proposed in this paper. It could simulate density distribution at different size fractions by using the density-size compound index and conflating the partition curves at different size fractions as one partition curve. The feasibility of three compound indexes, including mass index, settlement index and transformation index, were investigated. Specific forms of the improved model were also proposed. It is found that transformation index leads to the best fitting results, while the fitting error is only 1.75 according to the fitting partition curve.
基金Supported by the National Natural Science Foundation of China(11261025,11201412)the Natural Science Foundation of Yunnan Province(2011FB016)the Program for Middle-aged Backbone Teacher,Yunnan University
文摘Joint location and scale models of the skew-normal distribution provide useful ex- tension for joint mean and variance models of the normal distribution when the data set under consideration involves asymmetric outcomes. This paper focuses on the maximum likelihood estimation of joint location and scale models of the skew-normal distribution. The proposed procedure can simultaneously estimate parameters in the location model and the scale model. Simulation studies and a real example are used to illustrate the proposed methodologies.
基金The National Natural Science Foundation of China under contract Nos 41706056,91628301 and U1606401the Program of Chinese Academy of Sciences under contract Nos Y4SL021001,QYZDY-SSW-DQC005,YZ201325 and YZ201534+1 种基金the Natural Science Foundation of Guangdong Province of China under contract No.2017A030310066the China Ocean Mineral Resources R&D Association under contract No.DY135-S2-1-04
文摘The plate flexure and normal faulting characteristics along the Tonga, Japan, Izu-Bonin and Mariana Trenches are investigated by combining observations and modeling of elastoplastic deformation of the subducting plate. The observed average trench relief is found to be the smallest at the Japan Trench(3 km) and the largest at the Mariana Trench(4.9 km), and the average fault throw is the smallest at the Japan Trench(113 m) and the largest at the Tonga Trench(284 m). A subducting plate is modeled to bend and generate normal faults subjected to three types of tectonic loading at the trench axis: vertical loading, bending moment, and horizontal tensional force. It is inverted for the solutions of tectonic loading that best fit the observed plate flexure and normal faulting characteristics of the four trenches. The results reveal that a horizontal tensional force(HTF) for the Japan Trench is 33%, 50% and 60% smaller than those of the Mariana, Tonga and Izu-Bonin Trenches, respectively. The normal faults are modeled to penetrate to a maximum depth of 29, 23, 32 and 32 km below the sea floor for the Tonga,Japan, Izu-Bonin and Mariana Trenches, respectively, which is consistent with the depths of relocated normal faulting earthquakes in the Japan and Izu-Bonin Trenches. Moreover, it is argued that the calculated horizontal tensional force is generally positively correlated with the observed mean fault throw, while the integrated area of the reduction in the effective elastic thickness is correlated with the trench relief. These results imply that the HTF plays a key role in controlling the normal faulting pattern and that plate weakening can lead to significant increase in the trench relief.
基金Supported by the National Natural Science Foundation of China under Grant No 11774374the Natural Science Foundation of Shandong Province of China under Grant No ZR2016AL10
文摘We present an efficient three-dimensional coupled-mode model based on the Fourier synthesis technique. In principle, this model is a one-way model, and hence provides satisfactory accuracy for problems where the forward scattering dominates. At the same time, this model provides an efficiency gain of an order of magnitude or more over two-way coupled-mode models. This model can be applied to three-dimensional range-dependent problems with a slowly varying bathymetry or internal waves. A numerical example of the latter is demonstrated in this work. Comparisons of both accuracy and efficiency between the present model and a benchmark model are also provided.
基金the National Natural Science Foundation of China (Grant No. 19631040)
文摘This paper based on the essay [1], studies in case that replicated observations are available in some experimental points., the parameters estimation of one dimensional linear errors-in-variables (EV) models. Asymptotic normality is established.
文摘Conventional Akaike’s Information Criterion (AIC) for normal error models uses the maximum-likelihood estimator of error variance. Other estimators of error variance, however, can be employed for defining AIC for normal error models. The maximization of the log-likelihood using an adjustable error variance in light of future data yields a revised version of AIC for normal error models. It also gives a new estimator of error variance, which will be called the “third variance”. If the model is described as a constant plus normal error, which is equivalent to fitting a normal distribution to one-dimensional data, the approximated value of the third variance is obtained by replacing (n-1) (n is the number of data) of the unbiased estimator of error variance with (n-4). The existence of the third variance is confirmed by a simple numerical simulation.
文摘The distribution function is an important tool in the study of the stochastic variances. The normal distribution is very popular in the nature and our society. The idea of membership functions is the foundation of the fuzzy sets theory. While the fuzzy theory is widely used, the completely certain membership function that has no any fuzziness at all has been the bottleneck of the applications of this theory. Cloud models are the effective tools in transforming between qualitative concepts and their quantitative expressions. It can represent the fuzziness and randomness and their relations of uncertain concepts. Also cloud models can show the concept granularity in multi-scale spaces by the digital characteristic Entropy (En). The normal cloud model not only broadens the form conditions of the normal distribution but also makes the normal membership function be the expectation of the random membership degree. In this paper, the universality of the normal cloud model is proved, which is more superior and easier, and can fit the fuzziness and gentleness of human cognitive processing. It would be more applicable and universal in the representation of uncertain notions.
基金This work was conducted at the Physical Modeling and Centrifuge Laboratory of the Soil Mechanics and Foundation Engineering Department in the School of Civil Engineering,University of Tehran
文摘Seismic ground faulting is the greatest hazard for continuous buried pipelines.Over the years,researchers have attempted to understand pipeline behavior mostly via numerical modeling such as the finite element method.The lack of well-documented field case histories of pipeline failure from seismic ground faulting and the cost and complicated facilities needed for full-scale experimental simulation mean that a centrifuge-based method to determine the behavior of pipelines subjected to faulting is best to verify numerical approaches.This paper presents results from three centrifuge tests designed to investigate continuous buried steel pipeline behavior subjected to normal faulting.The experimental setup and procedure are described and the recorded axial and bending strains induced in a pipeline are presented and compared to those obtained via analytical methods.The influence of factors such as faulting offset,burial depth and pipe diameter on the axial and bending strains of pipes and on ground soil failure and pipeline deformation patterns are also investigated.Finally,the tensile rupture of a pipeline due to normal faulting is investigated.
基金the National Natural Science Foundation of China under contract Nos 42090043 and 41876074the National Basic Research Program(973 Program)of China under contract No.2014CB441502.
文摘Mesopelagic fish,the most important daily vertically migrating community in the oceans,are characterized by high lipid content which may obscure the interpretation of stable isotopes analysis.Demersal fish,which are important consumers in the food web dominated by mesopelagic fish,also have a high lipid content.Here we collected 127 fish samples from the South China Sea and evaluated the effect of lipid contents on△δ^(13)C of mesopelagic and demersal fish.In lipid-extracted mesopelagic fish,the C/N content ratio(<5.5)shows a clear correlation withΔδ^(13)C(the offset of bulk and lipid-extractedδ^(13)C values),especially in non-migratory and semi-migratory species;these values were less correlation in demersal fish.Based on our results,we suggest that mesopelagic and demersal fish in different regions of the South China Sea should be studied separately using appropriate correction models and less fit for the traditional model.Moreover,the C/N content ratio should be used cautiously for establishing the lipid normalization model,especially for the fish in migratory mesopelagic fish and demersal fish.Our results also reveal that mesopelagic fish across nearby regions could be analyzed together.The new models described here can be applied in future studies of mesopelagic and demersal fish in the South China Sea.
基金Supported by Zhejiang Provincial Philosophy and Social Science Planning Zhijiang Youth Project of China(Grant No.16ZJQN017YB)Ministry of Education of China,Humanities and Social Science Projects(Grant No.19YJA910006)+2 种基金Zhejiang Provincial Natural Science Foundation of China(Grant No.LY20A010019)Fundamental Research Funds for the Provincial Universities of Zhejiang(Grant No.GK199900299012-204)Zhejiang Provincial Statistical Science Research Base Project of China(Grant No.19TJJD08)
文摘In this paper,several properties of one-way classification model with skew-normal random effects are obtained,such as moment generating function,density function and noncentral skew chi-square distribution,etc.Based on the EM algorithm,we discuss the maximum likelihood(ML)estimation of unknown parameters.For testing problem of fixed effect,a parametric bootstrap(PB)approach is developed.Finally,some simulation results on the Type I error rates and powers of the PB approach are obtained,which show that the PB approach provides satisfactory performances on the Type I error rates and powers,even for small samples.For illustration,our main results are applied to a real data problem.
文摘A tube model to simulate the normal black smoker system has been built, where the Darcy flow equation, the Ergun equation and the turbulent pipe flow equation are used respectively to describe the dynamic process of different key areas in the hydrothermal circulation system. At the same time, a convection - diffuse Equation for the temperature field is used for describe the exchange of thermal energy and the law of temperature variation. Combining the above facts and using efficient mathematical algorithms and programming with the MatLab programming language, the variation curves of temperature, pressure and mass flow rate are determined, while also the dynamic heat equilibrium and pressure equilibrium within the black smoker system are analyzed. On the basis of the model of the normal black smoker system, a megaplume formation model is further built. For instance, the hydrothermal venting plume on the Juan de Fuca Ridge has been simulated and the simulation results are fairly consistent with Baker＇ s imputed data in 1986. On the basis of the above productive simulation, a series of factors for megaplume formation and the effectiveness of the main parameters of the periodicity of the megaplume formation, temperature and the maximum mass flow rate are systematically discussed. Main conclusions are as follows： The normal black smoker system can evolve into a megaplum eruption. In fact, the passageway of the hydrothermal discharge is blocked by the hydrothermal sediments during the black smoker period, which leads to a hydrothermal fluid accumulation, rise of temperature and increase of buoyancy pressure under the seabed. After a period of 2 - 3 a, the megaplume hydrothermal eruption will occur when the increasing buoyancy pressure is high enough to crack the blockage （cap）. Meanwhile, the temperature of the heat source must exceed 500 ℃, while the highest temperature of the eruption fluid may be high up to 413 ℃, which is fairly consistent with the surveying data. If the temperature of the heat source were to be higher than 500 ℃, then the critical period for the megaplume formation could be obviously curtailed to be less than 1 a, while the critical temperature and the maximum mass flow rate are nearly invariable. As the permeability increases, the maximum mass flow rate increases gradually close to a steady value.
基金supported by the Innovation Team Fund of the National Natural Science Foundation of China(11121202)
文摘This letter presents a theoretical model of the normal （head-on） collisions between two soft spheres for predicting the experimental characteristic of the coefficient of restitution dependent on impact velocity. After the contact force law between the contacted spheres during a collision is phenomenologically formulated in terms of the compression or overlap displacement under considera- tion of an elastic-plastic loading and a plastic unloading subprocesses, the coefficient of restitution is gained by the dynamic equation of the contact process once an initial impact velocity is input. It is found that the theoretical predictions of the coefficient of restitution varying with the impact velocity are well in agreement with the existing experimental characteristics which are fitted by the explicit formula.
文摘Generalized linear mixed models (GLMMs) are typically constructed by incorporating random effects into the linear predictor. The random effects are usually assumed to be normally distributed with mean zero and variance-covariance identity matrix. In this paper, we propose to release random effects to non-normal distributions and discuss how to model the mean and covariance structures in GLMMs simultaneously. Parameter estimation is solved by using Quasi-Monte Carlo (QMC) method through iterative Newton-Raphson (NR) algorithm very well in terms of accuracy and stabilization, which is demonstrated by real binary salamander mating data analysis and simulation studies.
文摘In this paper, we propose a robust mixture regression model based on the skew scale mixtures of normal distributions (RMR-SSMN) which can accommodate asymmetric, heavy-tailed and contaminated data better. For the variable selection problem, the penalized likelihood approach with a new combined penalty function which balances the SCAD and l<sub>2</sub> penalty is proposed. The adjusted EM algorithm is presented to get parameter estimates of RMR-SSMN models at a faster convergence rate. As simulations show, our mixture models are more robust than general FMR models and the new combined penalty function outperforms SCAD for variable selection. Finally, the proposed methodology and algorithm are applied to a real data set and achieve reasonable results.
基金Supported by the National Natural Science Foundation of China(11261025,11561075)the Natural Science Foundation of Yunnan Province(2016FB005)the Program for Middle-aged Backbone Teacher,Yunnan University
文摘Normal mixture regression models are one of the most important statistical data analysis tools in a heterogeneous population. When the data set under consideration involves asymmetric outcomes, in the last two decades, the skew normal distribution has been shown beneficial in dealing with asymmetric data in various theoretic and applied problems. In this paper, we propose and study a novel class of models： a skew-normal mixture of joint location, scale and skewness models to analyze the heteroscedastic skew-normal data coming from a heterogeneous population. The issues of maximum likelihood estimation are addressed. In particular, an Expectation-Maximization （EM） algorithm for estimating the model parameters is developed. Properties of the estimators of the regression coefficients are evaluated through Monte Carlo experiments. Results from the analysis of a real data set from the Body Mass Index （BMI） data are presented.
