Mathematical Researches

fa مقایسۀ کوچک نمونه‌ای برآورد تابع چگالی‌ با روش‌های هستۀ گاما و سری‌های متعامد Small-Sample Comparison of the Gamma Kernel and the Orthogonal Series Methods of Density Estimation آمار stat مقاله مستقل Original Manuscript برآوردگر هسته‌ای استاندارد دچار مشکل اریبی مرزی برای تابع‌های چگالی احتمال توزیع‌هایی با تکیه‌گاه اعداد حقیقی مثبت است. برآوردگرهای هسته‌ای گاما و سری‌های متعامد دو جای‌گزین برای برآوردگر هسته‌ای استاندارد هستند که &rlm;مبرا از اریبی مرز&lrm;&rlm;ی‌اند&lrm;. در این مقاله یک بررسی شبیه‌سازی به‌منظور مقایسۀ عملکرد کوچک‌نمونه‌ای برآوردگرهای هسته‌ای گاما و برآوردگر سری‌های متعامد با پایه‌های لاگِر و اِرمیت اجرا شده است. بر اساس این مطالعۀ شبیه‌سازی دریافتیم که اگر انتخاب پایه‌ برای برآوردگر سری‌های متعامد بر اساس شناخت جزئی از شکل تابع چگالی هدف باشد، برآوردگر سری‌های متعامد می‌تواند عملکردی بهتر از برآوردگر هسته‌ای گاما داشته باشد. &lrm; Introduction Estimation of a probability density function is an important area of nonparametric statistical inference that has received much attention in recent decades. The kernel method is widely used in nonparametric estimation of the probability density function of an absolutely continuous distribution with support on the whole real line. However, for a distribution with support on a subset of the real line, the kernel density estimator with fixed symmetric kernels encounters bias at the boundaries of the support, which is known as the boundary bias issue. This is due to smoothing data near the boundary points by the fixed symmetric kernel that leads to allocating probability density to outside of the distribution’s support (see Silverman, 1986). There are many applications, such as reliability, insurance and life testing, dealing with non-negative data and estimating the probability density function of distributions with support on the non-negative real line is the object of interest.  Using the kernel estimator with fixed symmetric kernels in these cases results in the boundary bias issue at the origin. A number of methods have been proposed to avoid the boundary bias issue at the origin. A simple remedy is to replace symmetric kernels by asymmetric kernels which never assign density to negative values. The Gamma kernels proposed by Chen (2000) are the effective asymmetric kernels to estimate the probability density function of distributions on the non-negative real line.<a name="__DdeLink__179_2839299288"></a> Orthogonal series estimators form another class of nonparametric probability density estimators, which go back to Cencov (1964). In this approach, as reviewed in Efromovich (2010), the target probability function is expanded in terms of a sequence of orthogonal basis functions. After selecting a suitable sequence of orthogonal basis functions, the observed data are used to estimated <a name="__DdeLink__190_2839299288">the coefficients of the expansion</a> in order to obtain the orthogonal series density estimator. <a name="__DdeLink__3076_2839299288">Similar to kernel estimators</a>, <a name="__DdeLink__3078_2839299288">under some mild conditions </a>the orthogonal series estimators have appealing large sample properties. Moreover, the boundary issue can be avoided by using orthogonal density estimators with suitable basis functions. Although small sample properties of asymmetric kernel estimators with the Gamma kernels and orthogonal series estimators <a name="__DdeLink__3080_2839299288">are well-studied</a> separately, but to the best of our knowledge, there have been no reports of comparing their performance in estimating the probability density function of distributions on the non-negative real line. In this paper&lrm;, a simulation study is conducted to compare the small-sample performance of the Gamma kernel estimators and orthogonal series estimators for a set of distributions on the positive real line.   Material and methods Following Malec and Schienle (2014), we consider six parameter settings for the generalized F distribution to obtain probability density functions with different shapes, near-origin behaviors and tail decays (Figure 2). Based on 5000 simulations from any of these density functions with sample sizes <img alt="" height="19" src="file:///C:Users1AppDataLocalTempmsohtmlclip11clip_image002.gif" width="104" >, we estimate the target density function using the type I and II Gamma kernel estimators and the orthogonal series estimators with Hermite and Laguerre basis functions and compute the <a name="__DdeLink__3000_239705243">mean integrated squared error (MISE)</a>. The bandwidth parameter in the Gamma kernel estimators and the cutoff and smoothing parameters in the orthogonal series estimators are significantly affect the performance of the estimators. We use optimal<a name="__DdeLink__3004_239705243"> bandwidths</a> for the Gamma kernel estimators and <a name="__DdeLink__2998_239705243">optimal cutoff and smoothing</a><a name="__DdeLink__3002_239705243"> parameters</a> of the orthogonal series estimators to avoid variations due to uncertainty of tuning parameters. To obtain the optimal tuning parameters for each target density, we compute and minimize the MISE with respect to the tuning parameters based on additional 5000 simulations from the true density function. Results and discussion  For each density function, the optimal tuning parameters and the MISE’s of the estimators are reported (Table 2). As expected from the large sample properties, increasing the sample size improves the performance of all estimators. The performances of estimators vary from cases to cases and there no considered estimator is the best in all cases. In all cases except one, the type II Gamma kernel estimator is superior to the type I Gamma kernel estimator, which is in agreement with Chen (2000) <a name="__DdeLink__2923_3606392128">suggestion of preferring </a>type II to type I Gamma kernel estimator. However, in <a name="_GoBack">one case the type I Gamma kernel estimator is better than all other estimators. In cases where the </a>shape and near-origin behavior of the target density is similar to the Hermit or Laguerre basis functions, the corresponding orthogonal series estimator outperforms all the other competing estimators. Conclusion The following conclusions were drawn from this research <a name="__DdeLink__2927_3606392128">to choose among the considered</a> estimators. <ul> <li>If the basis functions of the orthogonal series estimator <a name="__DdeLink__2925_3606392128">are chosen to have similar shape</a> and near-origin behavior to the target density function, then the corresponding orthogonal series estimator can outperform the Gamma kernel estimators.</li> <li>If there is no prior knowledge about the shape and near-origin behavior of the target density function and the sample size is relatively large (n=400),  then the type II Gamma kernel estimator can outperform the orthogonal series estimators.</li> </ul>     اریبی مرزی, میانگین جمع‌بستۀ توان دوم خطا, پارامتر هموارساز. ‎ کد موضوع‌بندی ریاضی (2010): 62G07, 62Gxx, 93E14 Boundary bias, Generalized F distribution, Mean integrated squared error, Smoothing parameter. 333 346 http://mmr.khu.ac.ir/browse.php?a_code=A-10-377-1&slc_lang=fa&sid=1 Muhyiddin Izadi محی الدین ایزدی m.izadi@razi.ac.ir 10031947532846003948 10031947532846003948 Yes Razi University دانشگاه رازی، گروه آمار Abdollah Jalilian عبدلاه جلیلیان jalilian@razi.ac.ir 10031947532846003949 10031947532846003949 No Razi University دانشگاه رازی، گروه آمار