Tag Archives: Hermun

Eitt, tvö eða þrjú vasaljós?

Það komst í fréttirnar í vikunni þegar þremur konum var bjargað úr helli. Í viðtali við Reykjavík síðdegis á Bylgjunni af því tilefni sagði Árni B. Stefánsson, augnlæknir og landskunnur hellakönnuður:

Það er regla hjá hellamönnum að fara aldrei í helli nema með þrjú ljós með sér. ... Af því að aðalljósið gæti klikkað og varaljósið gæti klikkað. En það er mjög ólíklegt að þriðja ljósið fari líka. Það er bara líkindareikningur.

Í framhaldinu ræddu útvarspmennirnir svo um búnað, aðstæður í hellaferðum og hverju fljótt göngumenn þreytast í hellum svo eitthvað sé nefnt. Spurningunni um það til hve mikilla bóta það er að bera aukaljós með sér er enn ósvarað.

Uppbygging og áreiðanleiki kerfa

Raðtengt kerfi með tveimur kerfishlutum.

Raðtengt kerfi með tveimur kerfishlutum.

Í greiningu á áreiðanleika þarf að skoða hvern hluta kerfis og uppbyggingu þess. Með raðtengingu kerfishluta er átt við að kerfið virki ef og aðeins ef allir hlutar þess virka. Kerfið á myndinni til hægri er sett saman úr tveimur kerfishlutum, A og B. Það mun aðeins virka á meðan báðir hlutar þess eru í lagi. Þegar annað hvort A eða B bilar þá stoppar kerfið.

Ef við lýsum líftíma hvors hvors kerfishluta um sig með slembibreytunum  og þá er líftími alls kerfisins jafn þeim hluta sem endist skemur, .

Samsíða kerfi með tveimur kerfishlutum.

Samsíða kerfi með tveimur kerfishlutum.

Það er talað um hliðtengt kerfi ef það virkar á meðan einhver hluti þess virkar. Kerfið á myndinni til vinstri er einnig sett samsett úr tveimur hlutum, A og B. Þessari framsetningu er ætlað að tákna að kerfið mun virka á meðan annar hvor hluti þess virkar. Þegar báðir hafa bilað þá hættir kerfið að virka.

Eins og í fyrra dæminu er líftíma hlutanna lýst með slembibreytunum  og . Í hliðtengdu kerfi er líftími kerfisins jafn líftíma þess hluta sem dugir lengst, .

Kerfi með varakerfi.

Kerfi með varakerfi.

Þriðja uppsetningin lýsir kerfi, sem notar A sem aðalkerfi og B er til vara. Í einfaldaðri uppsetningu gerum við ráð fyrir að varakerfið bili ekki á meðan aðalkerfið er í gangi. Táknmyndunum er ætlað að lýsa því að þegar aðalkerfið bilar mun varakerfið taka yfir. Þegar varakerfið bilar þá stoppar kerfið.

Líftími alls kerfisins er jafn summu líftíma hvors hluta um sig, .

Hve mörg vasaljós?

Þrjú hliðtengd vasaljós.

Þrjú hliðtengd vasaljós.

Gerum ráð fyrir að vasaljós hafi eina rafhlöðu og eina ljósaperu, sem við táknum með bókstöfunum R og L í raðtengdu kerfi. Bæði rafhlaðan og ljósaperan þurfa að vera í lagi svo að vasaljósið lýsi. Gerum svo ráð fyrir að hellakönnuðirnir noti þrjú vasaljós í tveimur ólíkum uppsetningum.

Ef þrír hellakönnuðir nota hver sitt vasaljós getum við táknað kerfið eins og sýnt er á bláu myndinni til vinstri. Við gerum ráð fyrir því að ef ljós eins hellakönnuðar bilar þá sé hann studdur áfram af öðrum í hópnum. Svo lengi sem a.m.k. eitt ljós er í lagi geta hellakönnuðirnir haldið för sinni áfram. Köllum þetta kerfi 1.

Eitt aðalljós og tvö varaljós.

Eitt aðalljós og tvö varaljós.

