• Journal of Korean Society for Quality Management
• /
• v.23 no.4
• /
• pp.28-41
• /
• 1995

We consider the problem of identifying multiple outliers in linear model. The available regression diagnostic methods often do not succeed in detecting multiple outliers because of the masking and swamping effect. Recently, among the various robust estimator of reducing the effect of outliers, LMS(Least Meadian Square) estimator has been to be a suitable method proposed to expose outliers and leverage points. However, as you know it, the data analysis method with LMS estimator is to be taken the median of the squared residuals in the sample which is extracted the sample space. Then this model causes the trouble, for the number of the chosen sample is nCp, i.e. as the size of sample space n is increasing, the number is increasing fastly. And the covariance matrix may be the singular matrix, so that matrix is approching collinearity. Thus we propose a procedure ELMS for the resampling in LMS method and study the size of the effective elementary set in this algorithm.

• The Journal of Asian Finance, Economics and Business
• /
• v.7 no.8
• /
• pp.161-173
• /
• 2020

This study examines the nonlinear relationship between financial development and economic growth in Pakistan using the threshold regression model for the period 1980-2017. We also employed quantile regression with 0.25, 0.50, and 0.75 quantiles of conditional distribution. The quantile regression is based on minimizing of sum of squared residuals. The result indicates that economic growth responds positively to financial development when the level of financial development surpasses the threshold value of 0.151. However, when financial development lies below the threshold value (that is, 0.151), its impact on economic growth is negative. Thus, when financial development of Pakistan surpasses the threshold level, it contributes more towards economic growth since greater level of financial development contributes more to boosts economic growth. This finding reveals that economic growth reacts differently to financial development, and the relationship between financial development and economic growth is U-shaped in Pakistan. Among the other variables, physical capital, labor force, and government expenditure exert a positive effect on economic growth. Furthermore, inflation rate and trade openness have an insignificant impact on economic growth. The results of quantile regression also confirm the non-linear relationship between financial development and economic growth in Pakistan. The finding of this study suggests revamping of financial sector policies in Pakistan.

• Journal of Korea Port Economic Association
• /
• v.31 no.1
• /
• pp.1-14
• /
• 2015

The standard GARCH model imposing symmetry on the conditional variance, tends to fail in capturing some important features of the data. This paper, hence, introduces the models capturing asymmetric effect. They are the EGARCH model and the GJR model. We provide the systematic comparison of volatility models focusing on the asymmetric effect of news on volatility. Specifically, three diagnostic tests are provided: the sign bias test, the negative size bias test, and the positive size bias test. This paper shows that there is significant evidence of GARCH-type process in the data, as shown by the test for the Ljung-Box Q statistic on the squared residual data. The estimated unconditional density function for squared residual is clearly skewed to the left and markedly leptokurtic when compared with the standard normal distribution. The observation of volatility clustering is also clearly reinforced by the plot of the squared value of residuals of export volume and values. The unconditional variance of both export volumes and export value indicates that large shocks of either sign tend to be followed by large shocks, and small shocks of either sign tend to follow small shocks. The estimated export volume news impact curve for the GARCH also suggests that $h_t$ is overestimated for large negative and positive shocks. The conditional variance equation of the GARCH model for export volumes contains two parameters ${\alpha}$ and ${\beta}$ that are insignificant, indicating that the GARCH model is a poor characterization of the conditional variance of export volumes. The conditional variance equation of the EGARCH model for export value, however, shows a positive sign of parameter ${\delta}$, which is contrary to our expectation, while the GJR model exhibits that parameters ${\alpha}$ and ${\beta}$ are insignificant, and ${\delta}$ is marginally significant. That indicates that the asymmetric volatility models are poor characterization of the conditional variance of export value. It is concluded that the asymmetric EGARCH and GJR model are appropriate in explaining the volatility of export volume, while the symmetric standard GARCH model is good for capturing the volatility.

• Journal of Industrial Technology
• /
• v.36
• /
• pp.71-76
• /
• 2016

