• Journal of Applied and Pure Mathematics
• /
• v.4 no.5_6
• /
• pp.263-268
• /
• 2022

We derive an alternative proof of Marsaglia's method for generating a pair of independent standard normal random variables.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.8 s.179
• /
• pp.2074-2082
• /
• 2000

In this study, to investigate and to predict the crack growth behavior under variable amplitude loading, crack growth tests are conducted on 7075-T6 aluminum alloy. The loading wave forms are generated by normal random number generator. All wave forms have same average and RMS(root mean square) value, but different standard deviation, which is to vary the maximum load in each wave. The modified Forman's equation is used as crack growth equation. Using the retardation coefficient D defined in previous study, the load interaction effect is considered. The variability in crack growth process is described by the random variable Z which was obtained from crack growth tests under constant amplitude loading in previous work. From these, a statistical model is developed. The curves predicted by the proposed model well describe the crack growth behavior under variable amplitude loading and agree with experimental data. In addition, this model well predicts the variability in crack growth process under variable amplitude loading.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.1162-1168
• /
• 2008

This study provides how the Dimension Reduction (DR) method as an efficient technique for reliability analysis can acquire its increased efficiency when it is applied to highly nonlinear problems. In the highly nonlinear engineering systems, 4N+1 (N: number of random variables) sampling is generally recognized to be appropriate. However, there exists uncertainty concerning the standard for judgment of non-linearity of the system as well as possibility of diverse degrees of non-linearity according to each of the random variables. In this regard, this study judged the linearity individually on each random variable after 2N+1 sampling. If high non-linearity appeared, 2 additional sampling was administered on each random variable to apply the DR method. The applications of the proposed sampling to the examples produced the constant results with increased efficiency.

• Journal of the Korean Statistical Society
• /
• v.30 no.1
• /
• pp.31-39
• /
• 2001

For estimating parameters of possibly nonlinear and/or non-stationary seasonal autoregressive(AR) processes, we introduce a new instrumental variable method which use the direction vector of the regressors in the same period as an instrument. On the basis of the new estimator, we propose new seasonal random walk tests whose limiting null distributions are standard normal regardless of the period of seasonality and types of mean adjustments. Monte-Carlo simulation shows that he powers of he proposed tests are better than those of the tests based on ordinary least squares estimator(OLSE).

• Communications for Statistical Applications and Methods
• /
• v.26 no.4
• /
• pp.371-383
• /
• 2019

Panel data sets have been developed in various areas, and many recent studies have analyzed panel, or longitudinal data sets. Maximum likelihood (ML) may be the most common statistical method for analyzing panel data models; however, the inference based on the ML estimate will have an inflated Type I error because the ML method tends to give a downwardly biased estimate of variance components when the sample size is small. The under estimation could be severe when data is incomplete. This paper proposes the restricted maximum likelihood (REML) method for a random effects panel data model with a censored dependent variable. Note that the likelihood function of the model is complex in that it includes a multidimensional integral. Many authors proposed to use integral approximation methods for the computation of likelihood function; however, it is well known that integral approximation methods are inadequate for high dimensional integrals in practice. This paper introduces to use the moments of truncated multivariate normal random vector for the calculation of multidimensional integral. In addition, a proper asymptotic standard error of REML estimate is given.

• Communications for Statistical Applications and Methods
• /
• v.16 no.5
• /
• pp.841-849
• /
• 2009

