The zero-energy variance principle can be exploited in variational quantum eigensolvers for solving general eigenstates but its capacity for obtaining a specified eigenstate,such as ground state,is limited as all eige...The zero-energy variance principle can be exploited in variational quantum eigensolvers for solving general eigenstates but its capacity for obtaining a specified eigenstate,such as ground state,is limited as all eigenstates are of zero energy variance.We propose a variance-based variational quantum eigensolver for solving the ground state by searching in an enlarged space of wavefunction and Hamiltonian.With a mutual variance-Hamiltonian optimization procedure,the Hamiltonian is iteratively updated to guild the state towards to the ground state of the target Hamiltonian by minimizing the energy variance in each iteration.We demonstrate the performance and properties of the algorithm with numeral simulations.Our work suggests an avenue for utilizing guided Hamiltonian in hybrid quantum-classical algorithms.展开更多
Doubled haploid(DH)plants have been widely used for breeding and biological research in crops.Pop ulus spp.have been used as model woody plant species for biological research.However,the induction of DH poplar plants ...Doubled haploid(DH)plants have been widely used for breeding and biological research in crops.Pop ulus spp.have been used as model woody plant species for biological research.However,the induction of DH poplar plants is onerous,and limited biological or breeding work has been carried out on DH individuals or populations.In this study,we provide an effective protocol for poplar haploid induction based on an anther culture method.A total of 96 whole DH plant lines were obtained using an F1hybrid of Populus simonii×P.nigra as a donor tree.The phenotypes of the DH population showed exceptionally high variance when compared to those of half-sib progeny of the donor tree.Each DH line displayed distinct features compared to those of the other DH lines or the donor tree.Additionally,some excellent homozygous lines have the potential to be model plants in genetic and breeding studies.展开更多
The inter-cycle correlation of fission source distributions(FSDs)in the Monte Carlo power iteration process results in variance underestimation of tallied physical quantities,especially in large local tallies.This stu...The inter-cycle correlation of fission source distributions(FSDs)in the Monte Carlo power iteration process results in variance underestimation of tallied physical quantities,especially in large local tallies.This study provides a mesh-free semiquantitative variance underestimation elimination method to obtain a credible confidence interval for the tallied results.This method comprises two procedures:Estimation and Elimination.The FSD inter-cycle correlation length is estimated in the Estimation procedure using the Sliced Wasserstein distance algorithm.The batch method was then used in the elimination procedure.The FSD inter-cycle correlation length was proved to be the optimum batch length to eliminate the variance underestimation problem.We exemplified this method using the OECD sphere array model and 3D PWR BEAVRS model.The results showed that the average variance underestimation ratios of local tallies declined from 37 to 87%to within±5%in these models.展开更多
[ Objective] The aim was to study variance type of capsule morphological characters in Platycodon grandiflorum population, and provide some theoretical basis for seeking to genetic markers which can differentiate diff...[ Objective] The aim was to study variance type of capsule morphological characters in Platycodon grandiflorum population, and provide some theoretical basis for seeking to genetic markers which can differentiate different P. grandiflorum and breeding new varieties. [ Method] According to shape morphological characters of capsule from the same population of perennial purple P. gandiflorum, seven types of distinct di- versity capsule were selected, variance analysis and multiple comparison on the length, diameter, length/diameter of the different types of capsule were carried out. [ Result] There is unicolor and bicolor, even trichrome, among main color was brown and purple. Capsule shape was main cone, furthermore, containing long roller type, spheroidicity and sphericity. [ Conclusion] P. gandiflorum capsule was divided into long form, short form and middle type from length/diameter size in perennial culture P. gandiflorum population.展开更多
This paper describes the application of the variance method for flux estimation over a mixed agricultural region in China. Eddy covariance and flux variance measurements were conducted in a near-surface layer over a n...This paper describes the application of the variance method for flux estimation over a mixed agricultural region in China. Eddy covariance and flux variance measurements were conducted in a near-surface layer over a non-uniform land surface in the central plain of China from 7 June to 20 July 2002. During this period, the mean canopy height was about 0.50 m. The study site consisted of grass (10% of area), beans (15%), corn (15%) and rice (60%). Under unstable conditions, the standard deviations of temperature and water vapor density (normalized by appropriate scaling parameters), observed by a single instrument, followed the Monin-Obukhov similarity theory. The similarity constants for heat (CT) and water vapor (Cq) were 1.09 and 1.49, respectively. In comparison with direct measurements using eddy covariance techniques, the flux variance method, on average, underestimated sensible heat flux by 21% and latent heat flux by 24%, which may be attributed to the fact that the observed slight deviations (20% or 30% at most) of the similarity "constants" may be within the expected range of variation of a single instrument from the generally-valid relations.展开更多
Peak ground acceleration(PGA) estimation is an important task in earthquake engineering practice.One of the most well-known models is the Boore-Joyner-Fumal formula,which estimates the PGA using the moment magnitude,t...Peak ground acceleration(PGA) estimation is an important task in earthquake engineering practice.One of the most well-known models is the Boore-Joyner-Fumal formula,which estimates the PGA using the moment magnitude,the site-to-fault distance and the site foundation properties.In the present study,the complexity for this formula and the homogeneity assumption for the prediction-error variance are investigated and an effi ciency-robustness balanced formula is proposed.For this purpose,a reduced-order Monte Carlo simulation algorithm for Bayesian model class selection is presented to obtain the most suitable predictive formula and prediction-error model for the seismic attenuation relationship.In this approach,each model class(a predictive formula with a prediction-error model) is evaluated according to its plausibility given the data.The one with the highest plausibility is robust since it possesses the optimal balance between the data fi tting capability and the sensitivity to noise.A database of strong ground motion records in the Tangshan region of China is obtained from the China Earthquake Data Center for the analysis.The optimal predictive formula is proposed based on this database.It is shown that the proposed formula with heterogeneous prediction-error variance is much simpler than the attenuation model suggested by Boore,Joyner and Fumal(1993).展开更多
This article studies the optimal proportional reinsurance and investment problem under a constant elasticity of variance (CEV) model. Assume that the insurer's surplus process follows a jump-diffusion process, the ...This article studies the optimal proportional reinsurance and investment problem under a constant elasticity of variance (CEV) model. Assume that the insurer's surplus process follows a jump-diffusion process, the insurer can purchase proportional reinsurance from the reinsurer via the variance principle and invest in a risk-free asset and a risky asset whose price is modeled by a CEV model. The diffusion term can explain the uncertainty associated with the surplus of the insurer or the additional small claims. The objective of the insurer is to maximize the expected exponential utility of terminal wealth. This optimization problem is studied in two cases depending on the diffusion term's explanation. In all cases, by using techniques of stochastic control theory, closed-form expressions for the value functions and optimal strategies are obtained.展开更多
In the present study, the authors investigated the relationship between the Arctic Oscillation (AO) and the high-frequency variability of daily sea level pressures in the Northern Hemisphere in winter (November throug...In the present study, the authors investigated the relationship between the Arctic Oscillation (AO) and the high-frequency variability of daily sea level pressures in the Northern Hemisphere in winter (November through March), using NCEP/NCAR reanalysis datasets for the time period of 1948/49-2000/01. High-frequency signals are defined as those with timescales shorter than three weeks and measured in terms of variance, for each winter for each grid. The correlations between monthly mean AO index and high-frequency variance are conducted. A predominant feature is that several regional centers with high correlation show up in the middle to high latitudes. Significant areas include mid- to high-latitude Asia centered at Siberia, northern Europe and the middle-latitude North Atlantic east of northern Africa. Their strong correlations can also be confirmed by the singular value decomposition analysis of covariance between mean SLP and high-frequency variance. This indicates that the relationship of AO with daily Sea Level Pressure (SLP) is confined to some specific regions in association with the inherent atmospheric dynamics. In middle-latitude Asia, there is a significant (at the 95% level) trend of variance of-2.26% (10 yr)-1. Another region that displays a strong trend is the northwestern Pacific with a significant rate of change of 0.80% (10 yr)-1. If the winter of 1948/49, an apparent outlier, is excluded, a steady linear trend of +1.51% (10 yr)-1 shows up in northern Europe. The variance probability density functions (PDFs) are found to change in association with different AO phases. The changes corresponding to high and low AO phases, however, are asymmetric in these regions. Some regions such as northern Europe display much stronger changes in high AO years, whereas some other regions such as Siberia show a stronger connection to low AO conditions. These features are supported by ECMWF reanalysis data. However, the dynamical mechanisms involved in the AO-high frequency SLP variance connection have not been well understood, and this needs further study.展开更多
Suppression effect in multiple regression analysis may be more common in research than what is currently recognized. We have reviewed several literatures of interest which treats the concept and types of suppressor va...Suppression effect in multiple regression analysis may be more common in research than what is currently recognized. We have reviewed several literatures of interest which treats the concept and types of suppressor variables. Also, we have highlighted systematic ways to identify suppression effect in multiple regressions using statistics such as: R2, sum of squares, regression weight and comparing zero-order correlations with Variance Inflation Factor (VIF) respectively. We also establish that suppression effect is a function of multicollinearity;however, a suppressor variable should only be allowed in a regression analysis if its VIF is less than five (5).展开更多
Bayes decision rule of variance components for one-way random effects model is derived and empirical Bayes (EB) decision rules are constructed by kernel estimation method. Under suitable conditions, it is shown that t...Bayes decision rule of variance components for one-way random effects model is derived and empirical Bayes (EB) decision rules are constructed by kernel estimation method. Under suitable conditions, it is shown that the proposed EB decision rules are asymptotically optimal with convergence rates near O(n-1/2). Finally, an example concerning the main result is given.展开更多
In order to enhance the robustness and contrast in the minimum variance(MV) beamformer, adaptive diagonal loading method was proposed. The conventional diagonal loading technique has already been used in the MV beamfo...In order to enhance the robustness and contrast in the minimum variance(MV) beamformer, adaptive diagonal loading method was proposed. The conventional diagonal loading technique has already been used in the MV beamformer, but has the drawback that its level is specified by predefined parameter and without consideration of input-data. To alleviate this problem, the level of diagonal loading was computed appropriately and automatically from the given data by shrinkage method in the proposed adaptive diagonal loaded beamformer. The performance of the proposed beamformer was tested on the simulated point target and cyst phantom was obtained using Field II. In the point target simulation, it is shown that the proposed method has higher lateral resolution than the conventional delay-and-sum beamformer and could be more robust in estimating the amplitude peak than the MV beamformer when acoustic velocity error exists. In the cyst phantom simulation, the proposed beamformer has shown that it achieves an improvement in contrast ratio and without distorting the edges of cyst.展开更多
The type 2 modified augmented design(MAD2) is an efficient unreplicated experimental design used for evaluating large numbers of lines in plant breeding and for assessing genetic variation in a population. Statistical...The type 2 modified augmented design(MAD2) is an efficient unreplicated experimental design used for evaluating large numbers of lines in plant breeding and for assessing genetic variation in a population. Statistical methods and data adjustment for soil heterogeneity have been previously described for this design. In the absence of replicated test genotypes in MAD2, their total variance cannot be partitioned into genetic and error components as required to estimate heritability and genetic correlation of quantitative traits, the two conventional genetic parameters used for breeding selection. We propose a method of estimating the error variance of unreplicated genotypes that uses replicated controls, and then of estimating the genetic parameters. Using the Delta method, we also derived formulas for estimating the sampling variances of the genetic parameters.Computer simulations indicated that the proposed method for estimating genetic parameters and their sampling variances was feasible and the reliability of the estimates was positively associated with the level of heritability of the trait. A case study of estimating the genetic parameters of three quantitative traits, iodine value, oil content, and linolenic acid content, in a biparental recombinant inbred line population of flax with 243 individuals, was conducted using our statistical models. A joint analysis of data over multiple years and sites was suggested for genetic parameter estimation. A pipeline module using SAS and Perl was developed to facilitate data analysis and appended to the previously developed MAD data analysis pipeline(http://probes.pw.usda.gov/bioinformatics_ tools/MADPipeline/index.html).展开更多
Started from the more general functional model and based on the work of Koch K R (1986) and QU Zi qiang (1989), marginal likelihood function of variance and covariance components is derived and is identical with the o...Started from the more general functional model and based on the work of Koch K R (1986) and QU Zi qiang (1989), marginal likelihood function of variance and covariance components is derived and is identical with the orthogonal complement likelihood function. Minimum norm quadratic unibiased estimator (MINQUE) is developed, which expands the formula by Rao C R (1973). It is proved that Helmert type estimation, MINQUE, BQUE(Best quadratic unibiased estimation) and maximum likelihood estimation are identical with one another. Besides, a universal formula for accuracy evalution is presented. Through these work, a universal theory of variance and covariance components is established.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.12005065)the Guangdong Basic and Applied Basic Research Fund(Grant No.2021A1515010317)。
文摘The zero-energy variance principle can be exploited in variational quantum eigensolvers for solving general eigenstates but its capacity for obtaining a specified eigenstate,such as ground state,is limited as all eigenstates are of zero energy variance.We propose a variance-based variational quantum eigensolver for solving the ground state by searching in an enlarged space of wavefunction and Hamiltonian.With a mutual variance-Hamiltonian optimization procedure,the Hamiltonian is iteratively updated to guild the state towards to the ground state of the target Hamiltonian by minimizing the energy variance in each iteration.We demonstrate the performance and properties of the algorithm with numeral simulations.Our work suggests an avenue for utilizing guided Hamiltonian in hybrid quantum-classical algorithms.
基金supported by the National Key R&D Program of China(2021YFD2200203)Heilongjiang Province Key R&D Program of China(GA21B010)+1 种基金Heilongjiang Touyan Innovation Team Program(Tree Genetics and Breeding Innovation Team)Heilongjiang Postdoctoral Financial Assistance(LBH-Z21097)。
文摘Doubled haploid(DH)plants have been widely used for breeding and biological research in crops.Pop ulus spp.have been used as model woody plant species for biological research.However,the induction of DH poplar plants is onerous,and limited biological or breeding work has been carried out on DH individuals or populations.In this study,we provide an effective protocol for poplar haploid induction based on an anther culture method.A total of 96 whole DH plant lines were obtained using an F1hybrid of Populus simonii×P.nigra as a donor tree.The phenotypes of the DH population showed exceptionally high variance when compared to those of half-sib progeny of the donor tree.Each DH line displayed distinct features compared to those of the other DH lines or the donor tree.Additionally,some excellent homozygous lines have the potential to be model plants in genetic and breeding studies.
基金supported by China Nuclear Power Engineering Co.,Ltd.Scientific Research Project(No.KY22104)the fellowship of China Postdoctoral Science Foundation(No.2022M721793).
文摘The inter-cycle correlation of fission source distributions(FSDs)in the Monte Carlo power iteration process results in variance underestimation of tallied physical quantities,especially in large local tallies.This study provides a mesh-free semiquantitative variance underestimation elimination method to obtain a credible confidence interval for the tallied results.This method comprises two procedures:Estimation and Elimination.The FSD inter-cycle correlation length is estimated in the Estimation procedure using the Sliced Wasserstein distance algorithm.The batch method was then used in the elimination procedure.The FSD inter-cycle correlation length was proved to be the optimum batch length to eliminate the variance underestimation problem.We exemplified this method using the OECD sphere array model and 3D PWR BEAVRS model.The results showed that the average variance underestimation ratios of local tallies declined from 37 to 87%to within±5%in these models.
基金supported by the Strategic Priority Research Program of the Chinese Academy of Sciences [grant number XDA20060500]the National Natural Science Foundation of China[grant numbers 41731173 and 42275035]+8 种基金the Natural Science Foundation of Guangdong ProvinceChina [grant number 2022A1515011967]the Science and Technology Program of GuangzhouChina [grant number 202002030492]the Open Fund Project of the Key Laboratory of Marine Environmental Information Technology,the Key Laboratory of Marine Science and Numerical Modeling,Ministry of Natural Resources of the People’s Republic of China [grant number 2020-YB-05]the MEL Visiting Fellowship [grant number MELRS2102]the Independent Research Project Program of the State Key Laboratory of Tropical Oceanography [grant number LTOZZ2005]the Key Special Project for the Introducing Talents Team of the Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou)[grant number GML2019ZD0306]the Innovation Academy of South China Sea Ecology and Environmental Engineering [grant number ISEE2018PY06]
文摘[ Objective] The aim was to study variance type of capsule morphological characters in Platycodon grandiflorum population, and provide some theoretical basis for seeking to genetic markers which can differentiate different P. grandiflorum and breeding new varieties. [ Method] According to shape morphological characters of capsule from the same population of perennial purple P. gandiflorum, seven types of distinct di- versity capsule were selected, variance analysis and multiple comparison on the length, diameter, length/diameter of the different types of capsule were carried out. [ Result] There is unicolor and bicolor, even trichrome, among main color was brown and purple. Capsule shape was main cone, furthermore, containing long roller type, spheroidicity and sphericity. [ Conclusion] P. gandiflorum capsule was divided into long form, short form and middle type from length/diameter size in perennial culture P. gandiflorum population.
文摘This paper describes the application of the variance method for flux estimation over a mixed agricultural region in China. Eddy covariance and flux variance measurements were conducted in a near-surface layer over a non-uniform land surface in the central plain of China from 7 June to 20 July 2002. During this period, the mean canopy height was about 0.50 m. The study site consisted of grass (10% of area), beans (15%), corn (15%) and rice (60%). Under unstable conditions, the standard deviations of temperature and water vapor density (normalized by appropriate scaling parameters), observed by a single instrument, followed the Monin-Obukhov similarity theory. The similarity constants for heat (CT) and water vapor (Cq) were 1.09 and 1.49, respectively. In comparison with direct measurements using eddy covariance techniques, the flux variance method, on average, underestimated sensible heat flux by 21% and latent heat flux by 24%, which may be attributed to the fact that the observed slight deviations (20% or 30% at most) of the similarity "constants" may be within the expected range of variation of a single instrument from the generally-valid relations.
基金Research Committee of University of Macao under Research Grant No.MYRG081(Y1-L2)-FST13-YKVthe Science and Technology Development Fund of the Macao SAR government under Grant No.012/2013/A1
文摘Peak ground acceleration(PGA) estimation is an important task in earthquake engineering practice.One of the most well-known models is the Boore-Joyner-Fumal formula,which estimates the PGA using the moment magnitude,the site-to-fault distance and the site foundation properties.In the present study,the complexity for this formula and the homogeneity assumption for the prediction-error variance are investigated and an effi ciency-robustness balanced formula is proposed.For this purpose,a reduced-order Monte Carlo simulation algorithm for Bayesian model class selection is presented to obtain the most suitable predictive formula and prediction-error model for the seismic attenuation relationship.In this approach,each model class(a predictive formula with a prediction-error model) is evaluated according to its plausibility given the data.The one with the highest plausibility is robust since it possesses the optimal balance between the data fi tting capability and the sensitivity to noise.A database of strong ground motion records in the Tangshan region of China is obtained from the China Earthquake Data Center for the analysis.The optimal predictive formula is proposed based on this database.It is shown that the proposed formula with heterogeneous prediction-error variance is much simpler than the attenuation model suggested by Boore,Joyner and Fumal(1993).
文摘This article studies the optimal proportional reinsurance and investment problem under a constant elasticity of variance (CEV) model. Assume that the insurer's surplus process follows a jump-diffusion process, the insurer can purchase proportional reinsurance from the reinsurer via the variance principle and invest in a risk-free asset and a risky asset whose price is modeled by a CEV model. The diffusion term can explain the uncertainty associated with the surplus of the insurer or the additional small claims. The objective of the insurer is to maximize the expected exponential utility of terminal wealth. This optimization problem is studied in two cases depending on the diffusion term's explanation. In all cases, by using techniques of stochastic control theory, closed-form expressions for the value functions and optimal strategies are obtained.
文摘In the present study, the authors investigated the relationship between the Arctic Oscillation (AO) and the high-frequency variability of daily sea level pressures in the Northern Hemisphere in winter (November through March), using NCEP/NCAR reanalysis datasets for the time period of 1948/49-2000/01. High-frequency signals are defined as those with timescales shorter than three weeks and measured in terms of variance, for each winter for each grid. The correlations between monthly mean AO index and high-frequency variance are conducted. A predominant feature is that several regional centers with high correlation show up in the middle to high latitudes. Significant areas include mid- to high-latitude Asia centered at Siberia, northern Europe and the middle-latitude North Atlantic east of northern Africa. Their strong correlations can also be confirmed by the singular value decomposition analysis of covariance between mean SLP and high-frequency variance. This indicates that the relationship of AO with daily Sea Level Pressure (SLP) is confined to some specific regions in association with the inherent atmospheric dynamics. In middle-latitude Asia, there is a significant (at the 95% level) trend of variance of-2.26% (10 yr)-1. Another region that displays a strong trend is the northwestern Pacific with a significant rate of change of 0.80% (10 yr)-1. If the winter of 1948/49, an apparent outlier, is excluded, a steady linear trend of +1.51% (10 yr)-1 shows up in northern Europe. The variance probability density functions (PDFs) are found to change in association with different AO phases. The changes corresponding to high and low AO phases, however, are asymmetric in these regions. Some regions such as northern Europe display much stronger changes in high AO years, whereas some other regions such as Siberia show a stronger connection to low AO conditions. These features are supported by ECMWF reanalysis data. However, the dynamical mechanisms involved in the AO-high frequency SLP variance connection have not been well understood, and this needs further study.
文摘Suppression effect in multiple regression analysis may be more common in research than what is currently recognized. We have reviewed several literatures of interest which treats the concept and types of suppressor variables. Also, we have highlighted systematic ways to identify suppression effect in multiple regressions using statistics such as: R2, sum of squares, regression weight and comparing zero-order correlations with Variance Inflation Factor (VIF) respectively. We also establish that suppression effect is a function of multicollinearity;however, a suppressor variable should only be allowed in a regression analysis if its VIF is less than five (5).
基金The project is partly supported by NSFC (19971085)the Doctoral Program Foundation of the Institute of High Education and the Special Foundation of Chinese Academy of Sciences.
文摘Bayes decision rule of variance components for one-way random effects model is derived and empirical Bayes (EB) decision rules are constructed by kernel estimation method. Under suitable conditions, it is shown that the proposed EB decision rules are asymptotically optimal with convergence rates near O(n-1/2). Finally, an example concerning the main result is given.
基金Project(2013GZX0147-3)supported by the Science and Technology Pillar Program of Sichuan Province,China
文摘In order to enhance the robustness and contrast in the minimum variance(MV) beamformer, adaptive diagonal loading method was proposed. The conventional diagonal loading technique has already been used in the MV beamformer, but has the drawback that its level is specified by predefined parameter and without consideration of input-data. To alleviate this problem, the level of diagonal loading was computed appropriately and automatically from the given data by shrinkage method in the proposed adaptive diagonal loaded beamformer. The performance of the proposed beamformer was tested on the simulated point target and cyst phantom was obtained using Field II. In the point target simulation, it is shown that the proposed method has higher lateral resolution than the conventional delay-and-sum beamformer and could be more robust in estimating the amplitude peak than the MV beamformer when acoustic velocity error exists. In the cyst phantom simulation, the proposed beamformer has shown that it achieves an improvement in contrast ratio and without distorting the edges of cyst.
基金partly supported by an A-base project funded by Agriculture and Agri-Food Canadathe TUFGEN project funded by Genome Canada and other stakeholdersfunds from the Western Grains Research Foundation
文摘The type 2 modified augmented design(MAD2) is an efficient unreplicated experimental design used for evaluating large numbers of lines in plant breeding and for assessing genetic variation in a population. Statistical methods and data adjustment for soil heterogeneity have been previously described for this design. In the absence of replicated test genotypes in MAD2, their total variance cannot be partitioned into genetic and error components as required to estimate heritability and genetic correlation of quantitative traits, the two conventional genetic parameters used for breeding selection. We propose a method of estimating the error variance of unreplicated genotypes that uses replicated controls, and then of estimating the genetic parameters. Using the Delta method, we also derived formulas for estimating the sampling variances of the genetic parameters.Computer simulations indicated that the proposed method for estimating genetic parameters and their sampling variances was feasible and the reliability of the estimates was positively associated with the level of heritability of the trait. A case study of estimating the genetic parameters of three quantitative traits, iodine value, oil content, and linolenic acid content, in a biparental recombinant inbred line population of flax with 243 individuals, was conducted using our statistical models. A joint analysis of data over multiple years and sites was suggested for genetic parameter estimation. A pipeline module using SAS and Perl was developed to facilitate data analysis and appended to the previously developed MAD data analysis pipeline(http://probes.pw.usda.gov/bioinformatics_ tools/MADPipeline/index.html).
文摘Started from the more general functional model and based on the work of Koch K R (1986) and QU Zi qiang (1989), marginal likelihood function of variance and covariance components is derived and is identical with the orthogonal complement likelihood function. Minimum norm quadratic unibiased estimator (MINQUE) is developed, which expands the formula by Rao C R (1973). It is proved that Helmert type estimation, MINQUE, BQUE(Best quadratic unibiased estimation) and maximum likelihood estimation are identical with one another. Besides, a universal formula for accuracy evalution is presented. Through these work, a universal theory of variance and covariance components is established.