In this study, we determined the precise coordinates of TBMs (Tidal Bench Marks), which used as the national reference points in coastal area of Korea, using a GPS data analysis SW for the academic and scientific applications, GAMIT/GLOBK. For accurate 3-D positioning of TBM locations, we performed the GPS point surveying according to the national surveying policy and also acquired the GPS data for 48 TBMs located in the western and southern coastal part of Korea. Considering the results of baseline analysis to each observation session obtained from GAMIT module, the baseline analysis was realized to be done precisely because the values of Normalized RMS (NRMS) were mostly less than ${\pm}0.20mm$. Before the network adjustment using GLOBK module, we evaluated the suitability of observations for each session by applying the chi-squared test (${\chi}^2$ test) to the degree of freedom in observed session. An overall distributions of ${\chi}^2$ test were less than 1.0 for all sessions, and the statistical of ${\chi}^2$ test showed the average, 0.267 with minimum and maximum value, 0.063 and 0.653, respectively. Finally, we analyzed the network adjustment for 48 TBMs to reduce the residuals of baseline analysis on each point by connecting with 42 permanent GPS stations in Korea. In the network adjustment procedure, we set up the weighted values of each permanent station to be allocated between 0.9 and 1.14, and also removed the observed points having residual exceeds 4-times of standard deviation ($4{\sigma}$).

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.35 no.1
• /
• pp.55-62
• /
• 2017

The SRP (Solar Radiation Pressure) model has always been an issue in the dynamic GPS (Global Positioning System) orbit determination. The widely used CODE (Center for Orbit Determination in Europe) model and its variants have nine parameters to estimate the solar radiation pressure from the Sun and to absorb the remaining forces. However, these parameters show a very high correlation with each other and, therefore, only several of them are estimated at most of the IGS (International GNSS Service) analysis centers. In this study, we attempted to numerically verify the correlation between the parameters. For this purpose, a bi-directional, multi-step numerical integrator was developed. The correlation between the SRP parameters was analyzed in terms of post-fit residuals of the orbit. The integrated orbit was fitted to the IGS final orbit as external observations. On top of the parametric analysis of the SRP parameters, we also verified the capabilities of orbit prediction at later time epochs. As a secondary criterion for orbit quality, the positional discontinuity of the daily arcs was also analyzed. The resulting post-fit RMSE (Root-Mean-Squared Error) shows a level of 4.8 mm on average and there is no significant difference between block types. Since the once-per-revolution parameters in the Y-axis are highly correlated with those in the B-axis, the periodic terms in the D- and Y-axis are constrained to zero in order to resolve the correlations. The 6-hr predicted orbit based on the previous day yields about 3 cm or less compared to the IGS final orbit for a week, and reaches up to 6 cm for 24 hours (except for one day). The mean positional discontinuity at the boundary of two 1-day arcs is on the level of 1.4 cm for all non-eclipsing satellites. The developed orbit integrator shows a high performance in statistics of RMSE and positional discontinuity, as well as the separations of the dynamic parameters. In further research, additional verification of the reference frame for the estimated orbit using SLR is necessary to confirm the consistency of the orbit frames.

• Journal of the Korean Data and Information Science Society
• /
• v.20 no.2
• /
• pp.293-299
• /
• 2009

In this article, we consider linear models such as regression, ARIMA (autoregressive integrated moving average), and regression+ARIMA (regression with ARIMA errors) for predicting hourly ozone concentration level in two areas of Daegu. Based on RASE(root average squared error), it is shown that the ARIMA is the best model in one area and that the regression+ARIMA model is the best in the other area. We further analyze the residuals from the optimal models, so that we might predict the ozone warning days where at least one of the hourly ozone concentration levels is over 120 ppb. Based on the training data in the years from 2000 to 2003, it is found that 35 ppb is a good cutoff value of residulas for predicting the ozone warning days. In on area of Daegu, our method predicts correctly one of two ozone warning days of 2004 as well as all of the remaining 364 non-warning days. In the other area, our methods predicts correctly all of one ozone warning days and 365 non-warning days of 2004.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.35 no.3
• /
• pp.216-222
• /
• 2002

Autocorrelation in time series data can affect statistical inference in correlation or regression analyses. To improve a regression model from which the residuals are autocorrelated, Yule-Walker method, nonlinear least squares estimation, maximum likelihood method and 'prewhitening' method have been used to estimate the parameters in a regression equation. This study reviewed on the estimation methods of preventing spurious correlation in the presence of autocorrelation and applied the former three methods, Yule-Walker, nonlinear least squares and maximum likelihood method, to a 20-year real data set. Monte carlo simulation was used to compare the three parameter estimation methods. However, the simulation results showed that the mean squared error distributions from the three methods simulated do not differ significantly.

