多重二元响应Probit模型的渐近有效估计

针对多重二元响应Probit模型提出了两步估计方法,第一步由边际似然得到参数√n相合的估计,第二步通过一步迭代得到渐近有效估计,由于只需一步迭代,因此在利用模拟方法计算信息阵时,可以增加模拟的次数,从而减少模拟所产生的扰动对估计的影响．     

H广义线性模型中的参数估计的比较

本文研究H广义线性模型中未知参数的两种估计方法,一种是边际似然函数法,另一种是Lee和Nelder提出来的L-N法．对于一类具有两个随机效应的典型的Poisson-Gamma类模型,在一些正则性条件之下,我们已经证明了其中固定效应卢的L-N估计的强相合性及渐近正态性,并得到了其收敛于真值的速度．针对这类模型,本文进一步给出了其边际似然函数的解析表达式,并且通过Monte Carlo模拟,对模型中固定效应β的边际似然估计和L—N估计进行了比较,模拟表明L—N估计比边际似然估计在拟Poisson-Gamma模型中有着更加优良的表现,具有更高的精度。     

物种总数的一类Bayes区间估计的改进

通常所使用的由传统方法得到的物种总数的一类Bayes置信区间不是最短的,就此意义而言也不是最优的。本文得到改进后的Bayes区间估计,并在各种有代表性的自由度情形下通过对比指出它相对于传统的Bayes区间估计的优越性所在,由分析可见,改进后的区间估计更能满足生产实践的需要,得到的物种总数的区间估计更为精确。     

利用多重抽样估计稀有物种的种类数

利用矩估计和二个稳健估计方法（jackknife估计,bootstrap估计）来处理野外生态学工作者的调查数据,在假定已经发现一些稀有物种的情形下,通过统计推断得到那些未被发现的物种的种类数。利用本文所提出的方法调查水稻水稻田的昆虫群落和林地的在面植被群落的稀有种是十分有效的。     

一类神经传导方程的变网格有限元方法及数值分析

研究在神经传播过程中的一类非线性拟双曲方程的初边值问题,对二维情形应用常规变换,提出了两种变网格有限元格式,最后通过细致的分析和估计得到了最佳阶的H1和L2 模误差估计结果,并且第二种格式使时间精度提高一阶,最后对第一种格式作了数值实验,指明方法是高效可行的．     

相依模型参数估计的相对效率

本文考虑相依模型,对其未知参数向量给出了其最佳线性无偏估计相对于协方差改进估计的四种相对效率,同时,还给出了最小二乘估计相对于协方差改进估计的三种相对效率,在不同条件下,分别给出了相对效率的上界与下界。     

非线性拟双曲方程的交替方向变网格有限元方法

研究在神经传播过程中的一类非线性拟双曲方程的初边值问题,提出了一种三维交替方向变网格有限元格式,应用微分方程先验估计的理论和技巧,得到了最佳阶的L^2模误差估计结果。并作了数值实验。指明方法是高效可行的．     

关于广义Potthoff—Roy估计

本文考察了生长曲线模型的定义形式,并因此建立了相应的广义Potthoff-Roy估计,在最小范数准则下,给出了估计的最佳选择并且讨论了协变量以及改进估计的方法,尤其当设计阵病态时,给出了两类新的岭型Potthoff-Roy估计。     

一类神经传播方程的特征差分方法和最佳阶l^2误差估计

给出矩形域上一类神经传播方程的特征差分,利用沿特征线方向构造差分逼近格式的方法和技巧．对给定的模型进行离散数值逼近和数值分析．特别是在沿特征线方向构造离散差分格式的过程中,可能会出现离散点在定义域之外的问题．本文提供了一个新的有效的差分逼近的处理方法,得到了该方法的三。一模误差估计．     

纵向资料回归系数的新估计

本文研究纵向资料中回归系数的估计,借助于相依回归模型中的信息迭加法,给出了一个新估计,它改进了最小二乘估计．     

Unbiasedness of the Estimator of the Function of Expected Value in the Mixed Linear Model

