no no l'indice di convenienza non c'entra , secondo un libro
di Leontino Gorgia di cui dispongo il termine Aspettabilità
fu introdotto da un altro famoso studioso del lotto detto "Il Fisico"
Allorquando una statistica teorica si confronta con una reale si notano
delle modeste deviazioni che al crescere del numero delle prove tendono a diminuire
mentre altre possono nascere e svilupparsi ...
il termine aspettabiltà è equivalente a quello introdotto dallo stesso
Gorgia Attendibilità.
La formula generica per il calcolo dell'aspettabilità riferita
ad una qualsiasi statistica è
Codice:
Aspettabilità = CasiReali / (CasiReali + Casi Teorici)
Ovviamente per utilizzare la formula si deve disporre di una statistica teorica.
Ad esempio se calcolassimo l'aspettabilità riferita alla statistica
del ritardo dovremmo inanzi tutto sapere quanti numeri si dovrebbero
trovare ancora in piedi a quel dato ritardo , ovvero dovremmo conoscere il valore teorico di numeri che a quel ritardo non sono ancora caduti.
Disponendo del suddetto valore che si calcola con la formula
ValTeoricoNumInPiedi = 90 * exp (-(Ritardo-1)/k)
dove k è la costante di decadimento per l'estratto semplice pari a 17,4952
quindi se troviamo un numero al ritardo 160 possiamo
calcolare quanti numeri dovrebbero esserci teoricamente al ritardo 160 con : 90 * exp(-(159)/17,4952) che fa 0 ,010
quindi sapendo che realmente ce n'è 1 possiamo calcolare l'attendibilita usando la formula gia detta
1/(1 +0,010) = 0,990
per finire un semplice script che prende in esame l'aspettabilita
sul ritardo dell'estratto , si possono avere anche numerosi altri tipi di aspettabilità ma il cocetto è sempre quello di confrontare
la statistica reale con quella teorica attraverso la formula generica
Codice:
Option Explicit
Sub Main
Dim r, n, Rit, Asp
Dim nEstrAnalisi
nEstrAnalisi =EstrazioneFin
ReDim aTitoli (4)
aTitoli(1) = "Ruota"
aTitoli(2) = "Numero"
aTitoli(3) = "Ritardo"
aTitoli(4) = "Aspettabilità"
Call InitTabella (aTitoli)
For r = 1 To 10
Call Messaggio (NomeRuota (r))
ReDim aRitardi (90)
For n = 1 To 90
Rit = RitardoEstratto(n ,nEstrAnalisi,r)
aRitardi (n) = Rit
Next
For n = 1 To 90
Rit = aRitardi(n)
Asp = CalcolaAspettabilita(Rit , QuantitaNumeriAlRit(Rit , aRitardi))
ReDim aColonne (4)
aColonne(1) = NomeRuota (r)
aColonne(2) = n
aColonne(3) = Rit
aColonne(4) = Asp
Call AddRigaTabella (aColonne)
Call AvanzamentoElab( 1 ,90 ,n)
Next
Next
Call AvanzamentoElab( 1 ,90 ,1)
Call Messaggio ("Ordinamento in corso")
Call CreaTabella (4)
End Sub
Function QuantitaNumeriAlRit(Rit , aRitardi)
Dim k , q
For k = 1 To 90
If aRitardi(k) = Rit Then
q = q +1
End If
Next
QuantitaNumeriAlRit = q
End Function
Function CalcolaAspettabilita(Rit, QNumeriAlRit)
Dim ValTeorico
ValTeorico =Round( 90 * Exp(-(Rit-1)/17.4952),4)
CalcolaAspettabilita =Round( ValTeorico / (ValTeorico + QNumeriAlRit),4)
End Function