13.07.2015 Views

Rett og riktig - En gjennomgang av Statens ... - Concept - NTNU

Rett og riktig - En gjennomgang av Statens ... - Concept - NTNU

Rett og riktig - En gjennomgang av Statens ... - Concept - NTNU

SHOW MORE
SHOW LESS
  • No tags were found...

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

63Vedlegg-A Hvordan finne Alfa,Beta <strong>og</strong> forskyvning fratrippelanslagFor å kunne benytte algoritmene som er oppgitt her må man tilgang på en funksjon for åregne ut den inverse kumulative tetthetsfunksjonen til en gammafordeling, slik som funksjonenGAMMAINV i Excel. Denne er d<strong>og</strong> ikke særlig egnet den ikke klarer verdier for alfastørre enn ca 300. Eller finnes kildekode til funksjoner som regner ut den inverse kumulativetetthetsfunksjonen til en gammafordeling fritt tilgjengelig på internett.Denne algoritmen klare å regne ut alfa for skjevhetsforhold opp til litt over 13 000 <strong>og</strong> takleri tilegg spesialtilfellet der n = s =eller s = ø. Algoritmen takler ikke tilfellet der n = s = ø dadette tilfellet ikke gir noen mening for en gammafordeling. Algoritmen forutsetter en høyreskjevkurve. Hvis man har en venstreskjev kurve må tallene transformeres først. Se“Transformasjon <strong>av</strong> trippelanslag” på side 19.Algoritme for å finne AlfaInngangsparametre: n, s <strong>og</strong> ø (nedre, mest sannsynlig <strong>og</strong> øvre anslag)1. Hvis n = s eller ø = s returner alfa = 1,15632. Regn ut skjevhetsforholdet SFnsø fra n, s <strong>og</strong> ø.3. Hvis SFnsø = 1 returner alfa = 1 000 000 000 000. Semerknad 14. Hvis SFnsø < 1, SFnsø = 1/SFnsø5. Gjør et grovt anslag på alfa. Se merknad 26. Sett alfa = 10, beta = 107. Regn ut skjevhetsforholdet SFab fra alfa <strong>og</strong> beta. Semerknad 38. Regn ut den interne forholdet mellomskjevhetsforholdene. i = SFab / SFnsø9. Hvis i = 1 ± 0,0001 <strong>av</strong>slutt <strong>og</strong> returner alfa10. Hvis vi antall iterasjoner hittil er større enn 50reduser i. Se merknad 411. Hvis i > 1 alfa = alfa * iGå tilbake til punkt 312. Hvis i < 1 alfa = alfa / (2 - i)13. Gå tilbake til punkt 3.<strong>Rett</strong> <strong>og</strong> <strong>riktig</strong> - <strong>En</strong> <strong>gjennomgang</strong> <strong>av</strong> <strong>Statens</strong> Vegvesens analysemodell

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!