The traditional method for estimating the linear function of fixed parameters in mixed linear model is a two-stage procedure. In the first stage of this procedure the variance components estimators are calculated and next in the second stage these estimators are taken as true values of variance components to estimating the linear function of fixed parameters according to generalized least squares method. In this paper the general mixed linear model is considered in which a matrix related to fixed parameters and or/a dispersion matrix of observation vector may be deficient in rank. It is shown that the estimators of a set of functions of fixed parameters obtained in second stage are unbiased if only the observation vector is symmetrically distributed about its expected value and the estimators of variance components from first stage are translation-invariant and are even functions of the observation vector.     

A general approach for two-stage analysis of multilevel clustered non-Gaussian data

In this article, we propose a two-stage approach to modeling multilevel clustered non-Gaussian data with sufficiently large numbers of continuous measures per cluster. Such data are common in biological and medical studies utilizing monitoring or image-processing equipment. We consider a general class of hierarchical models that generalizes the model in the global two-stage (GTS) method for nonlinear mixed effects models by using any square-root-n-consistent and asymptotically normal estimators from stage 1 as pseudodata in the stage 2 model, and by extending the stage 2 model to accommodate random effects from multiple levels of clustering. The second-stage model is a standard linear mixed effects model with normal random effects, but the cluster-specific distributions, conditional on random effects, can be non-Gaussian. This methodology provides a flexible framework for modeling not only a location parameter but also other characteristics of conditional distributions that may be of specific interest. For estimation of the population parameters, we propose a conditional restricted maximum likelihood (CREML) approach and establish the asymptotic properties of the CREML estimators. The proposed general approach is illustrated using quartiles as cluster-specific parameters estimated in the first stage, and applied to the data example from a collagen fibril development study. We demonstrate using simulations that in samples with small numbers of independent clusters, the CREML estimators may perform better than conditional maximum likelihood estimators, which are a direct extension of the estimators from the GTS method.     

Estimating Population Size and Density Using Line Transect Sampling

ANDERSON and POSPAHALA (1970) investigated the estimation of wildlife population size using the belt or line transect sampling method and devised a correction for bias, thus leading to a class of estimators with desirable characteristics. This work was given a basic and rigorous mathematica framework by BURNHAM and ANDERSON (1976). In the present article we use this mathematical framework to develop an estimator of population size and density using weighted least squares. The approach is a two-stage Method.     

Semiparametric additive time-varying coefficients model for longitudinal data with censored time origin

Statistical analysis of longitudinal data often involves modeling treatment effects on clinically relevant longitudinal biomarkers since an initial event (the time origin). In some studies including preventive HIV vaccine efficacy trials, some participants have biomarkers measured starting at the time origin, whereas others have biomarkers measured starting later with the time origin unknown. The semiparametric additive time-varying coefficient model is investigated where the effects of some covariates vary nonparametrically with time while the effects of others remain constant. Weighted profile least squares estimators coupled with kernel smoothing are developed. The method uses the expectation maximization approach to deal with the censored time origin. The Kaplan–Meier estimator and other failure time regression models such as the Cox model can be utilized to estimate the distribution and the conditional distribution of left censored event time related to the censored time origin. Asymptotic properties of the parametric and nonparametric estimators and consistent asymptotic variance estimators are derived. A two-stage estimation procedure for choosing weight is proposed to improve estimation efficiency. Numerical simulations are conducted to examine finite sample properties of the proposed estimators. The simulation results show that the theory and methods work well. The efficiency gain of the two-stage estimation procedure depends on the distribution of the longitudinal error processes. The method is applied to analyze data from the Merck 023/HVTN 502 Step HIV vaccine study.     

Closed Form Estimators of Latent Cell Frequencies

One-stage and two-stage closed form estimators of latent cell frequencies in multidimensional contingency tables are derived from the weighted least squares criterion. The first stage estimator is asymptotically equivalent to the conditional maximum likelihood estimator and does not necessarily have minimum asymptotic variance. The second stage estimator does have minimum asymptotic variance relative to any other existing estimator. The closed form estimators are defined for any number of latent cells in contingency tables of any order under exact general linear constraints on the logarithms of the nonlatent and latent cell frequencies.     

