SIR μοντέλο, Chi-Square και Στατιστική Ανάλυση

ψήφοι
0

Για μια επιδημία με πραγματικά δεδομένα .

Infected 3,6,25,73,222,294,258,237,191,125,69,27,11,4
Day      1,2,3,4,5,6,7,8,9,10,11,12,13,14 

Μπορώ να χρησιμοποιήσω το μοντέλο SIR να βρείτε 3 παραμέτρους: βήτα, γάμμα, δέλτα

Χρησιμοποιώ: Gauss-Newton, και Levenberg-Marquardt (προσαρμογή καμπύλες, μη γραμμική παλινδρόμηση) Προσομοίωση στη γλώσσα προγραμματισμού Python, SciPy βιβλιοθήκη. Η βελτιστοποίηση τιμή για παραμέτρους: βήτα = 0,0022, γάμμα = 0.4477, Ν = 800.00

Ο αριθμός των μολυσμένων ανά ημέρα που παραδίδει το θεωρητικό μοντέλο με τις βέλτιστες παραμέτρους είναι

Infected [3.0, 9.99280328, 32.06524767,   92.00309752,  201.83944433,
          291.31202646,  288.30631978,  228.92526434,  163.30817847,  
          110.59508837, 72.9160331,    47.37239594,   30.51768855,  
          19.56086609]
Day      1,2,3,4,5,6,7,8,9,10,11,12,13,14 

Δύο ερώτηση:

1-) Chi-square test της καλοσύνης-of-fit, Χρήση: scipy.stats.chisquare (f_obs, f_exp = Κανένα, ddof = 0, άξονας = 0)

Where:
f_obs=real data
f_exp=theoretical model
ddof=0, axis=0
So
X^2 =52.991019453928587
and
The p-value is 1.84146232217709e-11 

Πραγματικά δεν ξέρω αν κάνω καλά και να υποστηρίξει ή να απορρίψει τη μηδενική υπόθεση, και πόσοι βαθμοί ελευθερίας έχουν το σύστημα

2) Δοκιμή χωρίς γραμμικό καλοσύνη-of-fit: Ποιες άλλες εξετάσεις μπορώ να κάνω αίτηση για την επικύρωση του μοντέλου;

Δημοσιεύθηκε 03/08/2017 στις 02:28
πηγή χρήστη
Σε άλλες γλώσσες...                            

Cookies help us deliver our services. By using our services, you agree to our use of cookies. Learn more