• The Korean Journal of Applied Statistics
• /
• v.29 no.3
• /
• pp.471-485
• /
• 2016

Least squares (LS) regression is a classic method for regression that is optimal under assumptions of regression and usual observations. However, the presence of unusual data in the LS method leads to seriously distorted estimates. Therefore, various robust estimation methods are proposed to circumvent the limitations of traditional LS regression. Among these, there are M-estimators based on maximum likelihood estimation (MLE), L-estimators based on linear combinations of order statistics and R-estimators based on a linear combinations of the ordered residuals. In this paper, robust regression estimators with high breakdown point and/or with high efficiency are compared under several simulated situations. The paper analyses and compares distributions of estimates as well as relative efficiencies calculated from mean squared errors (MSE) in the simulation study. We conclude that MM-estimators or GR-estimators are a good choice for the real data application.

• The Journal of Fisheries Business Administration
• /
• v.42 no.2
• /
• pp.69-83
• /
• 2011

This paper estimates demand functions of oyster as Kimchi's ingredients of capital area, other areas excluding a capital area, and a whole area in Korea to forecast its demand quantities in 2011~2015. To estimate oyster demand function, this paper uses pooled data produced from Korean housewives over 30 years old in 2009 and 2010. Also, this paper adopts several econometrics methods such as Ordinary Least Squares and Feasible Generalized Least Squares. First of all, to choose appropriate variables of oyster demand functions by area, this paper carries out model's specification with joint significance test. Secondly, to remedy heteroscedasticity with pooled data, this paper attempts residual plotting between estimated squared residuals and estimated dependent variable and then, if it happens, undertakes White test to care the problem. Thirdly, to test multicollinearity between variables with pooled data, this paper checks correlations between variables by area. In this analysis, oyster demand functions of a capital area and a whole area need price of the oyster, price of the cabbage for Gimjang, and income as independent variables. The function on other areas excluding a capital area only needs price of the oyster and income as ones. In addition, the oyster demand function of a whole area needed White test to care a heteroscedasticity problem and demand functions of the other two regions did not have the problem. Thus, first model was estimated by FGLS and second two models were carried out by OLS. The results suggest that oyster demand quantities per a household as Kimchi's ingredients are going to slightly increase in a capital area and a whole area, but slightly decrease in other areas excluding a capital area in 2011~2015. Also, the results show that oyster demand quantities as kimchi's ingredients for total household targeting housewives over 30 years old are going to slightly increase in three areas in 2011~2015.

• Animal Bioscience
• /
• v.35 no.5
• /
• pp.648-658
• /
• 2022

Objective: The identification of nonlinear mixed models that describe the growth trajectory of New Zealand rabbits was performed based on weight records and carcass measures obtained using ultrasonography. Methods: Phenotypic records of body weight (BW) and loin eye area (LEA) were collected from 66 animals raised in a didactic-productive module of cuniculture located in the southern Piaui state, Brazil. The following nonlinear models were tested considering fixed parameters: Brody, Gompertz, Logistic, Richards, Meloun 1, modified Michaelis-Menten, Santana, and von Bertalanffy. The coefficient of determination (R²), mean squared error, percentage of convergence of each model (%C), mean absolute deviation of residuals, Akaike information criterion (AIC), and Bayesian information criterion (BIC) were used to determine the best model. The model that best described the growth trajectory for each trait was also used under the context of mixed models, considering two parameters that admit biological interpretation (A and k) with random effects. Results: The von Bertalanffy model was the best fitting model for BW according to the highest value of R² (0.98) and lowest values of AIC (6,675.30) and BIC (6,691.90). For LEA, the Logistic model was the most appropriate due to the results of R² (0.52), AIC (783.90), and BIC (798.40) obtained using this model. The absolute growth rates estimated using the von Bertalanffy and Logistic models for BW and LEA were 21.51g/d and 3.16 cm², respectively. The relative growth rates at the inflection point were 0.028 for BW (von Bertalanffy) and 0.014 for LEA (Logistic). Conclusion: The von Bertalanffy and Logistic models with random effect at the asymptotic weight are recommended for analysis of ponderal and carcass growth trajectories in New Zealand rabbits. The inclusion of random effects in the asymptotic weight and maturity rate improves the quality of fit in comparison to fixed models.