فهرست مطالب

نشریه اندیشه آماری
سال بیست و یکم شماره 2 (پیاپی 42، پاییز و زمستان 1395)

  • تاریخ انتشار: 1395/12/12
  • تعداد عناوین: 10
|
  • عباس پرچمی صفحات 1-10
    در این مقاله، مرور و مقایسه ای بر دو بسته نرم افزاری FuzzyNumbers و Calculator.LR.FNs انجام شده است. این بسته ها قابلیت نصب بر روی نرم افزار R را دارند و در واقع ابزار و توابعی در اختیار کابران شان قرار می دهند که به وسیله آنها رسم و انجام عملیات حسابی بر روی اعداد فازی LR به سادگی میسر می شود. برای درک بهتر خوانندگان، در این مقاله سعی شده است تا دستور العمل ها و مقایسات به وسیله مثال های عددی زیادی مطرح و توضیح داده شوند و پیشنهاد می شود که خوانندگان دستورات مثال ها را در نرم افزار R عملا اجرا و دنبال کنند.
    کلیدواژگان: بسته نرم افزاری R، عدد فازی LR، عملگر حسابی، اصل گسترش
  • علی آقامحمدی، سکینه محمدی صفحات 11-23
    مدل های داده های پانلی پویا قسمت مهمی از مطالعات حوزه های پزشکی، اجتماعی و اقتصادی را شامل می شوند.
    ویژگی بارز این مدل ها وجود متغیر وابسته تاخیری به عنوان متغیر توصیفی است. مشکل برآورد در این مدل ها از همبستگی بین متغیر وابسته تاخیری و مولفه خطای فعلی ناشی می شود. اخیرا رگرسیون چندکی تاوانیده برای تحلیل داده های پانلی پویا مورد توجه قرار گرفته است. در این مقاله نخست مدل رگرسیون چندکی با ایجاد تاوان لاسو سازوار روی اثر های تصادفی برای داده های پانلی پویا با فرض وابستگی اثر های تصادفی و مشاهدات اولیه ارائه می شود. همچنین این مدل با فرض استقلال بین اثر های تصادفی و مشاهدات اولیه نیز بررسی خواهد شد. هر دو مدل از دیدگاه آمار بیزی بیان شده، مورد تحلیل قرار می گیرند. چون در این دو روش، توزیع پسین پارامترها به شکل بسته قابل حصول نیست، توزیع های پسین شرطی کامل پارامترها محاسبه و از الگوریتم نمونه گیری گیبز برای استنباط استفاده می شود. برای مقایسه کارایی روش های بیزی ارائه شده با روش های متداول، مطالعه شبیه سازی انجام شده و در پایان نیز روش استفاده از مدل ها در قالب مثال کاربردی شرح داده خواهد شد.
    کلیدواژگان: استنباط بیزی، تاوان لاسو سازوار، توزیع لاپلاس نامتقارن، داده های پانلی پویا، رگرسیون چندکی، نمونه گیری گیبز
  • مسعود قاسمی بهجانی، میلاد اسدزاده صفحات 25-31
    در این مقاله شیوه تعیین اندازه نمونه مطلوب بر اساس تابع زیان نامتقارن لینکس کراندار به روش بیزی برای توزیع های نرمال، نمایی و پوآسن بیان شده است. اندازه نمونه مطلوب با استفاده از روش عددی محاسبه شده است. در روش عددی، نخست میانگین ریسک پسین را به دست آورده و آن را به تابع هزینه خطی لیندلی اضافه می کنیم تا میانگین هزینه کل به دست آید. سپس نمودار اندازه نمونه را در مقابل میانگین هزینه کل رسم می کنیم و در نهایت، اندازه نمونه مطلوب را که هزینه را مینیمم می کند، به دست می آوریم.
    کلیدواژگان: تابع هزینه خطی لیندلی، توزیع پسین، توزیع نرمال و ریسک پسین
  • مجید جانفدا، داود شاهسونی صفحات 33-44
    امکان مطالعه بسیاری از پدیده های علمی و طبیعی در شرایط آزمایشگاهی میسر نیست و لذا آنها را در قالب مدل های ریاضی بیان و با کدهای کامپیوتری پیچیده، شبیه سازی می کنند. اجرای مدل کامپیوتری با ورودی های متفاوت را آزمایش کامپیوتری می نامند. مباحث آماری، در آزمایش های کامپیوتری از جایگاه ویژه ای برخوردار است.
    در این مقاله ضمن تشریح ساختار این مدل ها، به معرفی وبیان اهمیت تحلیل حساسیت واریانس مبنا می پردازیم. تحلیل حساسیت، مجموعه روش هایی است که اثربخش بودن پارامترهای ورودی را بر خروجی مدل با شاخص های حساسیت تعیین می کنند. این شاخص ها بر اساس مفاهیم واریانس های شرطی بیان می شوند. از آن جا که شکل ریاضی این مدل ها، به صورت صریح مشخص نیست؛ مسئله برآورد این شاخص ها با روش های مونته کارلو مبنا مطرح می گردد. از طرف دیگر زمان اجرا، چالش جدی در مدل های کامپیوتری محسوب می شود. طرح نقاط آزمایش ویژه ای مبتنی بر اعداد شبه تصادفی، برای کاهش زمان اجرای مدل پیشنهاد شده است. به منظور پرداختن به جنبه عملی، از مدل INCA-N که میزان آلایندگی نیتروژن ورودی به آب رودخانه ها را شبیه سازی می کند، استفاده شده است تا بتوان با شاخص های حساسیت متغیرهای تاثیرگذار بر این عامل تهدید کننده سلامت انسان و محیط زیست را شناسایی کرد.
    کلیدواژگان: مدل تعیینی، تحلیل حساسیت واریانس مبنا، اعداد شبه تصادفی، مدل INCA-N
  • مرضیه باغبان سیچانی صفحات 45-52
    این مقاله مطالعه ای در اندازه اهمیت توام،JRI ، اجزای در سیستم های منسجم k از n با اجزای مستقل از هم است. JRI نرخ بهبود قابلیت اعتماد سیستم به سبب بهبود قابلیت اعتماد دو یا چند جزء است. علامت و مقدار JRI نشان دهنده نوع و میزان تاثیر متقابل اجزای بر حسب قابلیت اعتماد سیستم است. البته در حالاتی خاص می توان علامت JRI را بدون محاسبه تعیین کرد. از این اندازه در تعیین اولویت اجزای برای تعمیر و نگهداری پیشگیرانه سیستم ها و تعیین سطح افزونگی سیستم های k از n استفاده می شود.
    کلیدواژگان: افزونگی، اندازه اهمیت توام، تعمیر، سیستم های k از n
  • معصومه فرخی، محمدرضا ربیعی، محمد آرشی صفحات 53-61
    در این مقاله روش رگرسیون ستیغی فازی موزون weighted جدیدی برای مجموعه ای از داده های ورودی دقیق و خروجی فازی مثلثی پیشنهاد شده است. برای این منظور، برآوردگر ستیغی پارامترهای فازی مدل رگرسیونی به دست آمده و خطای پیشگوی predicted error آن با استفاده از نرم فازی موزون به ازای ضرایب ستیغی دقیق، متفاوت ارزیابی شده است. برای ارزیابی مدل رگرسیونی پیشنهادی، ضریب تعیین تعمیم یافته را معرفی می کنیم. سپس با ارائه مثال های عددی به مقایسه مدل رگرسیون فازی معمولی و رگرسیون ستیغی فازی به کمک مقدار میانگین توان دوم خطاهای پیشگویی و ضریب تعیین فازی می پردازیم. در صورت وجود هم خطی بین متغیرهای پیش بین، نتایج به دست آمده از مدل رگرسیون ستیغی فازی بهتر از نتایج رگرسیون فازی معمولی خواهد بود.
    کلیدواژگان: اعداد فازی مثلثی، رگرسیون فازی، برآوردگر ستیغی، رگرسیون ستیغی فازی موزون، ضریب تعیین تعمیم یافته
  • رامین کاظمی، الهه نادری صفحات 63-69
    ترکیبیات تحلیلی تلاشی برای توانمند ساختن پیش بینی های کمی ویژگی های ساختارهای ترکیبیاتی بزرگ است. این نظریه در دهه های اخیر به عنوان پایه ای برای تحلیل الگوریتم ها و مطالعه مدل های علمی در بسیاری از رشته ها شامل نظریه احتمال، فیزیک آماری، زیست شناسی محاسباتی و نظریه اطلاع ظاهر شده است.
    با یک ترکیب دقیق روش های ارزیابی نمادین، آنالیز مختلط، توابع مولد و تحلیل نقطه زینی، این نظریه برای مطالعه ساختارهای پایه ای نظیر جایگشت ها، دنباله ها، رشته ها، قدم زدن، مسیرها، درخت ها، گراف ها و نقشه ها به کار گرفته می شود. هدف این مقاله ، معرفی گام های ترتیبی یک ترکیبیات تحلیلی است.
    کلیدواژگان: رده، تابع مولد، تبدیل ملین، تبدیل پوآسون، تحلیل نقطه زینی
  • مریم شکری ساز، حمیدرضا نواب پور صفحات 71-81
    در بسیاری از آمارگیری ها، برخی از واحدهای نمونه به تعدادی از پرسش ها یا همه آنها پاسخ نمی دهند. این امر موجب بروز بی پاسخی می شود. اریبی و تورم واریانس، دو اثر مهم بی پاسخی بر آماره های آمارگیری هستند. اگرچه افزایش اندازه نمونه از تورم واریانس براوردها جلوگیری می کند، لزوما اثری بر کاهش اریبی آماره ها ندارد. از این رو روش های مختلفی برای تعدیل اریبی بی پاسخی مورد استفاده قرار می گیرد. هنگامی که ساختار گم شدگی تصادفی باشد، تعدیل موزون برای جبران اثر بی پاسخی واحد آماری مناسب است. یکی از روش های وزن دهی، روش امتیاز تمایل است. تخصیص وزن در روش امتیاز تمایل، بر مبنای براورد احتمال پاسخ واحدهای نمونه ای انجام می شود. این براوردها با برازش مدل های پارامتری مناسب به دست می آیند. در این مقاله روش امتیاز تمایل و براوردگرهای تعدیل شده حاصل از آن معرفی می شوند. سپس به مقایسه عمل کرد سه براوردگر تعدیل شده امتیاز تمایل پرداخته می شود. در آخر با استفاده از مجموعه داده های آمارگیری هزینه و درامد خانوارهای شهری مرکز آمار ایران در سال 1390 ، براوردگرهای تعدیل شده امتیاز تمایل از نظر معیارهای ریشه دوم میانگین توان دوم خطای نسبی و کارایی نسبی با هم مقایسه می شوند.
    کلیدواژگان: رگرسیون خطی فاصله ای، فاصله های فازی، شمول کل، تشخیص مدل، برنامه ریزی خطی
  • شهرستانی شهرام یعقوب زاده صفحات 83-87
    در این مقاله برآورد E- بیزی پارامترهای توزیع نمایی دوپارامتری تحت تابع زیان درجه دوم به دست می آید. سپس با استفاده از شبیه سازی مونته کارلو برآورد E- بیزی پارامترها با برآوردهای بیزی مقایسه می شوند.
    کلیدواژگان: توزیع نمایی دوپارامتری، برآورد E، بیزی، شبیه سازی مونته کارلو
  • زینب بهباش، غلامعلی پرهام صفحات 89-100
    توابع مفصل می توانند ساختار وابستگی بین متغیرها را به صورت مدل نشان دهند. تابع مفصل مناسب برای یک کاربرد خاص، تابعی است که به بهترین نحو ممکن وابستگی بین داده ها را نشان دهد. آزمون های نیکویی برازش از دیدگاه نظری بهترین روش انتخاب تابع مفصل می باشند. روش های آزمون نیکویی برازش متعددی برای توابع مفصل پیشنهاد شده است. در این مقاله در پی انتخاب روش های آزمون نیکویی برازش مناسب، سه روش متفاوت آزمون نیکویی برازش برای تابع مفصل را مورد بررسی و مقایسه عددی قرار داده، به بررسی نقاط قوت و ضعف آن ها نسبت به یکدیگر می پردازیم. در نهایت روش های آزمون مطرح شده را با استفاده از داده های شاخص بورس اوراق بهادار تهران تحلیل قرار می کنیم.
    کلیدواژگان: تابع مفصل، آزمون نیکویی برازش
