Weighting variables in Kohonen competitive learning algorithms

This paper presents a new variable weight method, called the singular value decomposition (SVD) approach, for Kohonen competitive learning (KCL) algorithms based on the concept of Varshavsky et al. [18 R. Varshavsky, A. Gottlieb, M. Linial, and D. Horn, Novel unsupervised feature filtering of bilogical data, Bioinformatics 22 (2006), pp. 507–513.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]]. Integrating the weighted fuzzy c-means (FCM) algorithm with KCL, in this paper, we propose a weighted fuzzy KCL (WFKCL) algorithm. The goal of the proposed WFKCL algorithm is to reduce the clustering error rate when data contain some noise variables. Compared with the k-means, FCM and KCL with existing variable-weight methods, the proposed WFKCL algorithm with the proposed SVD's weight method provides a better clustering performance based on the error rate criterion. Furthermore, the complexity of the proposed SVD's approach is less than Pal et al. [17 S.K. Pal, R.K. De, and J. Basak, Unsupervised feature evaluation: a neuro-fuzzy approach, IEEE. Trans. Neural Netw. 11 (2000), pp. 366–376.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]], Wang et al. [19 X.Z. Wang, Y.D. Wang, and L.J. Wang, Improving fuzzy c-means clustering based on feature-weight learning, Pattern Recognit. Lett. 25 (2004), pp. 1123–1132.[Crossref], [Web of Science ®] , [Google Scholar]] and Hung et al. [9 W. -L. Hung, M. -S. Yang, and D. -H. Chen, Bootstrapping approach to feature-weight selection in fuzzy c-means algorithms with an application in color image segmentation, Pattern Recognit. Lett. 29 (2008), pp. 1317–1325.[Crossref], [Web of Science ®] , [Google Scholar]].     

Is Each NPMLE of a Continuous Bivariate Distribution Function with Singly Right-Censored Data Really Inconsistent?

We consider non-parametric estimation of a continuous cdf of a random vector (X ₁, X ₂). With bivariate RC data, it is stated in van der Laan (1996 Van der Laan , M. J. ( 1996 ) Efficient estimation in the bivariate censoring model and repairing NPMLE . Ann. Statist. 24 : 596 – 627 .[Crossref], [Web of Science ®] , [Google Scholar], p. 598¹⁰, Ann. Statist.), Quale et al. (2006 Quale , C. M. , van der Laan , M. J. , Robins , J. R. ( 2006 ). Locally efficient estimation with bivariate right-censored data . JASA. 101 : 1076 – 1084 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar], JASA) etc. that “it is well known that the NPMLE for continuous data is inconsistent (Tsai et al. (1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar])).” The claim is based on a result in Tsai et al. (1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar], p.1352, Ann. Statist.) that if X ₁ is right censored but not X ₂, then common ways for defining one NPMLE lead to inconsistency. If X ₁ is right censored and X ₂ is type I right-censored (which includes the case in Tsai et al.), we present a consistent NPMLE. The result corrects a common misinterpretation of Tsai's example (Tsai et al., 1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar], Ann. Statist.).     

Additive Mixed Effect Model for Clustered Doubly Censored Data

Cai and Zeng (2011 Cai, J. and Zeng, D. 2011. Additive mixed effect model for clustered failure time data. Biometrics, 67(4): 1340–1351. [Crossref], [PubMed] , [Google Scholar]) proposed an additive mixed effect model to analyze clustered right-censored data. In this article, we demonstrate that the approach of Cai and Zeng (2011 Cai, J. and Zeng, D. 2011. Additive mixed effect model for clustered failure time data. Biometrics, 67(4): 1340–1351. [Crossref], [PubMed] , [Google Scholar]) can be extended to clustered doubly censored data. Furthermore, when both left- and right-censoring variables are always observed, we propose alternative estimators using the approach of Cai and Cheng (2004 Cai, T. and Cheng, S. C. 2004. Semiparametric regression analysis for doubly censored data. Biometrika, 91: 277–290. [Crossref], [Web of Science ®] , [Google Scholar]). A simulation study is conducted to investigate the performance of the proposed estimators.     

A Generalized Time-Dependent Conditional Linear Model with Left-Truncated and Right-Censored Data

