• Communications for Statistical Applications and Methods
• /
• v.24 no.6
• /
• pp.663-671
• /
• 2017

Some genetic association tests include an unidentifiable nuisance parameter under the null hypothesis of no association. When the mode of inheritance (MOI) is not specified in a case-control design, the Cochran-Armitage (CA) trend test contains an unidentifiable nuisance parameter. The transmission disequilibrium test (TDT) in a family-based association study that includes the unaffected also contains an unidentifiable nuisance parameter. The hypothesis tests that include an unidentifiable nuisance parameter are typically performed by taking a supremum of the CA tests or TDT over reasonable values of the parameter. The p-values of the supremum test statistics cannot be obtained by a normal or chi-square distribution. A common method is to use a Davies's upper bound of the p-value instead of an exact asymptotic p-value. In this paper, we provide a unified sine-cosine process expression of the CA trend test that does not specify the MOI and the TDT that includes the unaffected. We also present a closed form expression of the exact asymptotic formulas to calculate the p-values of the supremum tests when the score function can be written as a linear form in an unidentifiable parameter. We illustrate how to use the derived formulas using a pharmacogenetics case-control dataset and an attention deficit hyperactivity disorder family-based example.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.201-201
• /
• 2016

The asymptotic extreme value distributions of maxima are a natural choice when designing against future extreme events like flood peaks or wave heights, given a stationary time series. The generalized extreme value distribution (GEV) is often utilised in this context because it is seen as a convenient single expression for extreme event analysis. However, the GEV has a drawback because the location of the distribution bound relative to the data is a discontinuous function of the GEV shape parameter. That is, for annual maxima approximated by the Gumbel distribution, the data is also consistent with a GEV distribution with an upper bound (no lower bound) or a GEV distribution with a lower bound (no upper bound). A more consistent single extreme value expression for design purposes is proposed as the Weibull distribution of smallest extremes, as applied to transformed annual maxima. The Weibull distribution limit holds here for sufficiently large sample sizes, irrespective of the extreme value domain of attraction applicable to the untransformed maxima. The Gumbel, Type 2, and Type 3 extreme value distributions thus become redundant, together with the GEV, because in reality there is only a single asymptotic extreme value distribution required for design purposes - the Weibull distribution of minima as applied to transformed maxima. An illustrative synthetic example is given showing transformed maxima from the normal distribution approaching the Weibull limit much faster than the untransformed sample maxima approach the normal distribution Gumbel limit. Some New Zealand examples are given with the Weibull distribution being applied to reciprocal transformations of annual flood maxima, where the untransformed maxima follow apparently different extreme value distributions.

• Journal of Korea Water Resources Association
• /
• v.42 no.3
• /
• pp.227-234
• /
• 2009

In this study, cross correlations among sample data at each site are calculated to obtain the asymptotic cross correlations among probability weighted moments at each site using Monte Carlo simulation. As a result, the relations between the asymptotic cross correlations among probability weighted moments and the inter-site dependence among sample data at each site are nearly a linear relation with slope 1. The smaller ratio of concurrent data size to entire sample size is, the weaker the relationship grows. Simple power function which the correction term in power function accounts for the differences of the sample size between two sites was fitted to each case to estimate the parameter. It is noted that this result can be used in the various researches which include the estimation of the variance of quantile considering cross correlations.

• Journal of the Korean Data and Information Science Society
• /
• v.16 no.3
• /
• pp.591-601
• /
• 2005

In this paper we consider detection of a discontinuity point in the variance function. When the mean function is discontinuous at a point, the variance function is usually discontinuous at the point. In this case, we had better estimate the location of the discontinuity point with the mean function rather than the variance function. On the other hand, the variance function only has a discontinuity point. The target function in order to estimate the location can be used the second moment function since the variance function and the second moment function have the same location and jump size of the discontinuity point. We propose a nonparametric detection method of the discontinuity point with the second moment function. We give the asymptotic results of these estimators. Computer simulation demonstrates the improved performance of the method over the existing ones.

• Journal of the Korean Data and Information Science Society
• /
• v.4
• /
• pp.1-11
• /
• 1993

