<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJS</journal-id><journal-title-group><journal-title>Open Journal of Statistics</journal-title></journal-title-group><issn pub-type="epub">2161-718X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojs.2016.65060</article-id><article-id pub-id-type="publisher-id">OJS-70799</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Sample Selection Model with Bootstrap (BPSSM) Approach: Case Study of the Malaysian Population and Family Survey
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Muhamad</surname><given-names>Safiih Lola</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Wan</surname><given-names>Saliha Wan Alwi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Nurul</surname><given-names>Hila Zainuddin</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>School of Informatics and Applied Mathematics, Universiti Malaysia Terengganu, Kuala Terengganu, Terengganu, Malaysia</addr-line></aff><pub-date pub-type="epub"><day>22</day><month>09</month><year>2016</year></pub-date><volume>06</volume><issue>05</issue><fpage>741</fpage><lpage>748</lpage><history><date date-type="received"><day>June</day>	<month>30,</month>	<year>2016</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>September</month>	<year>19,</year>	</date><date date-type="accepted"><day>September</day>	<month>22,</month>	<year>2016</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Heckman Sampel Selection Model (PSSM) has been adopted widely in the study of labour work. This model contains exogenous, endogenous and standard error variables. However, this model is constantly exposed to high inaccuracy of estimation result. Therefore, to obtain an accurate and precise estimation, the bootstrap approach is introduced. The bootstrap approach will be hybrid with PSSM model known as BPSSM to achieve estimation result that is more precise. Then, the BPSSM is applied to Malaysian Population and Family Survey 1994 (MPFS-1994) data. The results showed that BPSSM provide a smaller standard error and shorter confidence intervals.
 
