gemact-code icon indicating copy to clipboard operation
gemact-code copied to clipboard

Negative Multinomial aggiornata

Open Davide-Ruffini opened this issue 9 months ago • 3 comments

Ciao a tutti,

ho aggiornato la Negative Multinomial seguendo l'implementazione di wikipedia "https://en.wikipedia.org/wiki/Negative_multinomial_distribution". In passato incontrai dei problemi sull'implementazione del metodo RVS (usare il metodo solito tramite cdf non mi risultava semplicissimo essendo una multivariata, ma magari mi sono arenato io ed esiste una formula una via semplice), che mi portò a cambiare il sistema di parametrizzazione seguendo quando implementato in R nel pacchetto MGLM, ora dovrei aver trovato una soluzione che dovrebbe funzionare:

Potremmo vedere la nostra distribuzione Negative Multinomial usando: -Una binomiale negative per generare il numero totale -Una multinomiale per ripartire i conteggi

In questo modo mi sembra proprio che i momenti della distribuzione simulata si parlino con quelli teorici, e l'implementazione è semplice. Fatemi sapere cosa ne pensate

Vi allego due semplici test, nel Test1 verifico che funzioni tutto a dovere. Nel Test2 invece ho fatto qualche confronto con l'implementazione di MGLM. In effetti esiste qualche legame tra le due parametrizzazioni e si parlano quando il parametro di dispersione è 1 (formula non dispersa).

Test 2- confronto con MGLM.txt Test1.txt

Davide-Ruffini avatar May 04 '25 10:05 Davide-Ruffini

Ciao sembra che tutto sia okay adesso! Cosa ne pensi @EdoLu ?

gpitt71 avatar May 09 '25 14:05 gpitt71

Ciao! Direi che va tutto bene. @Davide-Ruffini grazie per quanto hai fatto. Ti direi come ultima cosa di modificare nelle nuove distribuzioni multivariate che hai creato il metodo statico "category" come segue

@staticmethod def category(): return {''}

in modo che non ritorni 'frequency', in quanto, al momento, non possono essere usate come componente Frequency nella classe LossModel. @gpitt71 concordi?

Edoardo

EdoLu avatar May 10 '25 15:05 EdoLu

Ciao a tutti. Ok chiarissimo, aspetto che Gabri ti risponda e procedo Davide

Davide-Ruffini avatar May 11 '25 07:05 Davide-Ruffini

@gpitt71 puoi dare il nullaosta a @Davide-Ruffini ?

EdoLu avatar May 16 '25 15:05 EdoLu

Ciao a tutti, sono d'accordo

gpitt71 avatar May 19 '25 06:05 gpitt71

Ciao a tutti! Ho applicato la modifica alle distribuzioni multivariate

Davide

Davide-Ruffini avatar May 19 '25 19:05 Davide-Ruffini

Se siamo tutti allineati possiamo fare il merge. @EdoLu

gpitt71 avatar May 20 '25 12:05 gpitt71

Scusate il ritardo, sì procederei con la PR!

Il giorno mar 20 mag 2025 alle 14:21 Gabriele Pittarello < @.***> ha scritto:

gpitt71 left a comment (gpitt71/gemact-code#12) https://github.com/gpitt71/gemact-code/pull/12#issuecomment-2894201928

Se siamo tutti allineati possiamo fare il merge. @EdoLu https://github.com/EdoLu

— Reply to this email directly, view it on GitHub https://github.com/gpitt71/gemact-code/pull/12#issuecomment-2894201928, or unsubscribe https://github.com/notifications/unsubscribe-auth/ASWH45FOANKWFBJBGCIBCKT27MM33AVCNFSM6AAAAAB4MXHV46VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDQOJUGIYDCOJSHA . You are receiving this because you were mentioned.Message ID: @.***>

EdoLu avatar May 22 '25 10:05 EdoLu

Merge completato. Grazie a tutti per il lavoro. Chiudo questa thread

gpitt71 avatar May 22 '25 11:05 gpitt71