The Cox's regression model is frequently used for covariate effects in survival data analysis, But, much of the statistical work has focused on asymptotic behavior so the small sample evaluation has been neglected. In this paper, we compare the small or moderate sample performances of the survival function estimators for the Cox's regression model using bootstrap method. The smoothed PL type estimator and the Link estimator are slightly better than corresponding the PL type estimator and the Nelson type estimator in the sense of the achieved error rates.

• Journal of the Korean Statistical Society
• /
• v.28 no.4
• /
• pp.469-478
• /
• 1999

We propose a Kolmogorov-Smirnov-type test for independence of paired failure times in the presence of independent censoring times. This independent censoring mechanism is often assumed in case-control studies. To do this end, we first introduce a process defined as the difference between the bivariate survival function estimator proposed by Wang and Wells (1997) and the product of the product-limit estimators (Kaplan and Meier (1958)) for the marginal survival functions. Then, we derive its asymptotic properties under the null hypothesis of independence. Finally, we assess the performance of the proposed test by simulations, and illustrate the proposed methodology with a dataset for remission times of 21 pairs of leukemia patients taken from Oakes(1982).

• Communications for Statistical Applications and Methods
• /
• v.22 no.3
• /
• pp.233-239
• /
• 2015

This study compares two generalized Lindley distributions and assesses consistency between theoretical and analytical results. Data (complete and censored) assumed to follow the Lindley distribution are generated and analyzed using two generalized Lindley distributions, and maximum likelihood estimates of parameters from the generalized distributions are obtained. Size and power of tests of hypotheses on the parameters are assessed drawing on asymptotic properties of the maximum likelihood estimates. Results suggest that whereas size of some of the tests of hypotheses based on the considered generalized distributions are essentially ${\alpha}$-level, some are possibly not; power of tests of hypotheses on the Lindley distribution parameter from the two distributions differs.

• Journal of the Korean Mathematical Society
• /
• v.54 no.2
• /
• pp.479-491
• /
• 2017

We investigate the zeros of Epstein zeta functions associated with positive definite quadratic forms with rational coefficients in the vertical strip ${\sigma}_1$ < ${\Re}s$ < ${\sigma}_2$, where 1/2 < ${\sigma}_1$ < ${\sigma}_2$ < 1. When the class number h of the quadratic form is bigger than 1, Voronin gave a lower bound and Lee gave an asymptotic formula for the number of zeros. Recently Gonek and Lee improved their results by providing a new upper bound for the error term when h > 3. In this paper, we consider the cases h = 2, 3 and provide an upper bound for the error term, smaller than the one for the case h > 3.

• Bulletin of the Korean Mathematical Society
• /
• v.48 no.4
• /
• pp.713-724
• /
• 2011

We study the dynamics of transcendental entire functions with Siegel disks whose singular values are just two points. One of the two singular values is not only a superattracting fixed point with multiplicity more than two but also an asymptotic value. Another one is a critical value with free dynamics under iterations. We prove that if the multiplicity of the superattracting fixed point is large enough, then the restriction of the transcendental entire function near the Siegel point is a quadratic-like map. Therefore the Siegel disk and its boundary correspond to those of some quadratic polynomial at the level of quasiconformality. As its applications, the logarithmic lift of the above transcendental entire function has a wandering domain whose shape looks like a Siegel disk of a quadratic polynomial.

• Bulletin of the Korean Mathematical Society
• /
• v.54 no.1
• /
• pp.343-358
• /
• 2017

Let {$X_i$} be a sequence of stationary ${\alpha}-mixing$ random variables with probability density function f(x). The recursive kernel estimators of f(x) are defined by $$\hat{f}_n(x)={\frac{1}{n\sqrt{b_n}}{\sum_{j=1}^{n}}b_j{^{-\frac{1}{2}}K(\frac{x-X_j}{b_j})\;and\;{\tilde{f}}_n(x)={\frac{1}{n}}{\sum_{j=1}^{n}}{\frac{1}{b_j}}K(\frac{x-X_j}{b_j})$$, where 0 < $b_n{\rightarrow}0$ is bandwith and K is some kernel function. Under appropriate conditions, we establish the Berry-Esseen bounds for these estimators of f(x), which show the convergence rates of asymptotic normality of the estimators.