</p></abstract><kwd-group><kwd>Sampel Selection</kwd><kwd> Bootstrap</kwd><kwd> Standard Error</kwd><kwd> Confidence Intervals</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Sample selection model is part of the field of econometrics. The term “selection” or “select” is the term commonly used and it is mentioned in a number of different issues related with the econometric data. Sample selection was developed [<xref ref-type="bibr" rid="scirp.70799-ref1">1</xref>] . This model also has good interaction of relation and ideal to expert evaluation and quantitative information.</p><p>The earliest introduced model consists of female labour force and wage equality [<xref ref-type="bibr" rid="scirp.70799-ref2">2</xref>] , [<xref ref-type="bibr" rid="scirp.70799-ref3">3</xref>] . They have used this model widely in the model affects unions, occupational choices, schooling choices, options and also a residential area of the industry. According to [<xref ref-type="bibr" rid="scirp.70799-ref4">4</xref>] in the female labour supply model, the binary model choice whether a person’s functioning and conditions of work can be seen in working hours. Traditionally, this model has been estimated by maximum likelihood or using computerized methods called Heckman two-step method.</p><p>Model selection consists of two parts [<xref ref-type="bibr" rid="scirp.70799-ref5">5</xref>] . The first part is the selection equation (known as the participation or equality of results). In this section, not random sample taken into account for the observation of the process of entering the sample and determine a probability sample of the population. This combination can fix the not random sample and the estimated relation of the population. The second part is the structural equation. This is also known as equality of outcome or equality of salaries with relationship desired study population centres.</p><p>Therefore, to get the accurate and precise estimation, bootstrap approach was introduced. Through this method, Bootstrap approach will be hybridized with PSSM model called Bootstrap PSSM to get more accurate estimator result. The best model is the model contains of consistent and efficient.</p><p>Bootstrap introduced by [<xref ref-type="bibr" rid="scirp.70799-ref6">6</xref>] . Bootstrap is a common technique that builds confidence interval by resampling with replacement sample from finite. The main objective of this article, we propose to develop a sample selection model (PSSM) based on bootstrap approach to achieve estimation result that is more precise. After that, the model of BPSSM was applied for the Family and Population Data 1994.</p></sec><sec id="s2"><title>2. Methods</title>Bootstrapping the Base Model of BPSSM<p>In this study, a mean process is considered to monitor individual observations of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x2.png" xlink:type="simple"/></inline-formula> with assumption of dependent and uncorrelated distributed. Thus, the base model and residual of BPSSM can be given by:</p><disp-formula id="scirp.70799-formula2"><graphic  xlink:href="http://html.scirp.org/file/1-1240723x3.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.70799-formula3"><graphic  xlink:href="http://html.scirp.org/file/1-1240723x4.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.70799-formula4"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1240723x5.png"  xlink:type="simple"/></disp-formula><p>where</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x6.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x7.png" xlink:type="simple"/></inline-formula> = dependent variables</p><p>x and z = vectors of exogenous remaining variables</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x8.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x9.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x10.png" xlink:type="simple"/></inline-formula> = unknown parameter vectors</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x11.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x12.png" xlink:type="simple"/></inline-formula> = zero mean error terms</p><p>The standard approach is to assume that follow a bivarite normal distribution and then applied to the maximum likelihood estimation or a two-stage estimation procedure purposed by Heckman (1979). Firstly, how to estimate <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x13.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x14.png" xlink:type="simple"/></inline-formula> consistently from the data. In general, both the error terms are correlated, since that the regression of y on x for the selected sample will not give consistent estimates of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x15.png" xlink:type="simple"/></inline-formula>. It is well known that the consistency of those estimators depends on the assumption of bivariate normality. For a random sample from the population it is observed that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x16.png" xlink:type="simple"/></inline-formula> x and z. If and only if, observation of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x17.png" xlink:type="simple"/></inline-formula> then, we observed<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x18.png" xlink:type="simple"/></inline-formula>. This sample selection models in (1) consist of two equations or parts; the first structural part, embodying the desired population relationship or is the equation of primary interest and second, the selection part or is the reduced form takes account of the non-representative nature of the present non-random sample. From the literature (Martin, 2001), the identification purpose, the variable z contains at least one variable which does not appear the relation between an outcome in interest y and a vector of covariates x and the selection equation describing the relation between a binary participation decision <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x19.png" xlink:type="simple"/></inline-formula> and another vector of covariates z.</p><p>In this study, the hybridization of bootstrap approach in base model (1). This hybridization produced a hybrid control charts where the basic model named by Bootstrap PSSM (BPSSM). The algorithm for this hybrid process is as follows:</p><p>Step 1: A sample data, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x20.png" xlink:type="simple"/></inline-formula>generated by time <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x21.png" xlink:type="simple"/></inline-formula> from a dependent and colerated distribution. The data then will be used to estimate parameter of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x22.png" xlink:type="simple"/></inline-formula>, base model of PSSM <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x23.png" xlink:type="simple"/></inline-formula></p><p>Step 2: Find bootstrap replication, B(c) by using:</p><disp-formula id="scirp.70799-formula5"><graphic  xlink:href="http://html.scirp.org/file/1-1240723x24.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x25.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x26.png" xlink:type="simple"/></inline-formula> defines as average of expected of error and true value of error respectively. This estimation will be repeated several times to get a constant value of bias and this new value will be used to analyze data in Step 1.</p><p>Step 3: By continue from Step 2, the residual value will be used in sampling with replacement method to get a matrix of residual bootstrap, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x27.png" xlink:type="simple"/></inline-formula></p><disp-formula id="scirp.70799-formula6"><graphic  xlink:href="http://html.scirp.org/file/1-1240723x28.png"  xlink:type="simple"/></disp-formula><p>Step 4: For each residual in Step 3, compute new data <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x29.png" xlink:type="simple"/></inline-formula></p><p>Step 5: Compute average of column of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x30.png" xlink:type="simple"/></inline-formula> by using<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x31.png" xlink:type="simple"/></inline-formula>, with c defines as summation of bootstrap replication, for example c = 1000</p><p>Step 6: By using bootstrap data, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x32.png" xlink:type="simple"/></inline-formula>compute parameter of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x33.png" xlink:type="simple"/></inline-formula></p><p>Step 7: For complete PSSM model, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x34.png" xlink:type="simple"/></inline-formula>must be hybrid with bootstrap approach. Step 1 until Step 6 rapidly to compute bootstrap data of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x35.png" xlink:type="simple"/></inline-formula>.</p><p>Step 8: After the both equation have bootstrap data, Heckman Two Step estimation was used base on model BPSSM;</p><disp-formula id="scirp.70799-formula7"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1240723x36.png"  xlink:type="simple"/></disp-formula><p>In this study, we intend to examine the performance of BPSSM in terms of effectiveness and efficiency control. Numerical estimation was selected to be used in this study. For numerical estimation, basically it is used to examine effectiveness of base model where is evident in two kind of methods, i.e. confidence interval, bootstrap percentile (PB) and Biased Corrected and Accelerate (BCa). BP and BCa selection is motivated by the advantages of these two methods in which BP is the basic method for estimating bootstrap intervals while BCa is a method that can improve BP interval estimation [<xref ref-type="bibr" rid="scirp.70799-ref7">7</xref>] . The main reason for selecting different methods is to finds the differences in the effectiveness of the hybrid model when using those interval methods. In theory, a model that gives the shortest interval estimation is said to be more effective model. This is because; short interval giving the idea that model estimation is closer to real interval estimation. Therefore, the lower and upper limit <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x37.png" xlink:type="simple"/></inline-formula> can be given by:</p><p>Student’s-t:</p><disp-formula id="scirp.70799-formula8"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1240723x38.png"  xlink:type="simple"/></disp-formula><p>Bootstrap Persentile (BP):</p><disp-formula id="scirp.70799-formula9"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1240723x39.png"  xlink:type="simple"/></disp-formula><p>Bias Corrected and accelerated (Bca):</p><disp-formula id="scirp.70799-formula10"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1240723x40.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x41.png" xlink:type="simple"/></inline-formula> in Equation (3) represent value of percentage <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x42.png" xlink:type="simple"/></inline-formula> for student’s-t distribution with n − 1 degree of freedom.</p><p>In this study, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x43.png" xlink:type="simple"/></inline-formula>is valued as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x44.png" xlink:type="simple"/></inline-formula> and standard error estimation, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x45.png" xlink:type="simple"/></inline-formula>can be calculated using the discussion on [<xref ref-type="bibr" rid="scirp.70799-ref7">7</xref>] .</p><p>For Equation (4), length of this interval based on percentile of mean estimation of bootstrap replication, B of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x46.png" xlink:type="simple"/></inline-formula> where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x47.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x48.png" xlink:type="simple"/></inline-formula> [<xref ref-type="bibr" rid="scirp.70799-ref7">7</xref>] . In other word, upper and lower limit of BP refers to the interval length from 50-th through 950-th replication. While <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x49.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x50.png" xlink:type="simple"/></inline-formula> in Equation (5) refers to normal confidence interval with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x51.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x52.png" xlink:type="simple"/></inline-formula> respectively.</p><p>Mean Square Error (MSE) and Root Mean Square Error (RMSE) was used for error estimation in this study. In theory, a model that gives the smallest estimation value is said to be more efficient and automatically show the effectiveness of the model itself [<xref ref-type="bibr" rid="scirp.70799-ref8">8</xref>] . By taking the idea of small error in base model, it clearly shows that model is giving more accurate estimation. Thus, error estimation used in this study is given:</p><disp-formula id="scirp.70799-formula11"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1240723x53.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.70799-formula12"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1240723x54.png"  xlink:type="simple"/></disp-formula><p>where for both MSE and RMSE refers to differences of real error, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x55.png" xlink:type="simple"/></inline-formula>with expected of error estimation,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1240723x56.png" xlink:type="simple"/></inline-formula>.</p></sec><sec id="s3"><title>3. Results and Discussion</title>Application of Hybrid Model: The Malaysian Population and Family Survey 1994 Data<p>A comparison of performance of the real model, PSSM and hybrid model, BPSSM in terms of effectiveness or efficiency base on model estimation. The data set used for this study is from the Malaysian population and family survey 1994 (MPFS-1994). This survey was conducted by National Population and Family Development Board of Malaysia under Ministry of Women, Family and Community Development Malaysia. This survey was specifically for married women, providing relevant and significant information for the problem of married women status regarding wages, educational attainment, household composition and other socioeconomic characteristics. The original MPFS-94 sample data comprises 4444 married women.</p><p>The whole data sets used in this study consisted of 2792 married women. The selection rules [<xref ref-type="bibr" rid="scirp.70799-ref9">9</xref>] were applied to create the sample criteria of choosing for participant and non-participant married women on the basis of the MPFS-94 data set, which are as follows:</p><p>・ Married and aged below 60</p><p>・ Not in school or retired</p><p>・ Husband present in 1994</p><p>・ Husband reported positive earnings for 1994</p><p>The empirical results of the basic specification one are presented for the Heckman two-step approach. These approaches consider the probit estimates for the participation equation as a first step and OLS estimates for the wage equation as the second step. We discuss both the participation and wage equation on the estimated coefficient for interval method, the significant effect, and consistency and for PSSM, as well as BPSSM for comparison purposes.</p><p>Based on <xref ref-type="table" rid="table1">Table 1</xref>, the estimate for the standard intervals (Student’s-t confidence interval) using hybrid model gives a shorter interval. See for example, the length of BPSSM-t valued 0.0512791 compare to real model, PSSM-t which is more length 1.0898210. This result can be seen clearly in <xref ref-type="fig" rid="fig1">Figure 1</xref>. The differences shown in these two models proved that the bootstrap approach fixed the interval estimation and gives a good performance for hybrid model.</p><p>In <xref ref-type="fig" rid="fig1">Figure 1</xref> and <xref ref-type="fig" rid="fig2">Figure 2</xref>, a plot of interval estimation for Student’s-t method in <xref ref-type="fig" rid="fig1">Figure 1</xref> and also a plot of Bootstrap Percentile (BP) and Bias Corrected and accelerated (BCa) showed in <xref ref-type="fig" rid="fig2">Figure 2</xref>. Based on these two plots, we found that bootstrap interval method, i.e. BP and BCa, gave short length compare to the standard interval which show a length interval either for real or hybrid model. Significant difference on the length value shows that the bootstrap percentile interval method gave a good performance compare to Student’s-t estimation. By considering this result, we can say that BPSSM model gives a better performance when using bootstrap interval method. With this result, it’s clearly proved that hybrid model provides an effective estimation compare to the real model.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Interval estimation for the real and hybrid model</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Model</th><th align="center" valign="middle" >Lower</th><th align="center" valign="middle" >Upper</th><th align="center" valign="middle" >Length</th></tr></thead><tr><td align="center" valign="middle" >PSSM-t</td><td align="center" valign="middle" >22.8134200</td><td align="center" valign="middle" >23.9032400</td><td align="center" valign="middle" >1.0898210</td></tr><tr><td align="center" valign="middle" >BPSSM-t</td><td align="center" valign="middle" >23.2661200</td><td align="center" valign="middle" >23.3284100</td><td align="center" valign="middle" >0.0622905</td></tr><tr><td align="center" valign="middle" >BPSSM-PB</td><td align="center" valign="middle" >2.1398290</td><td align="center" valign="middle" >2.1136020</td><td align="center" valign="middle" >0.0262270</td></tr><tr><td align="center" valign="middle" >BPSSM-BCa</td><td align="center" valign="middle" >23.3015500</td><td align="center" valign="middle" >23.2929800</td><td align="center" valign="middle" >0.0085664</td></tr></tbody></table></table-wrap><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Interval estimation of standard interval: Student’s-t</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1240723x57.png"/></fig><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Interval estimation of bootstrap interval method</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1240723x58.png"/></fig><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Error value for real and hybrid model</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Model</th><th align="center" valign="middle" >MSE</th><th align="center" valign="middle" >RMSE</th></tr></thead><tr><td align="center" valign="middle" >PSSM</td><td align="center" valign="middle" >0.000574969</td><td align="center" valign="middle" >0.023978510</td></tr><tr><td align="center" valign="middle" >BPSSM</td><td align="center" valign="middle" >0.000394221</td><td align="center" valign="middle" >0.019855010</td></tr></tbody></table></table-wrap><p>Next, estimation of the effectiveness of the real and hybrid models is seen in the results of the MSE and RMSE and the good performance of the model is based on the theory of effectiveness estimation model, as discussed in the previous section. Therefore, the results of this error can be referred to <xref ref-type="table" rid="table2">Table 2</xref>.</p><p>Based on the results of <xref ref-type="table" rid="table2">Table 2</xref>, a model with the bootstrapping approach provides the smallest error value compared to the real model, PSSM. For example, the estimation result for MSE of hybrid model is 0.019855010 while real model gives large error estimation, i.e. 0.023978510. These results are plotted in <xref ref-type="fig" rid="fig3">Figure 3</xref> and <xref ref-type="fig" rid="fig4">Figure 4</xref> for an illustration comparing the results for both models of MSE and RMSE estimation value.</p><p>Based on such a significant error reduction in <xref ref-type="fig" rid="fig3">Figure 3</xref> and <xref ref-type="fig" rid="fig4">Figure 4</xref>, show the boot-</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Plot of error estimation, MSE</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1240723x59.png"/></fig><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Plot of error estimation, RMSE</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1240723x60.png"/></fig><p>strap approach in real base model of control charts fixing the estimation of real model and provide a more accurate estimation for the model. This small error values also indicate that the hybrid model is more effective and gives good performance compared to the real model, PSSM.</p></sec><sec id="s4"><title>4. Conclusion</title><p>In this study, a PSSM model was hybrid with bootstrap method using an alternative algorithm. Using an alternative algorithm, the hybrid process was involved the construction of a standard error of PSSM confidence interval and proposed a new hybrid model of BPSSM. The data set Malaysian population and Family survey 1994 (MPFS-1994) was used. Participation and wage equation on the estimated coefficient for interval method, the significant effect, and consistency and for PSSM, as well as BPSSM for comparison purposes was discussed. The differences shown in these two models proved that the bootstrap approach fixed the interval estimation and gives a good performance for hybrid model. Estimation of the effectiveness of the real and hybrid models is seen of the MSE and RMSE. This small error value also indicates that the hybrid model is more effective.</p></sec><sec id="s5"><title>Acknowledgements</title><p>A special gratitude for School of Informatics and Applied Mathematic (SIAM) and Research Management Centre (RMC), Universiti Malaysia Terengganu for supported this research paper.</p></sec><sec id="s6"><title>Cite this paper</title><p>Lola, M.S., Alwi, W.S.W. and Zainuddin, N.H. (2016) Sample Selection Model with Bootstrap (BPSSM) Approach: Case Study of the Malaysian Population and Family Survey. Open Journal of Statistics, 6, 741-748. http://dx.doi.org/10.4236/ojs.2016.65060</p></sec></body><back><ref-list><title>References</title><ref id="scirp.70799-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Heckman, J.J. (1974) Shadow Wages, Market Wages and Labor Supply. Econometrica, 42, 679-693. http://dx.doi.org/10.2307/1913937</mixed-citation></ref><ref id="scirp.70799-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Gronau, R. (1974) The Effect of Children on the Housewife’s Value of Time. Economics of the Family: Marriage, Children, and Human Capital, University of Chicago Press, 457-490.</mixed-citation></ref><ref id="scirp.70799-ref3"><label>3</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Heckman</surname><given-names> J.J. </given-names></name>,<etal>et al</etal>. (<year>1976</year>)<article-title>The Common Structure of Statistical Models of Truncation, Sample Selection and Limited Dependent Variables and a Simple Estimator for Such Models</article-title><source> Annals of Economic and Social Measurement</source><volume> 5</volume>,<fpage> 475</fpage>-<lpage>492</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.70799-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Zhou, C. (2010) The Extent of the Maximum Likelihood Estimator for the Extreme Value Index. Journal of Multivariate Analysis, 101, 971-983.  
http://dx.doi.org/10.1016/j.jmva.2009.09.013</mixed-citation></ref><ref id="scirp.70799-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Markus, F. (1998) Semi-Parametric Estimation of Selectivity Models. Unbublished Ph.D. Thesis, Konstanz University, Konstanz.</mixed-citation></ref><ref id="scirp.70799-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Efron, B. (1979) Bootstrap Methods: Another Look at the Jackknife. The Annals of Statistics, 7, 1-26. http://dx.doi.org/10.1214/aos/1176344552</mixed-citation></ref><ref id="scirp.70799-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Efron, B. (2003) Second Thoughts on the Bootstrap. Statistical Science, 18, 135-140.  
http://dx.doi.org/10.1214/ss/1063994968</mixed-citation></ref><ref id="scirp.70799-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Chou, S.-R., Chien, A., Changchien, C.-C. and Wu, C.-H. (2010) A Comparison of the Forecasting Volatility Performance on EWMA Family Models. International Research Journal of Finance and Economics, 54, 19-28.</mixed-citation></ref><ref id="scirp.70799-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Martins, F.M. (2001) Parametric and Semi-Parametric Estimation of Sample Selection Models: An Empirical Application to the Female Labour Force in Portugal. Journal of Applied Econometrics, 16, 23-39. http://dx.doi.org/10.1002/jae.572</mixed-citation></ref></ref-list></back></article>