Estimation of entropy and other functionals of a multivariate density

For a multivariate density f with respect to Lebesgue measure , the estimation of % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% Waa8qaaeaacaWGkbGaaiikaiaadAgacaGGPaGaamOzaiaadsgacqaH% 8oqBaSqabeqaniabgUIiYdaaaa!4404!\[\int {J(f)fd\mu } \], and in particular % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% Waa8qaaeaacaWGMbWaaWbaaSqabeaacaaIYaaaaOGaamizaiabeY7a% TbWcbeqab0Gaey4kIipaaaa!41E4!\[\int {f^2 d\mu } \] and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% Waa8qaaeaacaWGMbGaciiBaiaac+gacaGGNbGaamOzaiaadsgacqaH% 8oqBaSqabeqaniabgUIiYdaaaa!44AC!\[\int {f\log fd\mu } \], is studied. These two particular functionals are important in a number of contexts. Asymptotic bias and variance terms are obtained for the estimators % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% WaaybyaeqaleqabaGaey4jIKnaneaacaWGjbaaaOGaeyypa0Zaa8qa% aeaacaWGkbGaaiikamaawagabeWcbeqaaiabgEIizdqdbaGaamOzaa% aakiaacMcacaWGKbGaamOramaaBaaaleaacaWGobaabeaaaeqabeqd% cqGHRiI8aaaa!4994!\[\mathop I\limits^ \wedge = \int {J(\mathop f\limits^ \wedge )dF_N } \] and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% WaaybyaeqaleqabaGaeSipIOdaneaacaWGjbaaaOGaeyypa0Zaa8qa% aeaacaWGkbGaaiikamaawagabeWcbeqaaiabgEIizdqdbaGaamOzaa% aakiaacMcadaGfGbqabSqabeaacqGHNis2a0qaaiaadAgaaaGccaWG% KbGaeqiVd0galeqabeqdcqGHRiI8aaaa!4C40!\[\mathop I\limits^ \sim = \int {J(\mathop f\limits^ \wedge )\mathop f\limits^ \wedge d\mu } \], where % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% WaaybyaeqaleqabaGaey4jIKnaneaacaWGMbaaaaaa!3E9C!\[{\mathop f\limits^ \wedge }\] is a kernel density estimate of f and F _n is the empirical distribution function based on the random sample X ₁ ,..., X _n from f. For the two functionalsmentioned above, a first order bias term for Î can be made zero by appropriate choices of non-unimodal kernels. Suggestions for the choice of bandwidth are given; for % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% WaaybyaeqaleqabaGaey4jIKnaneaacaWGjbaaaOGaeyypa0Zaa8qa% aeaadaGfGbqabSqabeaacqGHNis2a0qaaiaadAgaaaGccaWGKbGaam% OramaaBaaaleaacaWGobaabeaaaeqabeqdcqGHRiI8aaaa!476C!\[\mathop I\limits^ \wedge = \int {\mathop f\limits^ \wedge dF_N } \], a study of optimal bandwidth is possible.This research was supported by an NSERC Grant and a UBC Killam Research Fellowship.     

On sequential procedures for the point estimation of the mean of a normal population

The sequential procedures developed by Starr (1966, Ann. Math. Statist., 37, 1173–1185) for estimating the mean of a normal population are further analyzed. Asymptotic properties of the regret and first two moments of the stopping rules are studied and second-order approximations are derived.     

Lifetime prediction of ABS polymers based on thermoanalytical data

The service life of ABS polymer, stabilized by 2-(3,5-di-tert-butyl-4-hydroxyanilino)-4,6-bis(octylthio)-1,3,5-triazine and containing 50% of a modifying rubber component, was estimated from oxidative induction times measured by DSC in isothermal mode in the temperature interval 140–170°C. The lifetime of ABS powder at the actual temperature of drying was predicted by linear extrapolation according to Arrhenius. However, the extrapolated value was much longer than the real lifetime determined from the long-term oven aging tests at 70 and 90°C, simulating the industrial drying process. The effect of changes in the apparent activation energy of oxidation due to antioxidant consumption during polymer aging is discussed.     

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.     

Interval estimation of the critical value in a general linear model

Summary This paper concerns interval estimation of the critical value θ which satisfies under the general linear model,Y _i =μ(x _i )+ε _i (i=1,2,···), where for and the functional forms off _j ^′ s are known. From an asymptotic expansion it is shown that, under reasonable conditions, the limiting distribution of is normal. Thus in the large-sample case a confidence interval for θ can be obtained. Such a result is useful when one is interested in carrying out a retrospective analysis rather than designing the experiment (as in the Kiefer-Wolfowitz procedure). In Section 3 a sequential procedure is considered for confidence intervals with fixed width 2d. It is shown that, for a given stopping variableN, is also asymptotically normal asd→0. Thus the coverage probability converges to 1−α (preassigned) asd→0. An example of application in estimating the phase parameter in circadian rhythms is given for the purpose of illustration. Research partially supported by the NSF Grant DMS-8502346.     