Consider the model φ(S(y | X)) = β(y) ^T X, where φ is a known link function, S(· | X) is the survival function of a response Y given a covariate X = (1, X, X ²,…, X ^p ), and β(y) is an unknown vector of time-dependent regression coefficients. The response Y is subject to left truncation and right censoring. We assume that given X, Y is independent of (C, T) where C and T are censoring and truncation variables with P(C ≥ T) = 1. In this article, with some modification of the assumptions in Lemmas 5 and 6 of Iglesias-Pérez and González-Manteiga (1999 Iglesias-Pérez , C. J. , González-Manteiga , W. G. ( 1999 ). Strong representation of a generalized product-limit estimator for truncated and censored data with some application . J. Nonparametric Statist. 10 : 213 – 244 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), we present an almost sure representation for the generalized product-limit estimator (GPL) of S(y | X). Based on the GPL and the approach of Teodorescu et al. (2010 Teodorescu , B. , Keilegom , I. V. , Cao , R. ( 2010 ). Generalized time-dependent conditional linear models under left truncation and right censoring . Ann. Instit. Statist. Math. 62 : 465 – 485 .[Crossref], [Web of Science ®] , [Google Scholar]), a least squares estimator of β(y) is obtained and a bootstrap procedure is proposed to choose the optimum bandwidth.     

Computing the Exact Fisher Information Matrix of Periodic State-Space Models

The purpose of this article is to develop algorithms for computing the exact Fisher information matrix of periodic time-varying state-space models. We first present a relatively simple recursive algorithm which computes the elements of the exact information matrix without involving numerical differentiation, since all required derivatives are analytically evaluated. The proposed algorithm extends the procedure due to Cavanaugh and Shumway (1996 Cavanaugh , J. E. , Shumway , R. H. ( 1996 ). On computing the expected Fisher information matrix for state-space model parameters . Statist. Probab. Lett. 26 : 347 – 355 .[Crossref], [Web of Science ®] , [Google Scholar]) to the periodic state-space framework. Exploiting the approach used in Klein et al. (2000 Klein , A. , Mélard , G. , Zahaf , T. ( 2000 ). Construction of the exact Fisher information matrix of Gaussian time series models by means of matrix differential rules . Linear Alg. Applic. 321 : 209 – 232 .[Crossref], [Web of Science ®] , [Google Scholar]), a second algorithm is proposed in order to obtain the exact information matrix as a whole instead of element by element. The algorithms are first developed in a general framework and then specialized to the case of a periodic Gaussian vector autoregressive moving-average (PVARMA) model.     

Semiparametric analysis of transformation models with doubly censored data

Double censoring arises when T represents an outcome variable that can only be accurately measured within a certain range, [L, U], where L and U are the left- and right-censoring variables, respectively. In this note, using Martingale arguments of Chen et al. [3 Chen, K., Jin, Z. and Ying, Z. 2002. Semiparametric analysis of transformation models with censored data. Biometrika, 89: 659–668. [Crossref], [Web of Science ®] [Google Scholar]], we propose an estimator (denoted by ?β) for estimating regression coefficients of transformation model when L is always observed. Under Cox proportional hazards model, the proposed estimator is equivalent to the partial likelihood estimator for left-truncated and right-censored data if the left-censoring variables L were regarded as left-truncated variables. In this case, the estimator ?β can be obtained by the standard software. A simulation study is conducted to investigate the performance of ?β. For the purpose of comparison, the simulation study also includes the estimator proposed by Cai and Cheng [2 Cai, T. and Cheng, S. 2004. Semiparametric regression analysis for doubly censored data. Biometrika, 91: 277–290. [Crossref], [Web of Science ®] [Google Scholar]] for the case when L and U are always observed.     

Aspects of Dependence in Generalized Farlie-Gumbel-Morgenstern Distributions