Let {$X_n$, n ${\ge}$ 1} be a negatively associated sequence of identically distributed random variables with mean zeros and positive finite variances. Set $S_n$ = ${\Sigma}^n_{i=1}\;X_i$. Suppose that 0 < ${\sigma}^2=EX^2_1+2{\Sigma}^{\infty}_{i=2}\;Cov(X_1,\;X_i)$ < ${\infty}$. We prove that, if $EX^2_1(log^+{\mid}X_1{\mid})^{\delta}$ < ${\infty}$ for any 0< ${\delta}{\le}1$, then $\lim_{{\epsilon}\downarrow0}{\epsilon}^{2{\delta}}\sum_{{n=2}}^{\infty}\frac{(logn)^{\delta-1}}{n^2}ES^2_nI({\mid}S_n{\mid}\geq{\epsilon}{\sigma}\sqrt{nlogn}=\frac{E{\mid}N{\mid}^{2\delta+2}}{\delta}$, where N is the standard normal random variable. We also prove that if $S_n$ is replaced by $M_n=max_{1{\le}k{\le}n}{\mid}S_k{\mid}$ then the precise rate still holds. Some results in Fu and Zhang (2007) are improved to the complete moment case.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2004.03b
• /
• pp.319-326
• /
• 2004

Uncertainty is inevitably involved in rock slope engineering since the rock masses are formed by natural process and subsequently the geotechnical characteristics of rock masses cannot be exactly obtained. Therefore the reliability analysis method has been suggested to deal properly with uncertainty. The reliability analysis method can be divided into level I, II and III on the basis of the approach for consideration of random variable and probability density function of reliability function. The level II approach, which is focused in this study, assumes the probability density function of random variables as normal distribution and evaluates the probability of failure with statistical moments such as mean and standard deviation. This method has the advantage that can be used the problem which the Monte Carlo simulation approach cannot be applied since the complete information on the random variables are not available. In this study, the analysis results of level II reliability approach compared with the analysis results of level III approach to verify the appropriateness of the level II approach. In addition, the results are compared with the results of the deterministic analysis.

• The Korean Journal of Applied Statistics
• /
• v.26 no.1
• /
• pp.37-47
• /
• 2013

A Q-Q plot is an effective and convenient graphical method to assess a distributional assumption of data. The primary step in the construction of a Q-Q plot is to obtain a closed-form expression to represent the relation between observed quantiles and theoretical quantiles to be plotted in order that the points fall near the line y = a + bx. In this paper, we introduce a Q-Q plot to assess goodness-of-fit for inverse Gaussian distribution. The procedure is based on the distributional result that a transformed random variable $Y={\mid}\sqrt{\lambda}(X-{\mu})/{\mu}\sqrt{X}{\mid}$ follows a half-normal distribution with mean 0 and variance 1 when a random variable X has an inverse Gaussian distribution with location parameter ${\mu}$ and scale parameter ${\lambda}$. Simulations are performed to provide a guideline to interpret the pattern of points on the proposed inverse Gaussian Q-Q plot. An illustrative example is provided to show the usefulness of the inverse Gaussian Q-Q plot.

• Journal of the Korean Mathematical Society
• /
• v.47 no.2
• /
• pp.263-275
• /
• 2010

Let {$X_n;n\;\geq\;1$} be a strictly stationary sequence of negatively associated random variables with mean zero and finite variance. Set $S_n\;=\;{\sum}^n_{k=1}X_k$, $M_n\;=\;max_{k{\leq}n}|S_k|$, $n\;{\geq}\;1$. Suppose $\sigma^2\;=\;EX^2_1+2{\sum}^\infty_{k=2}EX_1X_k$ (0 < $\sigma$ < $\infty$). We prove that for any b > -1/2, if $E|X|^{2+\delta}$(0<$\delta$$\leq$1), then $$lim\limits_{\varepsilon\searrow0}\varepsilon^{2b+1}\sum^{\infty}_{n=1}\frac{(loglogn)^{b-1/2}}{n^{3/2}logn}E\{M_n-\sigma\varepsilon\sqrt{2nloglogn}\}_+=\frac{2^{-1/2-b}{\sigma}E|N|^{2(b+1)}}{(b+1)(2b+1)}\sum^{\infty}_{k=0}\frac{(-1)^k}{(2k+1)^{2(b+1)}}$$ and for any b > -1/2, $$lim\limits_{\varepsilon\nearrow\infty}\varepsilon^{-2(b+1)}\sum^{\infty}_{n=1}\frac{(loglogn)^b}{n^{3/2}logn}E\{\sigma\varepsilon\sqrt{\frac{\pi^2n}{8loglogn}}-M_n\}_+=\frac{\Gamma(b+1/2)}{\sqrt{2}(b+1)}\sum^{\infty}_{k=0}\frac{(-1)^k}{(2k+1)^{2b+2'}}$$, where $\Gamma(\cdot)$ is the Gamma function and N stands for the standard normal random variable.

• Geomechanics and Engineering
• /
• v.24 no.2
• /
• pp.167-178
• /
• 2021

Spatial variability is an inherent uncertainty of soil properties. Current reliability analyses generally incorporate random field theory and Monte Carlo simulation (MCS) when dealing with spatial variability, in which the computational efficiency is a significant challenge. This paper proposes a KL-FORM algorithm to improve the computational efficiency. In the proposed KL-FORM, Karhunen-Loeve (KL) expansion is used for discretizing random fields, and first-order reliability method (FORM) is employed for reliability analysis. The KL expansion and FORM can be used in conjunction, through adopting independent standard normal variables in the discretization of KL expansion as the basic variables in the FORM. To illustrate the effectiveness of this KL-FORM, it is applied to a case study of a strip footing in spatially variable unsaturated soil under rainfall, in which the bearing capacity of the footing is computed by numerical simulation. This case study shows that the KL-FORM is accurate and efficient. The parametric analyses suggest that ignoring the spatial variability of the soil may lead to an underestimation of the reliability index of the footing.