nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

422976_fr-FR

liliconcom

Un modèle de série chronologique (ARIMA saisonnier) pour faire varier la variabilité selon les saisons

Existe-t-il un modèle ARIMA de série chronologique qui peut représenter une variation accrue (élargir les intervalles de confiance) des valeurs au cours d'une saison ?

À titre d'exemple, ci-joint un ensemble de données indiquant le nombre de décès aux États-Unis.Une analyse des séries chronologiques des décès avant mars 2020 (paramètres ci-dessous) donne une prédiction raisonnable des décès, mais sa précision est réduite de décembre à février (saison de la grippe et mort de l'hiver).Les barres d'erreur autour de ce modèle n'augmentent pas au cours de cette saison, donc alors que relativement peu de points sur l'ensemble du graphique sont en dehors de l'IC, 5 années sur 6 ont des points hors de la plage attendue en décembre ou en janvier, et aucun point n'est hors de gamme pendant toute autre saison.

Je pense que l'intervalle de confiance devrait être plus petit en été et plus large en hiver.Existe-t-il un meilleur modèle pour représenter ces données ? Une simple régression linéaire des valeurs de la semaine x des années précédentes montrerait cette variation accrue, mais j'ai l'impression que quelque chose dans la plate-forme de séries chronologiques devrait pouvoir aider.

Consultez les intervalles de confiance (lignes orange) au bas du graphique ci-dessous et notez qu'avant 2020, les points bleus (décès réels) sont souvent en dehors de cette plage vers le début de chaque année.

Paramètres du modèle de série chronologique utilisé pour tracer ce graphique :

Time Series(
 X( :Week Ending Date ),
 Y( :All Cause ),
 Autocorrelation Lags( 52 ),
 Forecast Periods( 52 ),
 Seasonal ARIMA( 0, 0, 0, 1, 1, 1, 52, Prediction Interval( 0.99 ) ),
 SendToReport(
  Dispatch( {}, "Time Series Basic Diagnostics", OutlineBox, {Close( 1 )} ),
  Dispatch( {"Model Comparison"}, "", ScrollBox, {Background Color( 117 )} )
 )
)

Vous pouvez recréer cette table de données, ou la voir pour un état individuel, en utilisant le script ci-dessous.Téléchargez les données (exportez au format csv et enregistrez-les sur votre ordinateur) à partir d'ici :

https://data.cdc.gov/NCHS/Weekly-Counts-of-Deaths-by-State-and-Select-Causes/3yf8-kanr
https://data.cdc.gov/NCHS/Weekly-Provisional-Counts-of-Deaths-by-State-and-S/muzy-jte6

Names default to here(1);

//Choose Juristiction , either 'United States' or a state name, match a value in that column in csv files.
Juristiction Filter = "United States";
//Juristiction Filter = "Illinois";

//Export data to csv files from the sites below:
//  https://data.cdc.gov/NCHS/Weekly-Counts-of-Deaths-by-State-and-Select-Causes/3yf8-kanr
//  https://data.cdc.gov/NCHS/Weekly-Provisional-Counts-of-Deaths-by-State-and-S/muzy-jte6