Polyparametric modification equation for estimating thermodynamic properties of energetic nitro compounds

Calculation methods, based on hybrid density-functional theory with the basis sets of B3LYP/ 6-311+G (2d, p)//B3LYP/6-31G(d, p)and B3LYP/6-31+G(d)//B3LYP/6-31G(d, p), were applied to determine the thermodynamic characteristics of various energetic nitro compounds. A parametric modification equation and the least-squares approach were used to identify 21 of the energetic research compounds. The atomization energies of these 21 compounds have an average relative error of 0.21–0.25% of the experimental values. The enthalpy (H _f) and the Gibbs energy (G _f) of formation have mean absolute errors of 10.8–11.4 kJ/ mol (2.6–2.7 kcal/mol) and 10.0–10.3 kJ/mol (2.4 kcal/ mol), respectively. The enthalpy and the Gibbs energy of formation obtained exceed those in the literature obtained by semiempirical calculations. The calibrated least-squares parameters and parametric equations were used to predict H _f and G _f for the five newly developed energetic nitro compounds for further applications.Acknowledgements. The authors would like to thank the National Science Council of the Republic of China for financial support of this work under grant no. NSC-91-2113-M-014-003. The National Center for High-Performance Computing providing the computation facility is also acknowledged.     

病态分析体系有偏估计的研究

运用广义岭估计和Ｌｉｕｋｅｊｉａｎ提出的有偏估计，对病态分析体系进行了数值模拟和实际光度测定，结果表明，广义岭估计显优于最小二乘估计，Ｌｉｕｋｅｊｉａｎ法有功效，可和为解析病态分析体系的化学计量学方法。     

Parameter identification of thermophilic anaerobic degradation of valerate

The considered mathematical model of the decomposition of valerate presents three unknown kinetic parameters, two unknown stoichiometric coefficients, and three unknown initial concentrations for biomass. Applying a structural identifiability study, we concluded that it is necessary to perform simultaneous batch experiments with differenitial conditions for estimating these parameters. Four simultaneous batch experiments were conducted at 55°C, characterized by four different initial acetate concentrations. Product inhibition of valerate degradation by acetate was considered. Practical identification was done optimizing the sum of the multiple determination coefficients for all measured state viariables and for all experiments simultaneously. The estimated values of kinetic parameters and stoichiometric coefficients were characterized by the parameter correlation matrix, the confidence interval, and the student's t-test at 9% significance level with positive results except for the saturation constant, for which more eperiments for improving its identifiability should be conducted. In this article, we discussekinetic parameter estimation methods.     

Preservation of transform orders under a semiparametric model and its mixture

In this paper, sufficient conditions for preservation of several transform orders under a typical family of semiparametric distributions are made. The preservation properties are developed to compare mixture semiparametric distributions. Possible applications of the achieved results to compare scale family of distributions and also compare coherent systems with dependent but identical components and series and parallel systems with randomized number of components are provided.     

Multiplicity of Solutions of the Weighted p-Laplacian with Isolated Singularity and Diffusion Suppressed by Convection

In this paper, the authors study the multiplicity of solutions to the weighted p-Laplacian with isolated singularity and di?usion suppressed by convection ?div(|x|^α|?u|^p?2?u) + λ 1/|x|^β |?u|^p?2?u · x = |x|^γg(|x|) in B \ {0} subject to nonlinear Robin boundary value condition |x|^α|?u|^p?2?u · ?n = A ? ρu on ?B,where λ > 0, B ? R^N(N ≥ 2) is the unit ball centered at the origin, α > 0, p > 1, β ∈ R, γ > ?N, g ∈ C([0,1]) with g(0) > 0, A ∈ R, ρ > 0 and ?n is the unit outward normal. The same problem with di?usion promoted by convection, namely λ ≤ 0, has already been discussed by the last two authors (Song-Yin (2012)), where the existence, nonexistence and classi?cation of singularities for solutions are presented. Completely di?erent from [Song, H. J. and Yin, J. X., Removable isolated singularities of solutions to the weighted p-Laplacian with singular convection, Math. Meth. Appl. Sci., 35, 2012, 1089–1100], in the present case λ > 0, namely the di?usion is suppressed by the convection, non-singular solutions are not only existent but also may be in?nite which vary according only to the values of solutions at the isolated singular point. At the same time, the singular solutions may exist only if the di?usion dominates the convection.