In this article, we obtained a dependence measure for generalized Farlie-Gumbel-Morgenstern (FGM) family in view of Kochar and Gupta (1987 Kochar , S. G. , Gupta , R. P. ( 1987 ). Competitors of Kendall-tau test for testing independence against PQD . Biometrika 74 ( 3 ): 664 – 669 .[Crossref], [Web of Science ®] , [Google Scholar]) and then compared this measure with Spearman's rho and Kendall's tau in FGM family. Moreover, we evaluated the empirical power of the class of distribution-free tests proposed by Kochar and Gupta (1987 Kochar , S. G. , Gupta , R. P. ( 1987 ). Competitors of Kendall-tau test for testing independence against PQD . Biometrika 74 ( 3 ): 664 – 669 .[Crossref], [Web of Science ®] , [Google Scholar], 1990 Kochar , S. G. , Gupta , R. P. ( 1990 ). Distribution-free tests based on sub-sample extrema for testing against positive dependence . Australian Journal of Statistics 32 : 45 – 51 .[Crossref] , [Google Scholar]) based on exact distribution of a U-statistics. This is derived via a simulation study for sample of sizes n = 6, 8, 10, 12, 16, and 20. Also, we compared our simulation results with those achieved by Amini et al. (2010 Amini , M. , Jabbari , H. , Mohtashami Borzadaran , G. R. , Azadbakhsh , M. ( 2010 ). Power comparison of independence test for the Farlie-Gumbel-Moregenstern family . Communications of the Korean Statistical Society 17 ( 4 ): 493 – 505 .[Crossref] , [Google Scholar]) and Güven and Kotz (2008 Güven , B. , Kotz , S. ( 2008 ). Test of independence for generalized Farlie-Gumbel-Morgenstern distributions . Journal of Computational and Applied Mathematics 212 : 102 – 111 .[Crossref], [Web of Science ®] , [Google Scholar]).     

Procedures for Selecting Good Exponential Populations

Suppose exponential populations π _i with parameters (μ_i , σ_i ) (i = 1,2,…,k) are given. This article discusses how to select “good” populations in the sense of [Lam (1986 Lam, K. 1986. A new procedure for selecting good populations. Biometrika, 73(1): 201–206. [Crossref], [Web of Science ®] , [Google Scholar]). A new procedure for selecting good populations. Biometrika 73(1):201–206]. Depending on whether the σ _i 's are known or unknown, several one-stage and a two-stage procedure of selection are proposed. The two-stage procedure can be replaced by a one-stage procedure if the second-stage sample is proved intangible. An attracting feature of these procedures is that they need no new statistical tables to implement.     

D-Optimal Chemical Balance Weighing Designs with N ≡ 0 (mod 4) and 3 Objects

In this article, the problem of estimation of the individual weights of three objects using a chemical balance weighing design is considered. We use the criterion of D-optimality. We assume that the covariance matrix of errors is the matrix of first-order autoregressive process. Such problems were discussed in Li and Yang (2005 Li , C. H. , Yang , S. Y. ( 2005 ). On a conjecture in D-optimal designs with n ≡ 0 (mod 4) . Lin. Alg. Applic. 400 : 279 – 290 .[Crossref], [Web of Science ®] , [Google Scholar]) and also in Yeh and Lo Huang (2005 Yeh , H. G. , Lo Huang , M. N. ( 2005 ). On exact D-optimal designs with 2 two-level factors and n autocorrelated observations . Metrika 61 : 261 – 275 .[Crossref], [Web of Science ®] , [Google Scholar]). We present some results of D-optimal designs in certain class of designs with the design matrix X  ∈ M _n×3(±1) such that each column of matrix X has at least one 1 and one ?1.     

More Efficient Tests Robust to Heteroskedasticity of Unknown Form

ABSTRACT