dtFinal = Open(
 pick file("Open final file (Weekly counts of deaths by states)", "$Downloads"),
 columns(
  New Column( "Jurisdiction of Occurrence", Character, "Nominal" ),
  New Column( "MMWR Year", Numeric, "Continuous", Format( "Best", 12 ) ),
  New Column( "MMWR Week", Numeric, "Continuous", Format( "Best", 12 ) ),
  New Column( "Week Ending Date",
   Numeric,
   "Continuous",
   Format( "m/d/y", 10 ),
   Input Format( "m/d/y" )
  ),
  New Column( "All  Cause", Numeric, "Continuous", Format( "Best", 12 ) ),
  New Column( "Natural Cause", Numeric, "Continuous", Format( "Best", 12 ) ),
  New Column( "Septicemia (A40-A41)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Malignant neoplasms (C00-C97)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Diabetes mellitus (E10-E14)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Alzheimer disease (G30)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Influenza and pneumonia (J10-J18)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Chronic lower respiratory diseases (J40-J47)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column(
   "Other diseases of respiratory system (J00-J06,J30-J39,J67,J70-J98)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column(
   "Nephritis, nephrotic syndrome and nephrosis (N00-N07,N17-N19,N25-N27)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column(
   "Symptoms, signs and abnormal clinical and laboratory findings, not elsewhere classified (R00-R99)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Diseases of heart (I00-I09,I11,I13,I20-I51)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Cerebrovascular diseases (I60-I69)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "flag_allcause", Character, "Nominal" ),
  New Column( "flag_natcause", Character, "Nominal" ),
  New Column( "flag_sept", Character, "Nominal" ),
  New Column( "flag_neopl", Character, "Nominal" ),
  New Column( "flag_diab", Character, "Nominal" ),
  New Column( "flag_alz", Character, "Nominal" ),
  New Column( "flag_inflpn", Character, "Nominal" ),
  New Column( "flag_clrd", Character, "Nominal" ),
  New Column( "flag_otherresp", Character, "Nominal" ),
  New Column( "flag_nephr", Character, "Nominal" ),
  New Column( "flag_otherunk", Character, "Nominal" ),
  New Column( "flag_hd", Character, "Nominal" ),
  New Column( "flag_stroke", Character, "Nominal" )
 ),
 Import Settings(
  End Of Line( CRLF, CR, LF ),
  End Of Field( Comma, CSV( 1 ) ),
  Strip Quotes( 1 ),
  Use Apostrophe as Quotation Mark( 0 ),
  Use Regional Settings( 0 ),
  Scan Whole File( 1 ),
  Treat empty columns as numeric( 0 ),
  CompressNumericColumns( 0 ),
  CompressCharacterColumns( 0 ),
  CompressAllowListCheck( 0 ),
  Labels( 1 ),
  Column Names Start( 1 ),
  Data Starts( 2 ),
  Lines To Read( "All" ),
  Year Rule( "20xx" )
 )
);

dtProv = Open(
 Pick File("Open provisional file (weekly provisional counts of deaths)", "$Downloads"),
 columns(
  New Column( "Data As Of",
   Numeric,
   "Continuous",
   Format( "m/d/y", 10 ),
   Input Format( "m/d/y" )
  ),
  New Column( "Jurisdiction of Occurrence", Character, "Nominal" ),
  New Column( "MMWR Year", Numeric, "Continuous", Format( "Best", 12 ) ),
  New Column( "MMWR Week", Numeric, "Continuous", Format( "Best", 12 ) ),
  New Column( "Week Ending Date",
   Numeric,
   "Continuous",
   Format( "yyyy-mm-dd", 10 ),
   Input Format( "yyyy-mm-dd" )
  ),
  New Column( "All Cause", Numeric, "Continuous", Format( "Best", 12 ) ),
  New Column( "Natural Cause", Numeric, "Continuous", Format( "Best", 12 ) ),
  New Column( "Septicemia (A40-A41)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Malignant neoplasms (C00-C97)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Diabetes mellitus (E10-E14)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Alzheimer disease (G30)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Influenza and pneumonia (J09-J18)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Chronic lower respiratory diseases (J40-J47)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column(
   "Other diseases of respiratory system (J00-J06,J30-J39,J67,J70-J98)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column(
   "Nephritis, nephrotic syndrome and nephrosis (N00-N07,N17-N19,N25-N27)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column(
   "Symptoms, signs and abnormal clinical and laboratory findings, not elsewhere classified (R00-R99)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Diseases of heart (I00-I09,I11,I13,I20-I51)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "Cerebrovascular diseases (I60-I69)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "COVID-19 (U071, Multiple Cause of Death)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "COVID-19 (U071, Underlying Cause of Death)",
   Numeric,
   "Continuous",
   Format( "Best", 12 )
  ),
  New Column( "flag_allcause", Character, "Nominal" ),
  New Column( "flag_natcause", Character, "Nominal" ),
  New Column( "flag_sept", Character, "Nominal" ),
  New Column( "flag_neopl", Character, "Nominal" ),
  New Column( "flag_diab", Character, "Nominal" ),
  New Column( "flag_alz", Character, "Nominal" ),
  New Column( "flag_inflpn", Character, "Nominal" ),
  New Column( "flag_clrd", Character, "Nominal" ),
  New Column( "flag_otherresp", Character, "Nominal" ),
  New Column( "flag_nephr", Character, "Nominal" ),
  New Column( "flag_otherunk", Character, "Nominal" ),
  New Column( "flag_hd", Character, "Nominal" ),
  New Column( "flag_stroke", Character, "Nominal" ),
  New Column( "flag_cov19mcod", Character, "Nominal" ),
  New Column( "flag_cov19ucod", Character, "Nominal" )
 ),
 Import Settings(
  End Of Line( CRLF, CR, LF ),
  End Of Field( Comma, CSV( 1 ) ),
  Strip Quotes( 1 ),
  Use Apostrophe as Quotation Mark( 0 ),
  Use Regional Settings( 0 ),
  Scan Whole File( 1 ),
  Treat empty columns as numeric( 0 ),
  CompressNumericColumns( 0 ),
  CompressCharacterColumns( 0 ),
  CompressAllowListCheck( 0 ),
  Labels( 1 ),
  Column Names Start( 1 ),
  Data Starts( 2 ),
  Lines To Read( "All" ),
  Year Rule( "20xx" )
 )
);

dtProv:" All Cause"n << Set Name("All Cause");
dtFinal:"All  Cause"n << Set Name("All Cause");

dtAll = dtFinal << Concatenate(
 dtProv,
 Create source column
);

dtAll << Select Where(:Jurisdiction of Occurrence == Juristiction Filter );

dtSub = dtAll << Subset(
 Selected Rows( 1 ),
 Selected columns only( 0 )
);

dtSub << Select where( :Week Ending Date > date dmy( 1,3,2020) );

//Hide new data
dtSub << Hide and Exclude(1);

ts = dtSub << Time Series(
 X( :Week Ending Date ),
 Y( :All Cause ),
 Autocorrelation Lags( 52 ),
 Forecast Periods( 52 ),
 Seasonal ARIMA( 0, 0, 0, 1, 1, 1, 52, Prediction Interval( 0.99 ) ),
 SendToReport(
  Dispatch( {}, "Time Series Basic Diagnostics", OutlineBox, {Close( 1 )} ),
  Dispatch( {"Model Comparison"}, "", ScrollBox, {Background Color( 117 )} )
 )
);

dtPred = (((ts << XPath("//OutlineBox[@helpKey='ARIMA Model']"))[1]) << Get Scriptable object) << Save Prediction Formula;

dtJoined = dtSub << Join(
 With( dtPred ),
 Merge Same Name Columns,
 By Row Number
);

dtJoined << New Column("Excess v Expected", Numeric, "Continuous", Format("Best", 12), Formula(:All Cause - :Predicted All Cause), Set Display Width(100));
dtJoined << New Column("Cumulative[Excess v Expected]", Numeric, "Continuous", Format("Best", 12), Formula(Col Cumulative Sum(:Excess v Expected)), Set Display Width(173));

last row with actuals = max(dtJoined << Get rows where( !Is Missing( :All Cause ) ))-1;
last date with actuals = dtJoined:Week Ending Date[last row with actuals];

dtJoined << Graph Builder(
 Size( 1174, 786 ),
 Show Control Panel( 0 ),
 Show Subtitle( 1 ),
 Variables(
  X( :Week Ending Date ),
  Y( :"Cumulative[Excess v Expected]"n ),
  Y( :Excess v Expected ),
  Y( :All Cause ),
  Y( :Predicted All Cause, Position( 3 ) ),
  Y( :"Upper CL (0.99) All Cause"n, Position( 3 ) ),
  Y( :"Lower CL (0.99) All Cause"n, Position( 3 ) )
 ),
 Relative Sizes( "Y", [229 204 387] ),
 Elements( Position( 1, 1 ), Points( X, Y, Legend( 18 ) ) ),
 Elements(
  Position( 1, 2 ),
  Points( X, Y, Legend( 20 ) ),
  Area( X, Y, Legend( 22 ), Interval Style( "Band" ) )
 ),
 Elements(
  Position( 1, 3 ),
  Points( X, Y( 1 ), Legend( 13 ) ),
  Line( X, Y( 2 ), Y( 3 ), Y( 4 ), Legend( 15 ) )
 ),
 Local Data Filter(
  Add Filter(
   columns( :Week Ending Date ),
   Where(
    :Week Ending Date >= 3471638400 & :Week Ending Date <= last date with actuals
   )
  )
 ),
 SendToReport(
  Dispatch(
   {},
   "Week Ending Date",
   ScaleBox,
   {Min( 3462203520 ), Max( last date with actuals + in days(7) ), Interval( "Year" ),
   Inc( 1 ), Minor Ticks( 0 ), Add Ref Line(
    last date with actuals, "Solid", "Black", Format( last date with actuals, "ddMonyyyy"), 1
   ), Label Row( Show Major Grid( 1 ) )}
  ),
  Dispatch(
   {},
   "Cumulative[Excess v Expected]",
   ScaleBox,
   {Min( 0 ), //Max( 819328.125380856 ), Inc( 200000 ),
   Minor Ticks( 1 ), Add Ref Line( 0, "Solid", "Black", "", 1 ),
   Label Row( Show Major Grid( 1 ) )}
  ),
  Dispatch(
   {},
   "Excess v Expected",
   ScaleBox,
   {Min( 0 ), //Max( 26451.6956907074 ), Inc( 10000 ),
   Minor Ticks( 1 ), Add Ref Line( 0, "Solid", "Black", "", 1 ),
   Label Row( {Show Major Grid( 1 ), Show Minor Grid( 1 )} )}
  ),
  Dispatch(
   {},
   "All Cause",
   ScaleBox,
   {Min( 0 ), //Max( 90470.7642581624 ), Inc( 10000 ),
   Minor Ticks( 4 ), Label Row(
    {Show Major Grid( 1 ), Show Minor Grid( 1 )}
   )}
  ),
  Dispatch(
   {},
   "400",
   ScaleBox,
   {Legend Model(
    22,
    Properties(
     0,
     {Fill Color( 20 ), Transparency( 0.5 )},
     Item ID( "Excess v Expected", 1 )
    )
   ), Legend Model(
    13,
    Properties( 0, {Line Color( 21 )}, Item ID( "All Cause", 1 ) )
   ), Legend Model(
    15,
    Properties(
     0,
     {Line Color( 32 )},
     Item ID( "Predicted All Cause", 1 )
    ),
    Properties(
     1,
     {Line Color( 70 )},
     Item ID( "Upper CL (0.99) All Cause", 1 )
    ),
    Properties(
     2,
     {Line Color( 70 )},
     Item ID( "Lower CL (0.99) All Cause", 1 )
    )
   )}
  ),
  Dispatch(
   {},
   "graph title",
   TextEditBox,
   {Set Text( Juristiction Filter || " Actual Deaths v Expected Deaths (Excess Deaths)" )}
  ),
  Dispatch(
   {},
   "graph 1 title",
   TextEditBox,
   {Set Text(
    "Source: https://data.cdc.gov/NCHS/Weekly-Counts-of-Deaths-by-State-and-Select-Causes/3yf8-kanr  and   https://data.cdc.gov/NCHS/Weekly-Provisional-Counts-of-Deaths-by-State-and-S/muzy-jte6"
   )}
  ),
  Dispatch( {}, "X title", TextEditBox, {Set Text( "" )} ),
  Dispatch(
   {},
   "Y title",
   TextEditBox,
   {Set Text( "Cumulative Excess Deaths" )}
  ),
  Dispatch(
   {},
   "Y 1 title",
   TextEditBox,
   {Set Text( "Excess v Expected\!NDeaths per Month" )}
  ),
  Dispatch( {}, "Y 2 title", TextEditBox, {Set Text( "Deaths per Month" )} ),
  Dispatch(
   {},
   "400",
   LegendBox,
   {Legend Position(
    {18, [0], 20, [1], 22, [6, -3], 13, [2], 15, [3, 4, 5]}
   )}
  ),
  Dispatch(
   {},
   "Graph Builder",
   FrameBox(3),
   {Add Graphics Script(
    3,
    Description( "Script" ),
    dt = Current Data Table();

    yMatrix = (dtJoined:All Cause << get values)[54::last row with actuals] |/
    Matrix( Reverse( As List( (dtJoined:Predicted All Cause << get values)[54::last row with actuals] ) ) );
    xMatrix = (dtJoined:Week Ending Date << get values)[54::last row with actuals] |/
    Matrix( Reverse( As List( (dtJoined:Week Ending Date << get values)[54::last row with actuals] ) ) );
    Fill Color( "Medium Dark Green" );
    Transparency( 0.5 );
    Polygon( xMatrix, yMatrix );
   ), Grid Line Order( 1 ), Reference Line Order( 3 )}
  )
 )
);

dtFinal << Close Window;
dtProv << Close Window;
dtAll << Close Window;
dtSub << Close Window;
dtPred << Close Window;

Les noms par défaut sont ici (1) ; //Choisissez Juristiction , soit « États-Unis » ou un nom d'État, correspond à une valeur dans cette colonne dans les fichiers csv. Filtre de juridiction = « États-Unis » ; //Filtre de juridiction = "Illinois"; //Exporter les données vers des fichiers csv à partir des sites ci-dessous : // https://data.cdc.gov/NCHS/Weekly-Counts-of-Deaths-by-State-and-Select-Causes/3yf8-kanr // https

/data.cdc.gov/NCHS/Weekly-Provisional-Counts-of-Deaths-by-State-and-S/muzy-jte6 dtFinal = Open( pick file ("Ouvrir le fichier final (nombre hebdomadaire de décès par états )", "$Téléchargements"), colonnes( Nouvelle colonne( "Juridiction d'occurrence", Caractère, "Nominal" ), Nouvelle colonne( "Année MMWR", Numérique, "Continu", Format( "Meilleur", 12 ) ) , Nouvelle colonne( "Semaine MMWR", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Date de fin de semaine", Numérique, "Continu", Format( "m/j/a", 10 ), Format d'entrée( "m/d/y" ) ), Nouvelle colonne( "Toutes causes", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Cause naturelle", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Septicémie (A40-A41)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Tumeurs malignes (C00 -C97)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Diabète sucré us (E10-E14)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Maladie d'Alzheimer (G30)", Numérique, "Continu", Format( "Meilleur", 12 ) ) , Nouvelle colonne( "Grippe et pneumonie (J10-J18)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Maladies chroniques des voies respiratoires inférieures (J40-J47)", Numérique, "Continu ", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Autres maladies du système respiratoire (J00-J06,J30-J39,J67,J70-J98)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Néphrite, syndrome néphrotique et néphrose (N00-N07,N17-N19,N25-N27)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Symptômes , signes et résultats anormaux cliniques et de laboratoire, non classés ailleurs (R00-R99)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Maladies cardiaques (I00-I09,I11,I13) ,I20-I51)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Maladies cérébrovasculaires (I60-I69)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "flag_allcaus e", Caractère, "Nominal" ), Nouvelle colonne( "flag_natcause", Caractère, "Nominal" ), Nouvelle colonne( "flag_sept", Caractère, "Nominal" ), Nouvelle colonne( "flag_neopl", Caractère, "Nominal" ), Nouvelle colonne( "flag_diab", Caractère, "Nominal" ), Nouvelle colonne( "flag_alz", Caractère, "Nominal" ), Nouvelle colonne( "flag_inflpn", Caractère, "Nominal" ), Nouvelle colonne( "flag_clrd" , Caractère, "Nominal" ), Nouvelle colonne( "flag_otherresp", Caractère, "Nominal" ), Nouvelle colonne( "flag_nephr", Caractère, "Nominal" ), Nouvelle colonne( "flag_otherunk", Caractère, "Nominal" ), Nouvelle colonne( "flag_hd", Caractère, "Nominal" ), Nouvelle colonne( "flag_stroke", Caractère, "Nominal" ) ), Paramètres d'importation ( Fin de ligne ( CRLF, CR, LF ), Fin de champ ( Virgule, CSV ( 1 ) ), Supprimer les guillemets ( 1 ), Utiliser l'apostrophe comme guillemet ( 0 ), Utiliser les paramètres régionaux ( 0 ), Analyser le fichier entier ( 1 ), Traiter les colonnes vides comme numériques ( 0 ), CompressNumericColumns ( 0 ), CompressCharacterColumns ( 0 ), CompressAllowListCheck( 0 ), Labels( 1 ), Début des noms de colonnes( 1 ), Début des données( 2 ), Lignes à lire( "Tous" ), Règle de l'année( "20xx" ) ) ); dtProv = Open( Pick File("Ouvrir le fichier provisoire (nombre provisoire hebdomadaire de décès)", "$Téléchargements"), colonnes( Nouvelle colonne( "Données à partir du", Numérique, "Continu", Format( "m/d/ y", 10 ), Format d'entrée( "m/d/y" ) ), Nouvelle colonne( "Juridiction of Occurrence", Caractère, "Nominal" ), Nouvelle colonne( "MMWR Year", Numérique, "Continu", Format ( "Meilleur", 12 ) ), Nouvelle colonne( "Semaine MMWR", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Date de fin de semaine", Numérique, "Continu", Format( "aaaa-mm-jj", 10 ), Format d'entrée( "aaaa-mm-jj" ) ), Nouvelle colonne( "Toutes causes", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne ( "Cause naturelle", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Septicémie (A40-A41)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Tumeurs malignes (C00-C97)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( " Diabète sucré (E10-E14)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Maladie d'Alzheimer (G30)", Numérique, "Co ntinuous", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Grippe et pneumonie (J09-J18)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Chronique des voies respiratoires inférieures maladies (J40-J47)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Autres maladies du système respiratoire (J00-J06,J30-J39,J67,J70-J98)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Néphrite, syndrome néphrotique et néphrose (N00-N07,N17-N19,N25-N27)", Numérique, "Continu", Format( " Meilleur", 12 ) ), Nouvelle colonne( "Symptômes, signes et résultats anormaux cliniques et de laboratoire, non classés ailleurs (R00-R99)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Maladies cardiaques (I00-I09,I11,I13,I20-I51)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "Maladies cérébrovasculaires (I60-I69)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "COVID-19 (U071, Cause de décès multiple)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "COVID-19 (U071, Underl ying Cause of Death)", Numérique, "Continu", Format( "Meilleur", 12 ) ), Nouvelle colonne( "flag_allcause", Caractère, "Nominal" ), Nouvelle colonne( "flag_natcause", Caractère, "Nominal" ) , Nouvelle colonne( "flag_sept", Caractère, "Nominal" ), Nouvelle colonne( "flag_neopl", Caractère, "Nominal" ), Nouvelle colonne( "flag_diab", Caractère, "Nominal" ), Nouvelle colonne( "flag_alz", Caractère, "Nominal" ), Nouvelle colonne( "flag_inflpn", Caractère, "Nominal" ), Nouvelle colonne( "flag_clrd", Caractère, "Nominal" ), Nouvelle colonne( "flag_otherresp", Caractère, "Nominal" ), Nouveau Colonne( "flag_nephr", Caractère, "Nominal" ), Nouvelle colonne( "flag_otherunk", Caractère, "Nominal" ), Nouvelle colonne( "flag_hd", Caractère, "Nominal" ), Nouvelle colonne( "flag_stroke", Caractère, "Nominal" ), Nouvelle colonne( "flag_cov19mcod", Caractère, "Nominal" ), Nouvelle colonne( "flag_cov19ucod", Caractère, "Nominal" ) ), Paramètres d'importation( Fin de ligne( CRLF, CR, LF ), Fin de Champ( Virgule, CSV( 1 ) ), Supprimer les guillemets( 1 ), Utiliser l'apostrophe comme guillemet( 0 ), Utiliser les paramètres régionaux( 0 ) , Scan Whole File( 1 ), Traiter les colonnes vides comme numérique( 0 ), CompressNumericColumns( 0 ), CompressCharacterColumns( 0 ), CompressAllowListCheck( 0 ), Labels( 1 ), Column Names Start( 1 ), Data Starts( 2 ), Lignes à lire( "Tous" ), Règle de l'année( "20xx" ) ) ); dtProv : "Toute cause"n << Définir le nom("Toute cause"); dtFinal:"Toute cause"n << Définir le nom("Toute cause"); dtAll = dtFinal << Concaténer( dtProv, Créer une colonne source ); dtAll << Sélectionnez Où(:Juridiction d'occurrence == Filtre de juridiction ); dtSub = dtAll << Sous-ensemble ( lignes sélectionnées ( 1 ), colonnes sélectionnées uniquement ( 0 ) ); dtSub << Sélectionnez où ( :Date de fin de semaine > date dmy( 1,3,2020) ); //Masquer les nouvelles données dtSub << Masquer et exclure (1) ; ts = dtSub << Séries temporelles( X( :Date de fin de semaine ), Y( :Toutes causes ), Délais d'autocorrélation( 52 ), Périodes de prévision( 52 ), ARIMA saisonnier( 0, 0, 0, 1, 1, 1, 52, Intervalle de prédiction (0.99 ) ), SendToReport( Dispatch( {}, "Diagnostics de base des séries temporelles", OutlineBox, {Close( 1 )} ), Dispatch( {"Comparaison de modèles"}, "", ScrollBox, {Couleur d'arrière-plan ( 117 )} ) ) ); dtPred = (((ts << XPath("//OutlineBox[@helpKey='ARIMA Model']"))[1]) << Obtenir un objet scriptable) << Enregistrer la formule de prédiction ; dtJoined = dtSub << Join( With( dtPred ), Fusionner les colonnes de même nom, Par numéro de ligne ); dtJoined << Nouvelle colonne("Excès v attendu", numérique, "Continu", Format("Meilleur", 12), Formule(:Toute cause - :Prédite toute cause), Définir la largeur d'affichage (100)); dtJoined << Nouvelle colonne("Cumulatif[Excès v Attendu]", Numérique, "Continu", Format("Meilleur", 12), Formule(Col Cumulative Sum(:Excess v Attendu)), Définir la largeur d'affichage (173)) ; dernière ligne avec actuals = max(dtJoined << Récupère les lignes où( !Is Missing( :All Cause ) ))-1; dernière date avec les chiffres réels = dtJoined:Week Ending Date [dernière ligne avec les chiffres réels] ; dtJoined << Graph Builder( Size( 1174, 786 ), Show Control Panel( 0 ), Show Subtitle( 1 ), Variables( X( :Week Ending Date ), Y( :"Cumulative[Excess v Expected]"n ), Y( :Excès v attendu ), Y( :Toute cause ), Y( :Toute cause prédite, Position( 3 ) ), Y( :"Upper CL (0,99)Toute cause"n, Position( 3 ) ), Y( :"Lower CL (0,99) toutes causes"n, Position( 3 ) ) ), Tailles relatives( "Y", [229 204 387] ), Elements( Position( 1, 1 ), Points( X, Y, Legend( 18 ) ) ), Éléments( Position( 1, 2 ), Points( X, Y, Légende( 20 ) ), Aire( X, Y, Légende( 22 ), Style d'intervalle( "Bande" ) ) ), Éléments( Position ( 1, 3 ), Points( X, Y( 1 ), Légende( 13 ) ), Ligne( X, Y( 2 ), Y( 3 ), Y( 4 ), Légende( 15 ) ) ), Filtre de données local ( Ajouter un filtre( colonnes( :Date de fin de semaine ), Où( :Date de fin de semaine >= 3471638400 & :Date de fin de semaine <= dernière date avec les chiffres réels ) ) ), SendToReport( Dispatch( {}, "Date de fin de semaine", ScaleBox , {Min( 3462203520 ), Max( dernière date avec chiffres réels + en jours(7) ), Interval( "Année" ), Inc( 1 ), Minor Ticks( 0 ), Ajouter une ligne de référence( dernière date avec les chiffres réels, "Solide", "Noir", Format( dernière date avec les chiffres réels, "ddMonyyyy"), 1 ), Ligne d'étiquette( Afficher la grille principale( 1 ) )} ), Dispatch( { }, "Cumulative[Excess v Expected]", ScaleBox, {Min( 0 ), //Max( 819328.125380856 ), Inc( 200000 ), Minor Ticks( 1 ), Add Ref Line( 0, "Solid", "Black" , "", 1 ), Label Row( Show Major Grid( 1 ) )} ), Dispatch( {}, "Excess v Expected", ScaleBox, {Min( 0 ), //Max( 26451.6956907074 ), Inc( 10000 ) , Coches mineures( 1 ), Ajouter une ligne de référence( 0, "Solide", "Noir", "", 1 ), Ligne d'étiquette( {Afficher la grille principale( 1 ), Afficher la grille mineure( 1 )} )} ), Répartir ( {}, "All Cause", ScaleBox, {Min( 0 ), //Max( 90470.7642581624 ), Inc( 10000 ), Minor Ticks( 4 ), Label Row( {Afficher la grille principale( 1 ), Afficher la grille secondaire( 1 )} )} ), Dispatch( {}, "400", ScaleBox, {Legend Model( 22, Properties( 0, {Fill Color( 20 ), Transparency( 0.5 )}, Item ID( "Excess v Expected", 1 ) ) ), Legend Model( 13, Properties( 0, {Line Color( 21 )}, Item ID( "All Cause", 1 ) ) ), Legend Model( 15, Prope rties( 0, {Couleur de la ligne( 32 )}, ID de l'article( "Toute cause prévue", 1 ) ), Propriétés( 1, {Couleur de la ligne( 70 )}, ID de l'article( "Upper CL (0,99) Toutes les causes", 1 ) ), Properties( 2, {Line Color( 70 )}, Item ID( "Lower CL (0,99) All Cause", 1 ) ) )} ), Dispatch( {}, "graph title", TextEditBox, {Set Texte( Filtre de juridiction || " Décès réels contre Décès attendus (décès excessifs)" )} ), Dispatch( {}, "titre du graphique 1", TextEditBox, {Set Text( "Source : https://data.cdc.gov/NCHS/Weekly-Counts -of-Deaths-by-State-and-Select-Causes/3yf8-kanr et https://data.cdc.gov/NCHS/Weekly-Provisional-Counts-of-Deaths-by-State-and-S/muzy -jte6" )} ), Dispatch( {}, "X title", TextEditBox, {Set Text( "" )} ), Dispatch( {}, "Y title", TextEditBox, {Set Text( "Cumulative Excess Deaths" )} ), Dispatch( {}, "Titre Y 1", TextEditBox, {Set Text( "Excess v Expected\!NDeaths per Month" )} ), Dispatch( {}, "Titre Y 2", TextEditBox, {Set Text( "Décès par mois" )} ), Dispatch( {}, "400", LegendBox, {Legend Position( {18, [0], 20, [1], 22, [6, -3], 13, [2], 15, [3, 4, 5]} )} ), Dispatch( {}, "Graph Builder", FrameBox(3), {Add Graphics Script( 3, Description( "Script" ), dt = Current Table de données (); yMatrix = (dtJoined:All Cause << get values)[54::last row with actuals] |/ Matrix( Reverse( As List( (dtJoined:Predicted All Cause << get values)[54: :dernière ligne avec les chiffres réels] ) ) ); xMatrix = (dtJoined:Week Ending Date << obtenir les valeurs)[54::dernière ligne avec les chiffres réels] |/ Matrix( Reverse( As List( (dtJoined:Week Ending Date << obtenir les valeurs)[54::dernière ligne avec les chiffres réels ] ) ) ); Couleur de remplissage( "Vert foncé moyen" ); Transparence( 0,5 ); Polygon( xMatrice, yMatrice ); ), Ordre des lignes de la grille ( 1 ), Ordre des lignes de référence ( 3 )} ) ) ); dtFinal << Fermer la fenêtre ; dtProv << Fermer la fenêtre ; dtAll << Fermer la fenêtre; dtSub << Fermer la fenêtre ; dtPred << Fermer la fenêtre ;

Analytique

Re : Un modèle de série chronologique (ARIMA saisonnier) pour une variabilité variable selon les saisons

Je ne peux pas être sûr qu'il s'agisse d'un problème de variation ou de biais ou des deux. Si les valeurs aberrantes sont toujours ou presque toujours proches de la partie supérieure, cela ressemble davantage à un problème de biais. Le modèle de fonction de transfert dans la plate-forme Time Series peut être utile. Par exemple, utilisez la semaine de l'année ou du mois ou tout autre indicateur de résolution de temps approprié comme variable d'entrée.
S'il s'agit bien d'un problème de variation, certains peuvent être traités dans la plate-forme. La modélisation de séries chronologiques transformées de Box-Cox, ajoutée pendant JMP16, est un outil. Mais ici, peut-être plus simple. Je vois que le modèle souhaité n'a pas de commandes non saisonnières. Cela devrait être techniquement équivalent à l'ajustement de 52 modèles distincts par semaine de l'année, mais en limitant la même variance, je parie. Si vous ajustez simplement 52 modèles séparément, vous obtenez des variations différentes. Mais assurez-vous d'avoir toujours 52 points de données par an.

Re : Un modèle de série chronologique (ARIMA saisonnier) pour une variabilité variable au cours des différentes saisons

Je retire les deux dernières phrases sur le montage de modèles séparés. Je n'avais pas réalisé que les coefficients devaient également être limités entre les modèles. Mais il semble que la transformation de Box-Cox soit pertinente, car la variabilité plus élevée est liée à une valeur de série plus élevée.

Comments

Quick Links