Structural Nested Mean Models for Assessing Time-Varying Effect Moderation

Summary .  This article considers the problem of assessing causal effect moderation in longitudinal settings in which treatment (or exposure) is time varying and so are the covariates said to moderate its effect.  Intermediate causal effects  that describe time-varying causal effects of treatment conditional on past covariate history are introduced and considered as part of Robins' structural nested mean model. Two estimators of the intermediate causal effects, and their standard errors, are presented and discussed: The first is a proposed two-stage regression estimator. The second is Robins' G-estimator. The results of a small simulation study that begins to shed light on the small versus large sample performance of the estimators, and on the bias–variance trade-off between the two estimators are presented. The methodology is illustrated using longitudinal data from a depression study.     

Joint analysis of longitudinal data with informative right censoring

Longitudinal data arise when subjects are followed over a period of time. A commonly encountered complication in the analysis of such data is the variable length of follow-up due to right censorship. This can be further exacerbated by the possible dependency between the censoring time and the longitudinal measurements. This article proposes a combination of a semiparametric transformation model for the censoring time and a linear mixed effects model for the longitudinal measurements. The dependency is handled via latent variables which are naturally incorporated. We show that the likelihood function has an explicit form and develops a two-stage estimation procedure to avoid direct maximization over a high-dimensional parameter space. The resulting estimators are shown to be consistent and asymptotically normal, with a closed form for the variance-covariance matrix that can be used to obtain a plug-in estimator. Finite sample performance of the proposed approach is assessed through extensive simulations. The method is applied to renal disease data.     

Estimating selective advantage of two alleles in discrete time

A class of estimators for the selective advantage, s, in a Wright-Fisher model with two alleles, variable population size, and genic selection is derived via martingale theory. Explicit expressions are given for these estimators which only involve simple computation. The optimal estimate among this class of estimators is obtained. Asymptotic results are readily established by an application of a martingale central limit theorem. The performance of this optimal estimator is compared to known estimators by means of a simulation study.     

Fast approximations of pseudo-observations in the context of right censoring and interval censoring

In the context of right-censored and interval-censored data, we develop asymptotic formulas to compute pseudo-observations for the survival function and the restricted mean survival time (RMST). These formulas are based on the original estimators and do not involve computation of the jackknife estimators. For right-censored data, Von Mises expansions of the Kaplan–Meier estimator are used to derive the pseudo-observations. For interval-censored data, a general class of parametric models for the survival function is studied. An asymptotic representation of the pseudo-observations is derived involving the Hessian matrix and the score vector. Theoretical results that justify the use of pseudo-observations in regression are also derived. The formula is illustrated on the piecewise-constant-hazard model for the RMST. The proposed approximations are extremely accurate, even for small sample sizes, as illustrated by Monte Carlo simulations and real data. We also study the gain in terms of computation time, as compared to the original jackknife method, which can be substantial for a large dataset.     

Statistical inference and power analysis for direct and spillover effects in two-stage randomized experiments

Two-stage randomized experiments become an increasingly popular experimental design for causal inference when the outcome of one unit may be affected by the treatment assignments of other units in the same cluster. In this paper, we provide a methodological framework for general tools of statistical inference and power analysis for two-stage randomized experiments. Under the randomization-based framework, we consider the estimation of a new direct effect of interest as well as the average direct and spillover effects studied in the literature. We provide unbiased estimators of these causal quantities and their conservative variance estimators in a general setting. Using these results, we then develop hypothesis testing procedures and derive sample size formulas. We theoretically compare the two-stage randomized design with the completely randomized and cluster randomized designs, which represent two limiting designs. Finally, we conduct simulation studies to evaluate the empirical performance of our sample size formulas. For empirical illustration, the proposed methodology is applied to the randomized evaluation of the Indian National Health Insurance Program. An open-source software package is available for implementing the proposed methodology.