Ef hellakönnuðir nota eitt aðalljós en bera tvö önnur með sér til vara þá getum við táknað kerfið eins og sýnt er á grænu myndinni hérna til hægri. Þegar aðalljósið bilar er gripið til fyrra varaljóssins. Þegar þegar það bilar er seinna varaljósið notað. Þetta er kerfi 2.

Líftími metinn með hermun

Til þess að meta ábatann af því að fjölga vasaljósum notum hermum við mögulegar útkomur og skoðum dreifingarnar. Gerum ráð fyrir að líftími rafhlaða lúti normaldreifingu með meðaltal 5 klukkustundir og staðalfrávik hálfa klukkustund. Látum líftíma ljósapera lúta veldisdreifingu með meðaltal 10 klukkustundir.

Uppsafnaðar líkur líftíma fyrir þrjú kerfi.

Uppsafnaðar líkur líftíma fyrir þrjú kerfi.

Niðurstöður hermunar sýna að eitt vasaljós dugir að meðaltali í 3,9 klst. og þrjú hliðtengd vasaljós í kerfi 1 að meðaltali í 5,1 klst. Líftími kerfis 2 með eitt aðalljós og tvö til vara dugar er að meðaltali í 11,8 klst. Það er u.þ.b. þrefaldur líftími þegar aðeins eitt ljós er notað.

Sér í lagi höfum við áhuga á lágum útkomum af því að það kemur göngumönnum í koll þegar ljósin bila. Í 5% prósent tilfella er líftími eins vasaljóss styttri en hálf klukkustund. Líftími kerfis nr. 1 er í 5% tilfella styttri en 4,1 klst. og 6,5 klst. með kerfi 2.

Sýni af normaldreifðu slembibreytunni

Í keppni um athygli bjóða fyrirtæki og aðrir aðilar nýjum vinum að tengjast sér í gegnum samfélagsmiðla og heita gjarnan verðlaunum til að örva aðsóknina. Af gamansemi var lofað að draga út milljón sýni af normaldreifðu slembibreytunni þegar Actuary.is næði því marki að eiga yfir eitthundrað vini á Facebook-síðunni. Það var allgóð aðsókn að vefsíðunni þegar frétta- og bloggvefir bentu á umfjöllun um nettó eigna- og skuldastöðu kynslóðanna fyrr í mánuðinum. Margir gerðust vinir síðunnar með því að smella á Like-hnappinn hér uppi til hægri. Vinir síðunnar eru núna orðnir gott betur en eitthundrað. Þeim, sem hafa komið efninu á framfæri við vini og venslafólk, er þökkuð aðstoðin við útbreiðsluna.

Öllu gamni fylgir nokkur alvara og nú er komið að efndum. Af praktískum ástæðum verða sýnin færri en milljón en í stað þess verður spyrt við umjöllunina hvernig þau eru framkölluð.

Eitt sinn var sögð saga frá árdögum tölvuútreikninga við spágerð í veðurfræði, sem kann að hafa byggt á e-s konar hermun. Það fylgdi sögunni að útreikningarnir hefðu krafist mikils reikniafls og veðrið hefði í reynd verið komið löngu á undan spánni. Með tilkomu aukins reikniafls í einkatölvum og tölvuklösum opnast ótal möguleikar. Margir þekkja möguleika með innbyggðum föllum og viðbótarpökkum í Excel töflureikninum. Þá bjóða forritunarmál eins og R-hugbúnaðurinn upp á öfluga hermunartól, sem eru þjál í notkun.

Hermun er notuð ýmsum viðfangsefnum. Í grunninn gengur hermun á því að búa til reikniverk þar sem inntakið eru óþekktar stærðir, sem við þekkjum betur sem slembibreytur. Óþekktu stærðirnar geta lýst hlutfallsbreytingu á virði verðbréfa á einum degi, fjöldi jarðsjálfta á ákveðnu tímabili eða matsfjárhæð skaða í líkams- eða eignatjóni. Einstakar útkomur slembibreytanna eru óþekktar en upplýsingar um mögulegar útkomur eru þekktar eða áætlaðar. Í þessum dæmum getur útkoman, sem við höfum áhuga á, lýst virði kaupréttar á verðbréf eða væntigildi bóta til þess sem kaupir sér tryggingu fyrir fjárskaða vegna tjóns.

