28.05.2013 Views

(Fac-simile frontespizio tesi di laurea specialistica) - Scor

(Fac-simile frontespizio tesi di laurea specialistica) - Scor

(Fac-simile frontespizio tesi di laurea specialistica) - Scor

SHOW MORE
SHOW LESS

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

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

Un'unica procedura per la stima del premio <strong>di</strong> rischio<br />

Presentiamo adesso un'unica procedura che, partendo dai dati della sinistrosità e della<br />

frequenza del low deductible, dai valori prescelti per OP, per CO e per DE, sia in<br />

grado <strong>di</strong> calcolare automaticamente il valore del parametro α , la stima della<br />

sinistrosità, quella del costo me<strong>di</strong>o e quin<strong>di</strong> il premio <strong>di</strong> rischio:<br />

[>procRP:=proc(sins::list,OP,CO,DE,freqOP) local sinsl,<br />

OPp, COp, DEp, freqOPp, alphastp, frequenzap,<br />

valoreme<strong>di</strong>osinp, vmsp, premio<strong>di</strong>rischiop, RLp;<br />

[>sinsl:=sins;OPp:=OP;COp:=CO;DEp:=DE;freqOPp:=freqOP;<br />

[>alphastp:=nops(sinsl)/sum(ln(sinsl[i]/OPp),<br />

i=1..nops(sinsl));<br />

[>frequenzap:=freqOPp*((OPp/DEp)^alphastp);<br />

[>valoreme<strong>di</strong>osinp:=proc(COq,DEq,alphastq) local a, c, d,<br />

vms, errelle;<br />

[>a:=alphastq;c:=COq;d:=DEq;errelle:=(c+d)/d;<br />

[>if a=1 then vms:=d*ln(errelle);<br />

[>else vms:=d*(errelle^(1-a)-1)/(1-a);<br />

[>fi;<br />

[>end;<br />

[>vmsp:=valoreme<strong>di</strong>osinp(COp,DEp,alphastp);<br />

[>premio<strong>di</strong>rischiop:=(f,v)->f*v;<br />

[>RLp:=premio<strong>di</strong>rischiop(frequenzap,vmsp);<br />

[>print(`Il parametro del modello è stato stimato in`);<br />

print (alphastp);<br />

[>print(`La frequenza me<strong>di</strong>a è stimata in`);<br />

print(frequenzap);<br />

[>print(`Il costo me<strong>di</strong>o <strong>di</strong> un sinistro trasferito è<br />

stimato in`); print(vmsp);<br />

[>print(`Infine,<br />

print(RLp);<br />

il premio <strong>di</strong> rischio è stimato in`);<br />

[>end;<br />

94

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

Saved successfully!

Ooh no, something went wrong!