Marginally oversaturated designs

A special class of supersaturated design, called marginally over saturated design (MOSD), in which the number of variables under investigation (k) is only slightly larger than the number of experimental runs (n), is presented. Several optimality criteria for supersaturated designs are discussed. It is shown that the resolution rank criterion is most appropriate for screening situations. The construction method builds on two major theorems which provide an efficient way to evaluate resolution rank. Examples are given for the cases n=8, 12, 16, and 20. Potential extensions for future work are discussed.     

Influence of human immunodeficiency virus infection on neurological impairment: an analysis of longitudinal binary data with informative drop-out

A study to investigate the human immunodeficiency virus (HIV) status on the course of neurological impairment, conducted by the HIV Center at Columbia University, followed a cohort of HIV positive and negative gay men for 5 years and assessed the presence or absence of neurological impairment every 6 months. Almost half of the subjects dropped out before the end of the study for reasons that might have been related to the missing neurological data. We propose likelihood-based methods for analysing such binary longitudinal data under informative and non-informative drop-out. A transition model is assumed for the binary response, and several models for the drop-out processes are considered which are functions of the response variable (neurological impairment). The likelihood ratio test is used to compare models with informative and non-informative drop-out mechanisms. Using simulations, we investigate the percentage bias and mean-squared error (MSE) of the parameter estimates in the transition model under various assumptions for the drop-out. We find evidence for informative drop-out in the study, and we illustrate that the bias and MSE for the parameters of the transition model are not directly related to the observed drop-out or missing data rates. The effect of HIV status on the neurological impairment is found to be statistically significant under each of the models considered for the drop-out, although the regression coefficient may be biased in certain cases. The presence and relative magnitude of the bias depend on factors such as the probability of drop-out conditional on the presence of neurological impairment and the prevalence of neurological impairment in the population under study.     

基于线性矩阵不等式的牵制控制策略分析

研究复杂动力学网络的同步牵制控制问题.通过引入牵制控制,对具有大规模节点网络中的小部分节点实施线性反馈控制策略,使得整个网络的所有节点均达到同步.针对一类一致连结的耦合复杂网络模型,采用线性矩阵不等式的方法给出了实施牵制控制达到同步的充分条件.该条件表明,选择牵制异质网络中度较大的节点可能更容易使得网络实现同步.特别地,对于特定的内部耦合矩阵和足够大的耦合强度,只要对网络中任一节点实施牵制控制就可以保证整个网络同步的实现.数值仿真验证了上述结论的有效性.     

A study of BIB designs through support matrices

This paper deals with the existence and nonexistence of BIB designs with repeated blocks. The approach is an algebraic one. The concept of a support matrix is introduced and some of its basic properties are noted. Some basic examples of support matrices are given when the block size is 3. The connection between full column rank proper support matrices and irreducible designs is explored and some examples of such matrices are given.     

Estimation of parameters in incomplete data models defined by dynamical systems

Parametric incomplete data models defined by ordinary differential equations (ODEs) are widely used in biostatistics to describe biological processes accurately. Their parameters are estimated on approximate models, whose regression functions are evaluated by a numerical integration method. Accurate and efficient estimations of these parameters are critical issues. This paper proposes parameter estimation methods involving either a stochastic approximation EM algorithm (SAEM) in the maximum likelihood estimation, or a Gibbs sampler in the Bayesian approach. Both algorithms involve the simulation of non-observed data with conditional distributions using Hastings–Metropolis (H–M) algorithms. A modified H–M algorithm, including an original local linearization scheme to solve the ODEs, is proposed to reduce the computational time significantly. The convergence on the approximate model of all these algorithms is proved. The errors induced by the numerical solving method on the conditional distribution, the likelihood and the posterior distribution are bounded. The Bayesian and maximum likelihood estimation methods are illustrated on a simulated pharmacokinetic nonlinear mixed-effects model defined by an ODE. Simulation results illustrate the ability of these algorithms to provide accurate estimates.     

AN ALGORITHMIC APPROACH TO BAYESIAN LINEAR MODEL CALCULATIONS

A framework is described for organizing and understanding the computations necessary to obtain the posterior mean of a vector of linear effects in a normal linear model, conditional on the parameters that determine covariance structure. The approach has two major uses; firstly, as a pedagogical tool in the derivation of formulae, and secondly, as a practical tool for developing computational strategies without needing complicated matrix formulae that are often unwieldy in complex hierarchical models. The proposed technique is based upon symbolic application of the sweep operator SWP to an appropriate tableau of means and covariances. The method is illustrated with standard linear model specifications, including the so-called mixed model, with both fixed and random effects.     

A Continuous Representation of the Family of Stable Law Distributions

Conventional parametric representations of stable law distributions do not allow all members of the family to be obtained as continuous limits of the parameters. Model building (or simulation) using such representations will be numerically unstable near such limits in consequence. Existing tables are not satisfactory near such limits as interpolation cannot be carried out. We show that these difficulties are overcome by using a new shifted Cartesian representation which characterizes the entire stable law family in a completely continuous way. Standardization is still possible with this representation so that tabulation, using just two bounded parameters, can be carried out. Its use is illustrated in a non-regular threshold estimation problem involving stable distributions which are discontinuous limits in conventional representations.     

一种基于冗余小波变换的DT网格运动估计和运动补偿方法

网格划分和网格顶点的运动估计是基于不规则网格的视频压缩技术的关键。为了进一步提高网格运动估计和运动补偿的效果,在综合比较现有冗余小波变换域运动估计方法和适用规则网格的EMRMC算法的基础上,提出了一种新的基于不规则网格的运动估计和运动补偿算法,即在冗余小波变换域提取特征点和运动潜在区,网格顶点的运动估计采用结合运动潜在区的在时域进行块匹配的运动估计和运动补偿方法,而运动补偿则通过三角形仿射变换完成。同时还给出了冗余小波变换域提取运动潜在区的计算模板。理论分析和实验结果表明,该算法在补偿效果方面较前两种方法得到了改进。     

二次型与(半)稳定矩阵的判定

对二次型与（半）稳定矩阵之间的相互关系进行了讨论，得到了二者之间的相互判定方法，且利用二次型的不同配方法得到判定（半）稳定矩阵的一些条件     

开放微带不连续电路的三端口网络描述

给出了开放微带不连续电路的三端口网络描述。将微带电路向空间的辐射作为一个端口来处理，成功地消除了由于忽略电路的空间辐射效应引入的电路数字模拟误差。利用时域有限差分法（ＦＤＴＤ）一次运算就可得出三端口网络的Ｓ参数矩阵的各元素，而无需任何预处理，大大地节省了ＣＰＵ时间。实际计算与已发表的数据吻合得很好，表明文中给出的三端口定义是合理有效的。此三端口微带电路的定义可容易地推广到其他非屏蔽传输线电路的分析，可广泛用于电磁兼容分析、微带天线及其他开放传输线电路的ＣＡＤ技术中。