Líkindafall og líkindadreififall fyrir útkomur sex hliða tengings

Líkindafall og líkindadreififall fyrir útkomur sex hliða tengings.

Slembibreytum er lýst með líkinda-, líkindadreifi- og líkindaþéttiföllum. Einfalt dæmi er útkoman þegar sex hliða tengini er kastað. Líkindafallið táknar líkur á hverri og einni útkomu en líkindaþéttifallið táknar líkur þess að útkoman sé minni en eða jafnt og tiltekið gildi. Líkindafall og líkindaþéttifall fyrir útkomur fullkomins sex hliða tenings er sýnt á myndinni til hægri. Útkoma teningsins er dæmi um stakræna slembibreytu. Út úr líkindafallinu á efri myndinni getum við lesið að líkur á hverri og einni útkomu eru 1/6, gjarnan táknað sem 0,1667 eða 16,67%. Af neðri myndinni má lesa líkindi þess að útkoman sé lægri en eða jöfn tölunni á lárétta ásnum.  Þannig eru líkurnar 0,50 eða 50% að útkoman sé lægri en eða jöfn tölunni 3 (þ.e.a.s. að útkoman sé 1, 2 eða 3).

Líkindadreififall og líkindaþéttifall normaldreifingar

Líkindadreififall og líkindaþéttifall normaldreifingar.

Samfelldum slembibreytum er lýst á áþekkan hátt og myndirnar til vinstri sýna líkindadreifi- og líkindaþéttifall normaldreifðu slembibreytunnar með meðaltal núll og einingastaðalfrávik. Til þess að finna líkur þess að útkoma sé að ákveðnu bili má reikna flatarmálið undir ferlinum innan þess bils. Líkur á einni einstakri útkomu eru núll. Að sama skapi sést að líkindadreififall normaldreifðu slembibreytunnar er ekki þreplaga líkt og dreififallið fyrir fullkomna teninginn. Við getum lesið af lóðrétta ás líkindadreififallsins líkur þess útkoma slembibreytunnar sé lægri eða jöfn gildinu á lárétta ásnum.

Notkun líkindadreififallsins er grunnurinn að þeirri hermunaraðferð, sem lagt er upp með hér og nefnist aðferð andhverfrar vörpunar (inverse method). Hana má nota þegar líkindadreififallið er andhverfanlegt. Aðferðin gengur út á að safna fyrst sýnum jafdreifðu slembibreytunnar og varpa þeim síðan yfir dreifinguna, sem sóst er eftir.

Hreyfimyndin hér fyrir neðan sýnir hermun tíu þúsund sýna af normaldreifðu slembibreytunni á myndrænan hátt. Við getum ímyndað okkur að fyrst rigni niður tíu þúsund sýnum, sem eru jafndreifð á bilinu 0 til 1 og við táknum það með bláu punktunum. Næst vörpum við þeim sýnum niður á talnalínu normaldreifingarinnar með því að finna skurðpunkt líkindadreififallsins við láréttra línu frá sérhverju sýni. Frá skurðpunktinum drögum við lóðrétta línu niður á talnalínu normaldreifingarinnar. Eftir vörpunina má sjá að dreifing sýnanna fellur vel að líkindaþéttifalli normaldreifingarinnar.

English summary:

Above you will find brief discussion about probability functions, probability density functions and probability distribution functions for random variables. As an example both probability function and probability distribution function are drawn for perfect dice to represent an arbitrary discrete random variable. As well, you will recognize the probability density function and probability distribution function for the normal random variable as example of frequently used discrete random variable.

The purpose of this post is to simulate ten thousand samples from the normal distribution using the inverse method. In the beginning, ten thousand samples from the uniform distribution are generated using random number generator within the R software package. Samples from the uniform distribution are represented with blue dots in the animated picture above. Each sample is transformed using the inverse of the normal probability distribution function. As the red samples stack up on top of previous samples we can recognize the shape of our imaginary pile of random numbers approximates the bell shaped probability density function for the normal random variable.