Inference on Clustered Survival Data Using Imputed Frailties

This article proposes a new method for fitting frailty models to clustered survival data that is intermediate between the fully parametric and nonparametric maximum likelihood estimation approaches. A parametric form is assumed for the baseline hazard, but only for the purpose of imputing the unobserved frailties. The regression coefficients are then estimated by solving an estimating equation that is the average of the partial likelihood score with respect to the conditional distribution of frailties given the observed data. We prove consistency and asymptotic normality of the resulting estimators and give associated closedform estimators of their variance. The algorithm is easy to implement and reduces to the ordinary Cox partial likelihood approach when the frailties have a degenerate distribution. Simulations indicate high efficiency and robustness of the resulting estimates. We apply our new approach to a study with clustered survival data on asthma in children in east Boston.     

中子针孔成像点扩展函数模拟研究

应用Geant4系统对中子针孔成像的点扩展函数进行了模拟研究.建立了点扩展函数的数学模型,并将其分布图像与能量沉积数据的分布图像作了比较,可看出点扩展函数与能量沉积数据吻合得较好.研究结果表明,在一定偏移量范围内,拟合得到的点扩展函数更好,能得到较好的图像分辨率.通过对不同入射能量在一定偏移量入射下的点扩展函数讨论认为,入射能量增大虽然会增强图像的对比度,但会降低分辨率.     

偏振光在多层散射介质中传输的蒙特卡罗模拟研究

基于米氏散射理论,运用斯托克斯-穆勒形式,结合子午面法和拒绝法,提出了偏振光在多层散射介质中传输的蒙特卡罗模型,并首次分析了当菲涅耳公式失效时,多层介质界面处的全反射行为.同时用该模型研究了微粒直径和介质折射率对漫反射和漫透射的影响.结果表明:无论是折射率匹配介质还是非匹配介质,当粒子浓度一致时,漫反射率都随微粒直径的增加而增加,随介质折射率的增加而减小,而漫透射率则相反.研究结果为生物组织的偏振光散射研究提供了理论验证模型.     

基于有效Monte Carlo模拟的外汇期权组合非线性VaR模型

为了克服极小概率事件发生概率估计的困难,提出了把重要抽样技术发展到外汇期权组合非线性VaR模型中,估计出组合损失概率。为了进一步达到减少模拟估计误差目的,在重要抽样技术基础上使用分层抽样技术,进行更有效的Monte Carlo模拟。数值结果表明,重要抽样技术算法比常用Monte Carlo模拟法的计算效率更有效;而重要抽样技术和分层抽样技术相结合算法比重要抽样技术算法更有效地减少模拟所要估计的组合损失概率的方差,有着更高的计算效率。     

Monte Carlo Simulation of Hard Hyperspheres in Six,Seven and Eight Dimensions for Low to Moderate Densities

The pair correlation function of hard hyperspheres in six, seven and eight dimensions is obtained from Monte Carlo simulations. The value of the pair correlation function at contact is compared with the results from molecular dynamics calculations and a variety of theoretical approaches. Remarkably good agreement is found with the simple, closed-form equations of Y. Song, E. A. Mason and R. M. Stratt, J. Phys. Chem., 93:6916–6919 (1989). The Monte Carlo results for the equation of state are compared with the theoretical expressions of M. Baus and J. L. Colot, Phys. Rev. A, 36:3912 (1987), M. Luban and J. P. J. Michels, Phys. Rev A, 41:6796 (1990), and high order virial expansions. In addition, in seven dimensions, comparisons are made with the exact PY solution provided by M. Robles, M. L. de Haro and A. Santos, J. Chem. Phys., 120:9113 (2004). Very good agreement was observed between theory and computer simulation in all dimensions.     

TEXONO中微子实验屏蔽效果的Monte Carlo研究

TEXONO实验采用CsI (T1)晶体量能器测量反应堆低能中微子的能谱和通量,确定中微子ν_e的反常磁矩并为中微子振荡实验做技术准备.详细描述了屏蔽体的安排,本底的主要来源,基于GEANT3软件包做了Monte Carlo研究,得到屏蔽体对本底的屏蔽效果.     

结构可靠度分析中的最小方差抽样

重要抽样技术是结构可靠度分析活跃的发展方向之一。本文研究了以抽样方差最小为条件的重要抽样问题，给出了正态抽样分布函数有关参数的计算公式，分析和实际模拟表明，在以正态分布进行了抽样的前提下，应用优化的参数，可使抽样效率进一步提高。     

平稳随机过程模拟的一个快速Hartley变换方法

由平稳随机过程谱表示定理导出了平稳随机过程蒙特卡罗模拟的一个快速Hartley变换方法。按照该方法，样本过程可直接由一正弦和余弦函数和的级数公式计算产生。可以证明，当级数项数N足够大时，模拟的样本过程可准确地反映平稳随机过程规定的性质；当样本过程足够多时，其总体均值和总体自相关函数均趋于相应目标函数。方法有如下一些特点：模拟的样本过程由相互正交的两个过程迭加而成，样本过程具备各态历经性质，样本过程随着级数项数N趋于无穷而渐近呈正态分布，可直接利用快速Hartley变换来大大提高模拟的计算效率。     

稳态光谱技术提取生物组织结构和功能互补信息研究

稳态光谱技术正被越来越多地用于组织结构和功能互补信息提取。通过仿体实验论证稳态光谱技术在光学参数和固有荧光提取方面的优越性。仿体实验中,血红蛋白被用作吸收子,聚苯乙烯小球被用作散射子,Stilbene被用来产生荧光。实验总共制备两组仿体,它们具有相同的散射系数和不同的吸收系数,在每组仿体中进行4次荧光素滴定。实验结果表明,组织光学参数和固有荧光均能够被精确提取,从而证明了稳态光谱技术在提取具有互补特性的组织结构和功能信息方面的优越性。