In the presence of heteroskedasticity of unknown form, the Ordinary Least Squares parameter estimator becomes inefficient, and its covariance matrix estimator inconsistent. Eicker (1963 Eicker , B. ( 1963 ). Limit theorems for regression with unequal and dependant errors . Ann. Math. Statist. 34 : 447 – 456 .[Crossref] , [Google Scholar]) and White (1980 White , H. ( 1980 ). A heteroskedasticity-consistent covariance matrix estimator and a direct test for heteroskedasticity . Econometrica 48 : 817 – 838 .[Crossref], [Web of Science ®] , [Google Scholar]) were the first to propose a robust consistent covariance matrix estimator, that permits asymptotically correct inference. This estimator is widely used in practice. Cragg (1983 Cragg , J. G. ( 1983 ). More efficient estimation in the presence of heteroskedasticity of unknown form . Econometrica 51 : 751 – 63 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed a more efficient estimator, but concluded that tests basd on it are unreliable. Thus, this last estimator has not been used in practice. This article is concerned with finite sample properties of tests robust to heteroskedasticity of unknown form. Our results suggest that reliable and more efficient tests can be obtained with the Cragg estimators in small samples.     

Bayesian Inference in Dynamic Disequilibrium Models: An Application to the Polish Credit Market

We propose a Bayesian approach for inference in a dynamic disequilibrium model. To circumvent the difficulties raised by the Maddala and Nelson (1974 Maddala , G. , Nelson , F. ( 1974 ). Maximum likelihood methods for models of markets in disequilibrium . Econometrica 42 ( 6 ): 1013 – 1030 .[Crossref], [Web of Science ®] , [Google Scholar]) specification in the dynamic case, we analyze a dynamic extended version of the disequilibrium model of Ginsburgh et al. (1980 Ginsburgh , V. , Tishler , A. , Zang , I. ( 1980 ). Alternative estimation methods for two regime models . European Economic Review 13 : 207 – 228 .[Crossref], [Web of Science ®] , [Google Scholar]). We develop a Gibbs sampler based on the simulation of the missing observations. The feasibility of the approach is illustrated by an empirical analysis of the Polish credit market, for which we conduct a specification search using the posterior deviance criterion of Spiegelhalter et al. (2002 Spiegelhalter , D. J. , Best , N. G. , Carlin , B. P. , Van Der Linde , A. ( 2002 ). Bayesian measures of complexity and fit . Journal of the Royal Statistical Society, Ser. B 64 : 583 – 639 .[Crossref] , [Google Scholar]).     

A Note on Importance Resampling for Multi-Dimensional Statistics

Johns (1988 Johns , M. V. (1988). Importance sampling for bootstrap confidence intervals. Journal of the American Statistical Association 83:709–714.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Davison (1988 Davison , A. C. ( 1988 ). Discussion of paper by D. V. Hinkley . Journal of the Royal Statistical Society Series B 50 : 356 – 357 . [Google Scholar]), and Do and Hall (1991 Do , K. A. , Hall , P. ( 1991 ). On importance sampling for the bootstrap . Biometrika 78 : 161 – 167 .[Crossref], [Web of Science ®] , [Google Scholar]) used importance sampling for calculating bootstrap distributions of one-dimensional statistics. Realizing that their methods can not be extended easily to multi-dimensional statistics, Fuh and Hu (2004 Fuh , C. D. , Hu , I. ( 2004 ). Efficient importance sampling for events of moderate deviations with applications . Biometrika 91 : 471 – 490 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed an exponential tilting formula for statistics of multi-dimension, which is optimal in the sense that the asymptotic variance is minimized for estimating tail probabilities of asymptotically normal statistics. For one-dimensional statistics, Hu and Su (2008 Hu , J. , Su , Z. ( 2008 ). Adaptive resampling algorithms for estimating bootstrap distributions . Journal of Statistical Planning and Inference 138 ( 6 ): 1763 – 1777 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed a multi-step variance minimization approach that can be viewed as a generalization of the two-step variance minimization approach proposed by Do and Hall (1991 Do , K. A. , Hall , P. ( 1991 ). On importance sampling for the bootstrap . Biometrika 78 : 161 – 167 .[Crossref], [Web of Science ®] , [Google Scholar]). In this article, we generalize the approach of Hu and Su (2008 Hu , J. , Su , Z. ( 2008 ). Adaptive resampling algorithms for estimating bootstrap distributions . Journal of Statistical Planning and Inference 138 ( 6 ): 1763 – 1777 .[Crossref], [Web of Science ®] , [Google Scholar]) to multi-dimensional statistics, which applies to general statistics and does not resort to asymptotics. Empirical results on a real survival data set show that the proposed algorithm provides significant computational efficiency gains.     

On a Generalization of Bivariate Cauchy Distribution

This paper addresses a generalization of the bivariate Cauchy distribution discussed by Fang et al. (1990 Fang , K. T. , Kotz , S. , Ng , K. W. ( 1990 ). Symmetric Multivariate and Related Distributions . London : Chapman and Hall .[Crossref] , [Google Scholar]), derived from a trivariate normal distribution with a general correlation matrix. We obtain explicit expressions for the joint distribution function and joint density function, and show that they reduce in a special case to the corresponding expressions of Fang et al. (1990 Fang , K. T. , Kotz , S. , Ng , K. W. ( 1990 ). Symmetric Multivariate and Related Distributions . London : Chapman and Hall .[Crossref] , [Google Scholar]). Finally, we show that this generalized distribution is useful in determining the orthant probability of a bivariate skew-normal distribution of Azzalini and Dalla Valle (1996 Azzalini , A. , Dalla Valle , A. ( 1996 ). The multivariate skew-normal distribution . Biometrika 83 : 715 – 726 .[Crossref], [Web of Science ®] , [Google Scholar]).     

Change Point Detection in a General Class of Distributions

ABSTRACT

We consider the change point problem in a general class of distributions, and derive a test statistic T _n which reduces to the statistic obtained by Kander and Zacks (1966 Kander , Z. , Zacks , S. ( 1966 ). Test procedure for possible changes in parameter of statistical distributions occurring at unknown time points. Ann. Math. Statist. 37 : 1196 – 1210 . [CSA] [Crossref] , [Google Scholar]) for the exponential family. Properties of the test, including its asymptotic distribution, are discussed.     

Improved Bar-Lev,Bobovitch, and Boukai Randomized Response Models

In this article, two new improved randomized response models have been proposed. The proposed models are found to be more efficient than the recent randomized response model studied by Bar-Lev et al. (2004 Bar-Lev , S. K. , Bobovitch , E. , Boukai , B. ( 2004 ). A note on randomized response models for quantitative data . Metrika 60 : 255 – 260 .[Crossref], [Web of Science ®] , [Google Scholar]). The relative efficiency of the proposed models has been studied with respect to the Bar-Lev et al. (2004 Bar-Lev , S. K. , Bobovitch , E. , Boukai , B. ( 2004 ). A note on randomized response models for quantitative data . Metrika 60 : 255 – 260 .[Crossref], [Web of Science ®] , [Google Scholar]) model under different situations.     

Upper Bounds for the Ruin Probabilities of the Entrance-Based Risk Model

Li et al. (2005 Li , Z. , Zhu , J. , Chen , F. ( 2005 ). Study of a risk model based on the entrance process . Statist. Probab. Lett. 72 ( 1 ): 1 – 10 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed a risk model based on the entrance process and studied the asymptotic behavior of the surplus when time goes to infinity. This article considers the ruin problem in that model. Some simple characteristics (stochastic intensity, compensator, mean process, etc.) of the risk process and other related processes are also considered. Under small claim condition, exponential upper bounds for the ruin probability are obtained.     

Optimum Mixture Design Using Deficiency Criterion

In a mixture experiment the measured response is assumed to depend only on the relative proportion of ingredients or components present in the mixture. Scheffe (1958 Scheffe , H. ( 1958 ). Experiments with mixtures . J. R. Statist. Soc. B 20 : 344 – 360 . [Google Scholar], 1963 Scheffe , H. ( 1963 ). Simplex – centroid design for experiments with mixtures . J. R. Statist. Soc. B 25 : 235 – 263 . [Google Scholar]) first systematically considered this problem and introduced different models and designs suitable in such situations. Optimum designs for the estimation of parameters of different mixture models are available in the literature. The problem of estimating the optimum proportion of mixture components is of great practical importance. Pal and Mandal (2006 Pal , M. , Mandal , N. K. ( 2006 ). Optimum designs for optimum mixtures . Statist. Probab. Lett. 76 ( 13 ): 1369 – 1379 .[Crossref], [Web of Science ®] , [Google Scholar], 2007 Pal , M. , Mandal , N. K. (2007). Optimum mixture design via equivalence theory. communicated. JCISS 32:107–126. [Google Scholar]) attempted to find a solution to this problem by adopting a pseudo-Bayesian approach and using the trace criterion. Subsequently, Pal and Mandal (2008 Pal , M. , Mandal , N. K. ( 2008 ). Minimax designs for optimum mixtures . To appear in Statist. Probab. Lett.  [Google Scholar]) solved the problem using minimax criterion. In this article, the deficiency criterion due to Chatterjee and Mandal (1981 Chatterjee , S. K. , Mandal , N. K. ( 1981 ). Response surface designs for estimating the optimal point . Calcutta Statist. Assoc. Bull. . 30 : 145 – 169 .[Crossref] , [Google Scholar]) has been used as a measure for comparing the performance of competing designs.     

Inverse Sampling with Unequal Selection Probabilities

In some real situations the population of interest is divided into two groups, of which one contains only a few units. In other cases, the population may be considered as subdivided into two group', for example, if only a few units display a value of the variable of interest which is highly different from zero, while all the other units show a value equal to or near zero. In both cases, inverse sampling is more efficient than classical fixed sample-size designs to obtain the parameter estimators for the whole population as well as for its groups (e.g., Salehi and Seber, 2004 Salehi , M. M. , Seber , G. A. F. ( 2004 ). A general inverse sampling scheme and its application to adaptive cluster sampling . Austral. NZ J. Statist. 46 : 483 – 494 .[Crossref], [Web of Science ®] , [Google Scholar]). In fact, in this design the procedure selection continues until a prefixed number of units with the characteristic of interest is sampled. Since it is not known a priori to which group the population units belong, the sample size is a random variable. Christman and Lan (2001 Christman , M. C. , Lan , F. ( 2001 ). Inverse adaptive cluster sampling . Biometrics 57 : 1096 – 1105 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Salehi and Seber (2001 Salehi , M. M. , Seber , G. A. F. ( 2001 ). A new proof of Murthy's estimator which applies to sequential sampling . Austral. NZ J. Statist. 43 : 281 – 286 . [Google Scholar] 2004 Salehi , M. M. , Seber , G. A. F. ( 2004 ). A general inverse sampling scheme and its application to adaptive cluster sampling . Austral. NZ J. Statist. 46 : 483 – 494 .[Crossref], [Web of Science ®] , [Google Scholar]) considered inverse sampling designs when all the population units have equal selection probabilities. In this article, we consider the general case in which the units may have unequal probabilities of being included in the sample. In fact, in many real situations different units may have different selection probabilities because of some inherent features of the sampling procedure, or in order to obtain better estimates. We derive unbiased estimators of the totals of the two groups, their variance and the corresponding unbiased variance estimators in inverse sampling with replacement. Finally, we derive similar results for more complex designs, where the selection procedure stops before observing the prefixed number of units from the rare group.     

Regression Analysis with Covariates Missing at Random: A Piece-wise Nonparametric Model for Missing Covariates

Statistical analysis for the regression model f _β(y | x, z) with missing values in the covariate vector X requires modeling of the covariate distribution g(x | z). Likelihood methods, including Ibrahim (1990 Ibrahim , J. G. ( 1990 ). Incomplete data in generalized linear models . J. Amer. Statist. Assoc. 85 : 765 – 769 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Chen (2004 Chen , H. Y. (2004). Nonparametric and semiparametric models for missing covariates in parametric regression. J. Amer. Statist. Assoc. 99:1176–1189.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), and Zhao (2005 Zhao , Y. ( 2005 ). Design and Efficient Estimation in Regression Analysis with Missing Data in Two-Phase Studies. Ph.D. thesis , University of Waterloo . [Google Scholar]), need either X or Z to be discrete. This article considers extending the likelihood methods to deal with cases where both X and Z may be continuous. We propose modeling the covariate distribution g(x | z) using a piece-wise nonparametric model, then a maximum likelihood estimate (MLE) of β can be computed following the maximum likelihood estimating procedure of Chen (2004 Chen , H. Y. (2004). Nonparametric and semiparametric models for missing covariates in parametric regression. J. Amer. Statist. Assoc. 99:1176–1189.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) or Zhao (2005 Zhao , Y. ( 2005 ). Design and Efficient Estimation in Regression Analysis with Missing Data in Two-Phase Studies. Ph.D. thesis , University of Waterloo . [Google Scholar]). The resulting estimation method is easy to implement and the asymptotic properties of the MLE follow under certain conditions. Extensive simulation studies for different models indicate that the proposed method is acceptable for practical implementation. A real data example is used to illustrate the method.     

Optimal Designs for Best Linear Unbiased Prediction in Diallel Crosses

ABSTRACT

Recently, Ghosh and Das (2003 Ghosh , H. , Das , A. ( 2003 ). Optimal diallel cross designs for estimation of heritability . J. Statist. Plann. Inference 116 : 185 – 196 . [Google Scholar]) considered the estimation of variance components and the variances of these estimates. While comparing the yielding capacities of the cross (i, j), Kempthorne and Curnow (1961 Kempthorne , O. , Curnow , R. N. ( 1961 ). The partial diallel cross . Biometrics 17 : 229 – 250 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed the estimation of the yielding capacity of any cross based on the least square estimators of the general combining ability effects and/or the mean yield of the cross (i, j). In this article, the problem of predicting the yielding capacity of the cross (i, j) from the sample of inbred lines has been considered. The properties of the best linear unbiased predictor for predicting the unobserved general combining ability effects together with general mean effect has been studied. We characterize A-optimal complete diallel cross designs and some efficient partial diallel cross designs under this setup.