Mathematical Researches

fa برآورد ناپارامتری مخاطره فضایی میدان تصادفی نامانای در میانگین Nonparametric Estimation of Spatial Risk for a Mean Nonstationary Random Field} آمار Stat مقاله استخراج شده از پایان نامه Research Paper اغلب روش‌های مرسوم برآورد مخاطره فضایی، برای میدان‌های تصادفی مانا ارائه شده‌اند و برای سادگی با پذیره‌ معلوم بودن توزیع داده‌ها یک مدل پارامتری به تابع تغییرنگار برازش داده میشود. در این مقاله یک روش فضایی ناپارامتری برای برآورد مخاطره فضایی ارائه می‌شود، که در آن  توابع روند و تغییرنگار با برآوردگرهای خطی موضعی مدل‌بندی و با تصحیح اریبی مانده‌ها  یک مدل ناپارامتری معتبر به تغییرنگار برازش داده خواهد شد. سپس با روش بوت‌استرپ، مخاطره فضایی در موقعیت‌های جدید محاسبه و پهنه‌بندی آن تهیه می‌شود.  روش فضایی ناپارامتری برای برآورد مخاطره شرطی سازوار می‌شود و با مخاطره حاصل از کریگیدن نشانگر مورد ارزیابی و مقایسه قرار می‌گیرد. به علاوه دقت روش فضایی ناپارامتری در مطالعات شبیه‌سازی و به کارگیری برای داده‌های دمای هوای ایران مورد بررسی و ارزیابی قرار می‌گیرد. Introduction Estimating the spatial hazard, or in other words, the probability of exceeding a certain boundary is one of the important issues in environmental studies that are used to control the level of pollution and prevent damage from natural disasters. Risk zoning provides useful information to decision-makers; For example, in areas where spatial hazards are high, zoning is used to design preventive policies to avoid adverse effects on the environment or harm to humans. Generally, the common spatial risk estimating methods are for stationary random fields. In addition, a parametric form is usually considered for the distribution and variogram of the random field. Whereas in practice, sometimes these assumptions are not realistic. For an example of these methods, we can point to the Indicator kriging, Disjunctive kriging, Geostatistical Markov Chain, and simple kriging.  In practice utilize the parametric spatial models caused unreliable results. In this paper, we use a nonparametric spatial model to estimate the unconditional probability or spatial risk: <m:omathpara><m:omath><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r>r</m:r></m:e><m:sub><m:r>c</m:r></m:sub></m:ssub><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>0</m:r></m:sub></m:ssub></m:e></m:d><m:r>=</m:r><m:r>P</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r>Z</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>0</m:r></m:sub></m:ssub></m:e></m:d><m:r>⩾</m:r><m:r>c</m:r></m:e></m:d><m:r>.                 (1)</m:r></m:omath></m:omathpara> Because the conditional distribution at points close to the observations has less variability than the unconditional probability, nonparametric spatial methods will be used to estimate the unconditional probability.  Material and methods Let <m:omath><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>Z</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>=</m:r><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r>Z</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>1</m:r></m:sub></m:ssub></m:e></m:d><m:r>,…,</m:r><m:r>Z</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><a name="_Hlk65962569"><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></a></m:e><m:sub><m:r>n</m:r></m:sub></m:ssub></m:e></m:d></m:e></m:d><m:ctrlpr></m:ctrlpr></m:e><m:sup><m:r>T</m:r></m:sup></m:ssup></m:omath> be an observation vector from the random field <m:omath><m:r>{</m:r><m:r>Z</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e></m:d><m:r>;</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r><m:r>∈</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>D</m:r><m:r>&sube;</m:r><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:r><m:rpr><m:scr m:val="double-struck"><m:sty m:val="i"></m:sty></m:scr></m:rpr>R</m:r></m:e><m:sup><m:r>d</m:r></m:sup></m:ssup><m:r>}</m:r></m:omath> which is decomposed as follows <m:omath><m:r>Z</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e></m:d><m:r>=</m:r><m:r>μ</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e></m:d><m:r>+</m:r><m:r>ε</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e></m:d><m:r>,</m:r></m:omath>                       (2) where <m:omath><m:r>μ</m:r><m:r>(</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r><m:r>)</m:r></m:omath> is the trend and <m:omath><m:r>ε</m:r><m:r>(</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r><m:r>)</m:r></m:omath> is the error term, that is a second-order stationary random field with zero mean and covariogram <m:omath><m:r>C</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>h</m:r></m:e></m:d><m:r>=</m:r><m:r>Cov</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r>ε</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e></m:d><m:r>,</m:r><m:r>ε</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r><m:r>+</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>h</m:r></m:e></m:d></m:e></m:d><m:r>.</m:r></m:omath> The local linear model for the trend is given by <m:omath><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r>μ</m:r></m:e></m:acc></m:e><m:sub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>H</m:r></m:sub></m:ssub><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r>s</m:r></m:e></m:d><m:r>= </m:r><m:ssubsup><m:ssubsuppr><m:ctrlpr></m:ctrlpr></m:ssubsuppr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>e</m:r></m:e><m:sub><m:r>1</m:r></m:sub><m:sup><m:r>T</m:r></m:sup></m:ssubsup><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssubsup><m:ssubsuppr><m:ctrlpr></m:ctrlpr></m:ssubsuppr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>S</m:r></m:e><m:sub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:sub><m:sup><m:r>T</m:r></m:sup></m:ssubsup><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>W</m:r></m:e><m:sub><m:r>s</m:r></m:sub></m:ssub><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>S</m:r></m:e><m:sub><m:r>s</m:r></m:sub></m:ssub></m:e></m:d></m:e><m:sup><m:r>-</m:r><m:r>1</m:r></m:sup></m:ssup><m:r> </m:r><m:ssubsup><m:ssubsuppr><m:ctrlpr></m:ctrlpr></m:ssubsuppr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>S</m:r></m:e><m:sub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:sub><m:sup><m:r>T</m:r></m:sup></m:ssubsup><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>W</m:r></m:e><m:sub><m:r>s</m:r></m:sub></m:ssub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>Z</m:r><m:r>&equiv;  </m:r><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:r>ϕ</m:r></m:e><m:sup><m:r>T</m:r></m:sup></m:ssup><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r>s</m:r></m:e></m:d><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>Z</m:r></m:omath>, where <m:omath><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>e</m:r></m:e><m:sub><m:r>1</m:r></m:sub></m:ssub></m:omath> is a vector with 1 in the first entry and all other entries 0, <m:omath><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>S</m:r></m:e><m:sub><m:r>s</m:r></m:sub></m:ssub></m:omath> is a matrix with ith row equal to (1, (si-s)T), <m:omath><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r>W</m:r></m:e><m:sub><m:r>s</m:r></m:sub></m:ssub><m:r> = </m:r><m:r>diag</m:r><m:r> </m:r><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r>{</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>K</m:r></m:e><m:sub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>H</m:r></m:sub></m:ssub><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>1</m:r></m:sub></m:ssub><m:r> – </m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e></m:d><m:r>,…,</m:r><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>K</m:r></m:e><m:sub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>H</m:r></m:sub></m:ssub><m:r>(</m:r><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>n</m:r></m:sub></m:ssub><m:r>-</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r><m:r>)}</m:r></m:omath>, <m:omath><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>K</m:r><m:ctrlpr></m:ctrlpr></m:e><m:sub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>H</m:r></m:sub></m:ssub><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>u</m:r></m:e></m:d><m:r>=</m:r><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:d><m:dpr><m:begchr m:val="|"><m:endchr m:val="|"><m:ctrlpr></m:ctrlpr></m:endchr></m:begchr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>H</m:r></m:e></m:d></m:e><m:sup><m:r>-</m:r><m:r>1</m:r></m:sup></m:ssup><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>K</m:r><m:r>(</m:r><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>H</m:r></m:e><m:sup><m:r>-</m:r><m:r>1</m:r></m:sup></m:ssup><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>u</m:r><m:r>)</m:r></m:omath>, <m:omath><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>K</m:r></m:omath> is a triple multiplicative multivariate kernel function and <m:omath><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>H</m:r></m:omath> is a nonsingular symmetric d×d bandwidth matrix. In this model, the bandwidth matrix obtained from a bias corrected and estimated generalized cross-validation (CGCV). From nonparametric residuals <m:omath><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r>ε</m:r></m:e></m:acc><m:r>(</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r><m:r>) = </m:r><m:r>Z</m:r><m:r>(</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r><m:r>) -</m:r><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>μ</m:r></m:e></m:acc><m:r>(</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r><m:r>)</m:r></m:omath> a local linear estimate of the variogram 2<m:omath><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr> </m:r><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>γ</m:r></m:e></m:acc><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>(⋅)</m:r></m:omath>is obtained as the solution of the following least-squares problem <table align="center" class="MsoTableGrid" dir="rtl" style="width:534px; border-collapse:collapse; border:none" width="534"> <tbody> <tr> <td style="width:30px; padding:0cm 7px 0cm 7px; height:40px" valign="top"></td> <td style="width:504px; padding:0cm 7px 0cm 7px; height:40px" valign="top"><m:omath><m:func><m:funcpr><m:ctrlpr></m:ctrlpr></m:funcpr><m:fname><m:limlow><m:limlowpr><m:ctrlpr></m:ctrlpr></m:limlowpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>min</m:r><m:ctrlpr></m:ctrlpr></m:e><m:lim><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>α</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>.</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>β</m:r><m:ctrlpr></m:ctrlpr></m:lim></m:limlow></m:fname><m:e><m:nary><m:narypr><m:chr m:val="∑"><m:limloc m:val="undOvr"><m:grow m:val="on"><m:ctrlpr></m:ctrlpr></m:grow></m:limloc></m:chr></m:narypr><m:sub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>i<j< m:r=""></j<></m:r></m:sub><m:sup><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>n</m:r></m:sup><m:e><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r>ε</m:r></m:e></m:acc></m:e><m:sub><m:r>i</m:r></m:sub></m:ssub><m:r>-</m:r><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r>ε</m:r></m:e></m:acc></m:e><m:sub><m:r>j</m:r></m:sub></m:ssub></m:e></m:d></m:e><m:sup><m:r>2</m:r></m:sup></m:ssup><m:r>-</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>α</m:r><m:r>-</m:r><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>β</m:r></m:e><m:sup><m:r>T</m:r></m:sup></m:ssup><m:r> </m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>i</m:r></m:sub></m:ssub><m:r>-</m:r><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>j</m:r></m:sub></m:ssub><m:r>-</m:r><m:r>u</m:r></m:e></m:d></m:e></m:d></m:e><m:sup><m:r>2</m:r></m:sup></m:ssup><m:r> </m:r><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r>K</m:r></m:e><m:sub><m:r>G</m:r></m:sub></m:ssub><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>i</m:r></m:sub></m:ssub><m:r>-</m:r><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>j</m:r></m:sub></m:ssub><m:r>-</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>u</m:r></m:e></m:d></m:e></m:nary></m:e></m:func></m:omath>,</td> </tr> </tbody> </table> where <m:omath><m:r>G</m:r></m:omath> is the corresponding bandwidth matrix, that obtained from minimizing cross-validation relative squared error of semi-variogram estimate. Algorithm1: Semiparametric Bootstrap <ol> <li style="text-align:justify">Obtain estimates of the error covariance and nonparametric residuals covariance.</li> <li style="text-align:justify">Generate bootstrap samples with the estimated spatial trend <m:omath><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r>μ</m:r></m:e></m:acc></m:e><m:sub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>H</m:r></m:sub></m:ssub><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r>s</m:r></m:e></m:d></m:omath> and adding bootstrap errors generated as a spatially correlated set of errors.</li> <li style="text-align:justify">Compute the kriging prediction <m:omath><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>Z</m:r></m:e></m:acc><m:ctrlpr></m:ctrlpr></m:e><m:sup><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>*</m:r></m:sup></m:ssup><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>0</m:r></m:sub></m:ssub></m:e></m:d></m:omath> at each unsampled location <m:omath><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>0</m:r></m:sub></m:ssub></m:omath> from the bootstrap sample <m:omath><m:d><m:dpr><m:begchr m:val="{"><m:endchr m:val="}"><m:ctrlpr></m:ctrlpr></m:endchr></m:begchr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr> </m:r><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>Z</m:r></m:e></m:acc><m:ctrlpr></m:ctrlpr></m:e><m:sup><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>*</m:r></m:sup></m:ssup><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>1</m:r></m:sub></m:ssub></m:e></m:d><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>,…,</m:r><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>Z</m:r></m:e></m:acc><m:ctrlpr></m:ctrlpr></m:e><m:sup><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>*</m:r></m:sup></m:ssup><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>n</m:r></m:sub></m:ssub></m:e></m:d></m:e></m:d></m:omath>.</li> <li style="text-align:justify">Repeat steps 2 and 3 a larger number times B. Therefore, for each un-sampled location <m:omath><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>0</m:r></m:sub></m:ssub></m:omath>, B bootstrap replications <m:omath><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>Z</m:r></m:e></m:acc><m:ctrlpr></m:ctrlpr></m:e><m:sup><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>*</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>(1)</m:r></m:sup></m:ssup><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>0</m:r></m:sub></m:ssub></m:e></m:d><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>.…. </m:r><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>Z</m:r></m:e></m:acc><m:ctrlpr></m:ctrlpr></m:e><m:sup><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>*</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>(B)</m:r></m:sup></m:ssup><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>(</m:r><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>0</m:r></m:sub></m:ssub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>)</m:r></m:omath> are obtained.</li> <li style="text-align:justify; margin-bottom:11px">Calculate (1) at position s_0 by calculating the relative frequency of Bootstrap repetition as follows to estimate the unconditional probability of excess of boundary c.</li> </ol> <table align="right" class="MsoTableGrid" dir="rtl" style="width:578px; margin-left:54px; border-collapse:collapse; border:none" width="578"> <tbody> <tr> <td style="width:116px; padding:0cm 7px 0cm 7px; height:39px" valign="top"></td> <td style="width:285px; padding:0cm 7px 0cm 7px; height:39px" valign="top"><m:omathpara><m:omath><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>r</m:r></m:e></m:acc></m:e><m:sub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>c</m:r></m:sub></m:ssub><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>0</m:r></m:sub></m:ssub></m:e></m:d><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>= </m:r><m:f><m:fpr><m:ctrlpr></m:ctrlpr></m:fpr><m:num><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>1</m:r></m:num><m:den><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>B</m:r></m:den></m:f><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr> </m:r><m:nary><m:narypr><m:chr m:val="∑"><m:limloc m:val="undOvr"><m:ctrlpr></m:ctrlpr></m:limloc></m:chr></m:narypr><m:sub><m:r>j</m:r><m:r>=1</m:r></m:sub><m:sup><m:r>B</m:r></m:sup><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>I</m:r></m:e><m:sub><m:d><m:dpr><m:begchr m:val="{"><m:endchr m:val="}"><m:ctrlpr></m:ctrlpr></m:endchr></m:begchr></m:dpr><m:e><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:acc><m:accpr><m:ctrlpr></m:ctrlpr></m:accpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>Z</m:r></m:e></m:acc><m:ctrlpr></m:ctrlpr></m:e><m:sup><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>*</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>j</m:r></m:e></m:d></m:sup></m:ssup><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e><m:sub><m:r>0</m:r></m:sub></m:ssub></m:e></m:d><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>≥ c</m:r></m:e></m:d></m:sub></m:ssub></m:e></m:nary></m:omath></m:omathpara></td> <td style="width:177px; padding:0cm 7px 0cm 7px; height:39px" valign="top"></td> </tr> </tbody> </table> Results and discussion To analysis the practical behavior of the proposed methods a simulation study is conducted under different scenarios. For N=150 samples and n=<m:omath><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>16×16</m:r></m:omath> were generated on a regular grid in the unit square following model (2), with mean function <m:omath><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>μ</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="bi"></m:sty></m:scr></m:rpr>s</m:r></m:e></m:d><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>=2.5</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr> </m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>+ sin( 2π </m:r><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>x</m:r></m:e><m:sub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>1</m:r></m:sub></m:ssub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>) + 4</m:r><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:ssub><m:ssubpr><m:ctrlpr></m:ctrlpr></m:ssubpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>x</m:r></m:e><m:sub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>2</m:r></m:sub></m:ssub><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr> - 0.5</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr> </m:r></m:e></m:d></m:e><m:sup><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>2</m:r></m:sup></m:ssup></m:omath>, and random errors normally distributed with zero mean and isotropic exponential covariogram  <m:omathpara><m:omath><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>C</m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="b"></m:sty></m:scr></m:rpr>h</m:r></m:e></m:d><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>= </m:r><m:r>0.04</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr> + 2.01 </m:r><m:d><m:dpr><m:ctrlpr></m:ctrlpr></m:dpr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>1- </m:r><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>exp</m:r></m:e><m:sup><m:f><m:fpr><m:ctrlpr></m:ctrlpr></m:fpr><m:num><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>-</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>3 ∥ </m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="b"></m:sty></m:scr></m:rpr>h</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>∥</m:r></m:num><m:den><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>0.5</m:r></m:den></m:f></m:sup></m:ssup></m:e></m:d><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>,   </m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="b"></m:sty></m:scr></m:rpr>h</m:r><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>∈ </m:r><m:ssup><m:ssuppr><m:ctrlpr></m:ctrlpr></m:ssuppr><m:e><m:r><m:rpr><m:scr m:val="double-struck"><m:sty m:val="p"></m:sty></m:scr></m:rpr>R</m:r></m:e><m:sup><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>2</m:r></m:sup></m:ssup><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>.</m:r></m:omath></m:omathpara> For comparing nonparametric spatial methods for estimate unconditional risk, conditional risk, and Indicator kriging, we considered 7 missing observations in certain situations. Empirical spatial risk and its estimates are presented in Table 1. The Indicator kriging is overestimating and estimate spatial risk larger than 1. Generally, an estimated risk with unconditional and conditional methods is near value to empirical value. Table 1. Empirical spatial risk and its estimates <table align="right" class="MsoTableGrid" dir="rtl" style="width:100.0%; border-collapse:collapse; border:none" width="100%"> <tbody> <tr> <td colspan="7" style="border-bottom:1px solid black; width:83%; padding:0cm 7px 0cm 7px; border-top:1px solid black; border-right:none; border-left:none" valign="top">Locations</td> <td style="border-bottom:none; width:16%; padding:0cm 7px 0cm 7px; border-top:1px solid black; border-right:none; border-left:none" valign="top"></td> </tr> <tr> <td style="border-bottom:1px solid black; width:10%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">(0.13, 0)</td> <td style="border-bottom:1px solid black; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">(0.87, 0.87)</td> <td style="border-bottom:1px solid black; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">(0.80, 0.20)</td> <td style="border-bottom:1px solid black; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">(0.94, 0.27)</td> <td style="border-bottom:1px solid black; width:10%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">(0, 0.47)</td> <td style="border-bottom:1px solid black; width:11%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">(074, .60)</td> <td style="border-bottom:1px solid black; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">(0.34, 0.60)</td> <td style="border-bottom:1px solid black; width:16%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">Methods</td> </tr> <tr> <td style="border-bottom:none; width:10%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.999</td> <td style="border-bottom:none; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.300</td> <td style="border-bottom:none; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.069</td> <td style="border-bottom:none; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.317</td> <td style="border-bottom:none; width:10%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.504</td> <td style="border-bottom:none; width:11%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.011</td> <td style="border-bottom:none; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.989</td> <td style="border-bottom:none; width:16%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">Empirical</td> </tr> <tr> <td style="border-bottom:none; width:10%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.998</td> <td style="border-bottom:none; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.351</td> <td style="border-bottom:none; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.054</td> <td style="border-bottom:none; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.347</td> <td style="border-bottom:none; width:10%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.494</td> <td style="border-bottom:none; width:11%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.057</td> <td style="border-bottom:none; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.954</td> <td style="border-bottom:none; width:16%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">Conditional</td> </tr> <tr> <td style="border-bottom:none; width:10%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">1.002</td> <td style="border-bottom:none; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.230</td> <td style="border-bottom:none; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.091</td> <td style="border-bottom:none; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.091</td> <td style="border-bottom:none; width:10%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.652</td> <td style="border-bottom:none; width:11%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.006</td> <td style="border-bottom:none; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.996</td> <td style="border-bottom:none; width:16%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">Indicator</td> </tr> <tr> <td style="border-bottom:1px solid black; width:10%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">1.000</td> <td style="border-bottom:1px solid black; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.388</td> <td style="border-bottom:1px solid black; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.418</td> <td style="border-bottom:1px solid black; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.481</td> <td style="border-bottom:1px solid black; width:10%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.602</td> <td style="border-bottom:1px solid black; width:11%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top"><m:omathpara><m:omath><m:r><m:rpr><m:scr m:val="roman"><m:sty m:val="p"></m:sty></m:scr></m:rpr>0.024</m:r></m:omath></m:omathpara></td> <td style="border-bottom:1px solid black; width:12%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">0.994</td> <td style="border-bottom:1px solid black; width:16%; padding:0cm 7px 0cm 7px; border-top:none; border-right:none; border-left:none" valign="top">Unconditional</td> </tr> </tbody> </table> The spatial risk mapping for the maximum temperature means of Iran in 364 stations in March 2018 is obtained. By applying Algorithm 1 final trend and semi-variogram estimates are smoother than the pilot version.  The conditional and unconditional spatial risk with 150 bootstrap replicates for two values of threshold 25 and 31 on a 75×75 grid are estimated. The unconditional risk estimate is smoother than the conditional risk estimate. Because of this in the unconditional version, biased residual unused directly in the spatial prediction but in the conditional risk estimating, original residuals and simple kriging used. Conclusion The spatial risk estimated with the nonparametric spatial method. For the trend and variability of the random field, modeling applied a local linear nonparametric model. In the simulation study, this method better results than Indicator kriging. Because the flexibility of the nonparametric spatial method could apply for the construction of confidence or prediction intervals and hypothesis testing. پهنای باند, برآوردگر خطی موضعی, بوت‌استرپ, تصحیح اریبی تغییرنگار, مخاطره فضایی. Bandwidth, Bias corrected of variogram, Bootstrap, Local linear estimator, Spatial risk. 242 265 http://mmr.khu.ac.ir/browse.php?a_code=A-10-41-22&slc_lang=fa&sid=1 Mohammad Moghadam محمد مقدم moghaddam1214@gmail.com 10031947532846005483 10031947532846005483 No Tarbiat Modares University دانشگاه تربیت مدرس Mohsen Mohammadzadeh محسن محمدزاده mohsen_m@modares.ac.ir 10031947532846005484 10031947532846005484 Yes Tarbiat Modares University دانشگاه تربیت مدرس