|
  • Abbas Parchami Pages 1-10
    ýThis paper has been discussed and reviewed two R packages FuzzyNumbers and Calculator.LR.FNsý. ýThese packages have the ability of installation on R softwareý, ýand in fact they propose some useful instruments and functions to the users for draw and easily using arithmetic operators on LR fuzzy numbersý. ýFor the convenience of the readersý, ýthe proposed methods and functions have been presented with several numerical examples in this paper which can help to better understandingý.
    Keywords: ýR packageý, ýLR fuzzy numberý, ýArithmetic operatorsý, ýExtension principleý
  • Ali Aghmohammadi, Sakine Mohammadi Pages 11-23
    ýDynamic panel data models include the important part of medicineý, ýsocial and economic studiesý. ýExistence of the lagged dependent variable as an explanatory variable is a sensible trait of these modelsý. ýThe estimation problem of these models arises from the correlation between the lagged depended variable and the current disturbanceý. ýRecentlyý, ýquantile regression to analyze dynamic panel data has been taken in to considerationý. ýIn this paperý, ýquantile regression model by adding an adaptive Lasso penalty term to the random effects for dynamic panel data is introduced by assuming correlation between the random effects and initial observationsý. ýAlsoý, ýthis model is illustrated by assuming that the random effects and initial values are independentý. ýThese two models are analyzed from a Bayesian point of viewý. ýSinceý, ýin these models posterior distributions of the parameters are not in explicit formý, ýthe full conditional posterior distributions of the parameters are calculated and the Gibbs sampling algorithm is used to deductioný. ýTo compare the performance of the proposed method with the conventional methodsý, ýa simulation study was conducted and at the endý, ýapplications to a real data set are illustratedý.
    Keywords: ýAdaptive Lasso penaltyý, ýAsymetrice Laplace distributioný, ýBayesian inferenceý, ýDynamic panel dataý, ýGibbs samplingý, ýQuantile regressioný
  • Masoud Ghasemi Behjani Pages 25-31
    ýIn this articleý, ýthe method of determining the optimal sample size is based on Linex asymmetric loss function and has been expressed through Bayesian method for normalý, ýPoisson and exponential distributionsý. ýThe desirable sample size has been calculated through numerical methodý. ýIn numerical methodý, ýthe average posterior risk is calculated and then it is added to the Lindley linear cost function to achieve the average of the total costý. ýThený, ýthe diagram of sample size is drawn in comparison to the average of total cost and eventuallyý, ýthe optimal sample size that minimizes the cost has been achieved.
    Keywords: ýLindely linear cost functioný, ýnormal distributioný, ýposterior distribitioný, ýposterior riský
  • Majid Janfada, Dr Davood Shahsavani Pages 33-44
    ýThe study of many scientific and natural phenomena in laboratory condition is sometimes impossibleý, ýtherefore theire expresed by mathemathical models and simulated by complex computer models (codes)ý. ýRunning a computer model with different inputs is called a computer exprimentý.
    ýStatistical issues allocated a wide range of applications for computer expriment to itselfý. ýIn this paperý, ýtheý ýstructure of computer models is describedý, ýand one of statistical applicationsý, ýthat is variance-based sensitivity analysis is expressedý. ýSensitivity analysisý, ýinvolves a set of methods that determine the effect on model inputs on the output by using sensitivity indicesý. ýThe indices are defined based on the conceptý ýof condition variance and the since explicit mathematical form of the model is unclearý, ýhence the essues monte carlo based them are proposedý.
    ýDue to the inherent complexity of the modelý, ýexecuation time is problem.Therefore a specifict design of exprimentý, ýbase on Quasi-random numberý, ýis proposedý ýto reduce the runnig costsý. ýAs an applicationý, ýthe INCA-N model that simulates amount of Nitrogen in riverý ýand underground sources was usedý. ýUsing the sensitivity indicesý, ýwe could found the effectiveý ývariable on this danger pollution that threaten human life and inviromentalý.
    Keywords: ýdeterministic modelý, ývariance based sensitivity analysisý, ýquasi-random numberý, ýINCA-N modelý
  • Pages 45-52
    ýIn this paperý, ýthe concept of joint reliability importance (JRI) of two or groups of components in a coherent system with independent components have been studiedý. ýThe JRI is defined as the rate at which the system reliability improves as the reliabilitiesý ýof the two or groups of components improveý.
    ýGenerallyý, ýthe sign and the value of the JRI represent the type and theý ýdegree of interactions between components with respect to systems reliabilityý.
    Keywords: ýredundancyý, ýjoint reliability importanceý, ýk-out-of-n systemsý
  • Pages 53-61
    ýIn this paper a new weighted fuzzy ridge regression method for a given set of crisp input and triangular fuzzy output values is proposedý. ýIn this regardý, ýridge estimator of fuzzy parameters is obtained for regression model and its prediction error is calculated by using the weighted fuzzy norm of crisp ridge coefficientsý. . ýTo evaluate the proposed regression modelý, ýwe introduce the fuzzy coefficient of determination (FCD)ý. ýFuzzy regression is compared with its ridge version by using mean predict error and FCDý, ýnumericallyý. ýIt is evident from comparison results the proposed fuzzy ridge regression is superior to the non-ridge counterpar.
    Keywords: ýTriangular fuzzy numbersý, ýfuzzy regressioný, ýridge estimatorý, ýweighted fuzzy ridge regressioný, ýcoefficient of determinationý
  • Pages 63-69
    ýAnalytic combinatorics aims to enable precise quantitative predictions of the properties of large combinatorial structuresý. ýThis theory has emerged over recent decades as essential both for the analysis of algorithms and for the study of scientific models in many disciplinesý, ýincluding probability theoryý, ýstatistical physicsý, ýcomputational biology and information theoryý. ýWith a careful combination of symbolic enumeration methodsý, ýcomplex analysisý, ýgenerating functions and saddle point analysisý, ýit can be applied to study of fundamental structures such as permutationsý, ýsequencesý, ýstringsý, ýwalksý, ýpathsý, ýtreesý, ýgraphs and mapsý. ýThis paper aims to introduce the order steps of an analytic combinatorics.
    Keywords: ýclassý, ýgenerating functioný, ýMellin transformý, ýPoisson transformý, ýsaddle point analysisý
  • Maryam Shekarisaz, Hamidreza Navvabpour Pages 71-81
    ýIn many statistical studies some units do not respond to a number or all of the questionsý. ýThis situation causes a problem called non-responseý. ýBias and variance inflation are two important consequences of non-response in surveysý. ýAlthough increasing the sample size can prevented variance inflationý, ýbut cannot necessary adjust for the non-response biasý. ýTherefore a number of methods are used for reducing non-response effectsý. ýIn the cases where missing mechanism is at randomý, ýweighting adjustment is an appropriate method for compensating the effects of unit non-responseý. ýPropensity score is a weighting method in which weight allocation is accomplished based on the estimates of response probabilitiesý. ýThese estimates are obtained by fitting suitable parametric modelsý. ýIn this paperý, ýthe propensity score method and its resulted adjusted estimators are introducedý. ýThen we compare the performance of three propensity score adjusted estimatorsý. ýFinallyý, ýdata on Household Income and Expenditure Survey for urban families conducted by Statistical Centre of Iran in spring 1390 are used to compare the adjusted propensity score estimators by two measures of comparisonsý, ýroot relative mean squared error and relative efficiencyý.
    Keywords: ýnon-responseý, ýweightingý, ýpropensity scoreý, ýroot relative mean squared errorý, ýaugmented response modelý, ýmissing at randomý
  • Shahrastani Shahram Yaghoobzadeh Pages 83-87
    ýIn this studyý, ýE-Bayesian of parameters of two parameter exponential distribution under squared error loss function is obtainedý. ýThe estimated and the efficiency of the proposed method has been compared with Bayesian estimator using Monte Carlo simulationý.
    Keywords: ýtwo parameter exponential distributioný, ýE-Bayesian estimationý, ýMonte Carlo simulationý
  • Pages 89-100
    ýCopula functions as a model can show the relationship between variablesý. ýAppropriate copula function for a specific application is a function that shows the dependency between data in a best wayý. ýGoodness of fit tests theoretically are the best way in selection of copula functioný. ýDifferent ways of goodness of fit for copula existý. ýIn this paper we will examine the goodness of fit tests from theoretical point of view and evaluate three different methods for comparing the copula functions as well as numerical comparison in order to show the advantage and weak points of each methodý. ýAt the end we will analyze the methods of discussed test by using the information from Tehran Stock Exchangeý.
    Keywords: ýcopula functionsý, ýgoodness of fit testsý