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_es-ES

liliconcom

Un modelo de serie temporal (ARIMA estacional) para variar la variabilidad en diferentes estaciones.

¿Existe un modelo ARIMA de serie temporal que pueda representar una mayor variación (ampliar los intervalos de confianza) en los valores durante una temporada?

Como ejemplo, se adjunta un conjunto de datos que muestra el número de muertes en los Estados Unidos.Un análisis de series de tiempo de muertes antes de marzo de 2020 (configuración a continuación) brinda una predicción razonable de muertes, sin embargo, su precisión se reduce durante diciembre-febrero (temporada de gripe y final del invierno).Las barras de error alrededor de este modelo no aumentan durante esa temporada, por lo que aunque relativamente pocos puntos en todo el gráfico están fuera del IC, 5 de cada 6 años tienen puntos fuera del rango esperado en diciembre o enero, y ningún punto está fuera de rango durante cualquier otra temporada.

Creo que el intervalo de confianza debería ser menor en verano y más amplio en invierno.¿Existe un modelo mejor para representar estos datos? Una simple regresión lineal de los valores de la semana x de años anteriores mostraría esta mayor variación, pero creo que algo en la plataforma de series de tiempo debería poder ayudar.

Vea los intervalos de confianza (líneas naranjas) en la parte inferior del cuadro a continuación y observe cómo antes de 2020 los puntos azules (muertes reales) a menudo están fuera de ese rango cerca del comienzo de cada año.

Configuración para el modelo de series de tiempo utilizado para dibujar este gráfico:

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 )} )
 )
)

Puede volver a crear esta tabla de datos, o verla para un estado individual, utilizando la secuencia de comandos a continuación.Descargue los datos (expórtelos como csv y guárdelos en su computadora) desde aquí:

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;

Los nombres están predeterminados aquí (1); // Elija Juristiction, ya sea 'Estados Unidos' o el nombre de un estado, haga coincidir un valor en esa columna en los archivos csv. Filtro de jurisdicción = "Estados Unidos"; // Filtro de jurisdicción = "Illinois"; // Exportar datos a archivos csv desde los siguientes sitios: // 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 (Recuentos semanales de muertes por estados ) "," $ Descargas "), columnas (Nueva columna (" Jurisdicción del suceso ", Carácter," Nominal "), Nueva columna (" Año MMWR ", Numérico," Continuo ", Formato (" Mejor ", 12)) , Nueva columna ("Semana MMWR", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("Fecha de finalización de la semana", Numérico, "Continuo", Formato ("m / d / y", 10), Formato de entrada ("m / d / y")), Nueva columna ("Todas las causas", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("Causa natural", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("Septicemia (A40-A41)", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("Neoplasias malignas (C00 -C97) ", Numérico," Continuo ", Formato (" Mejor ", 12)), Nueva columna (" Diabetes mellit us (E10-E14) ", Numérico," Continuo ", Formato (" Mejor ", 12)), Nueva columna (" Enfermedad de Alzheimer (G30) ", Numérico," Continuo ", Formato (" Mejor ", 12)) , Nueva columna ("Influenza y neumonía (J10-J18)", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("Enfermedades crónicas de las vías respiratorias inferiores (J40-J47)", Numérico, "Continuo ", Formato (" Mejor ", 12)), Nueva columna (" Otras enfermedades del sistema respiratorio (J00-J06, J30-J39, J67, J70-J98) ", Numérico," Continuo ", Formato (" Mejor ", 12)), Nueva columna ("Nefritis, síndrome nefrótico y nefrosis (N00-N07, N17-N19, N25-N27)", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("Síntomas , signos y hallazgos clínicos y de laboratorio anormales, no clasificados en otra parte (R00-R99) ", Numérico," Continuo ", Formato (" Mejor ", 12)), Nueva columna (" Enfermedades del corazón (I00-I09, I11, I13 , I20-I51) ", Numérico," Continuo ", Formato (" Mejor ", 12)), Nueva columna (" Enfermedades cerebrovasculares (I60-I69) ", Numérico," Continuo ", Formato (" Mejor ", 12) ), Nueva columna ("flag_allcaus e ", Carácter," Nominal "), Nueva columna (" flag_natcause ", Carácter," Nominal "), Nueva columna (" flag_sept ", Carácter," Nominal "), Nueva columna (" flag_neopl ", Carácter," Nominal " ), Nueva columna ("flag_diab", Carácter, "Nominal"), Nueva columna ("flag_alz", Carácter, "Nominal"), Nueva columna ("flag_inflpn", Carácter, "Nominal"), Nueva columna ("flag_clrd" , Carácter, "Nominal"), Nueva columna ("flag_otherresp", Carácter, "Nominal"), Nueva columna ("flag_nephr", Carácter, "Nominal"), Nueva columna ("flag_otherunk", Carácter, "Nominal"), Nueva columna ("flag_hd", carácter, "nominal"), nueva columna ("flag_stroke", carácter, "nominal")), configuración de importación (final de línea (CRLF, CR, LF), final de campo (coma, CSV) (1)), Quitar comillas (1), Usar apóstrofo como comillas (0), Usar configuración regional (0), Escanear todo el archivo (1), Tratar columnas vacías como numéricas (0), CompressNumericColumns (0), CompressCharacterColumns ( 0), CompressAllowListCheck (0), Etiquetas (1), Inicio de nombres de columna (1), Inicio de datos (2 ), Líneas para leer ("Todas"), Regla del año ("20xx"))); dtProv = Abrir (Seleccionar archivo ("Abrir archivo provisional (recuentos provisionales semanales de muertes)", "$ Descargas"), columnas (Nueva columna ("Datos a partir de", Numérico, "Continuo", Formato ("m / d / y ", 10), Formato de entrada (" m / d / y ")), Nueva columna (" Jurisdicción del suceso ", Carácter," Nominal "), Nueva columna (" Año MMWR ", Numérico," Continuo ", Formato ("Mejor", 12)), Nueva columna ("Semana MMWR", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("Fecha de finalización de la semana", Numérico, "Continuo", Formato ( "aaaa-mm-dd", 10), Formato de entrada ("aaaa-mm-dd")), Nueva columna ("Todas las causas", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("Causa natural", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("Septicemia (A40-A41)", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("Neoplasias malignas (C00-C97)", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("Diabetes mellitus (E10-E14)", Numérico, "Continuo", Formato ( "Mejor", 12)), Nueva columna ("Enfermedad de Alzheimer (G30)", Numérico, "Co n continuo ", Formato (" Mejor ", 12)), Nueva columna (" Influenza y neumonía (J09-J18) ", Numérico," Continuo ", Formato (" Mejor ", 12)), Nueva columna (" Crónico enfermedades (J40-J47) ", Numérico," Continuo ", Formato (" Mejor ", 12)), Nueva columna (" Otras enfermedades del sistema respiratorio (J00-J06, J30-J39, J67, J70-J98) ", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("Nefritis, síndrome nefrótico y nefrosis (N00-N07, N17-N19, N25-N27)", Numérico, "Continuo", Formato (" Mejor ", 12)), Nueva columna (" Síntomas, signos y hallazgos clínicos y de laboratorio anormales, no clasificados en otra parte (R00-R99) ", Numérico," Continuo ", Formato (" Mejor ", 12)), Nueva columna ( "Enfermedades del corazón (I00-I09, I11, I13, I20-I51)", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("Enfermedades cerebrovasculares (I60-I69)", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ("COVID-19 (U071, Causa de muerte múltiple)", Numérico, "Continuo", Formato ("Mejor", 12)), Nueva columna ( "COVID-19 (U071, Underl ying causa de muerte) ", Numérico," Continuo ", Formato (" Mejor ", 12)), Nueva columna (" flag_allcause ", Carácter," Nominal "), Nueva columna (" flag_natcause ", Carácter," Nominal ") , Nueva columna ("flag_sept", Carácter, "Nominal"), Nueva columna ("flag_neopl", Carácter, "Nominal"), Nueva columna ("flag_diab", Carácter, "Nominal"), Nueva columna ("flag_alz", Carácter, "Nominal"), Nueva columna ("flag_inflpn", Carácter, "Nominal"), Nueva columna ("flag_clrd", Carácter, "Nominal"), Nueva columna ("flag_otherresp", Carácter, "Nominal"), Nuevo Columna ("flag_nephr", Carácter, "Nominal"), Nueva columna ("flag_otherunk", Carácter, "Nominal"), Nueva columna ("flag_hd", Carácter, "Nominal"), Nueva columna ("flag_stroke", Carácter, "Nominal"), Nueva columna ("flag_cov19mcod", Carácter, "Nominal"), Nueva columna ("flag_cov19ucod", Carácter, "Nominal")), Configuración de importación (Fin de línea (CRLF, CR, LF), Fin de Campo (coma, CSV (1)), Quitar comillas (1), Usar apóstrofo como comillas (0), Usar configuración regional (0) , Escanear todo el archivo (1), Tratar columnas vacías como numéricas (0), CompressNumericColumns (0), CompressCharacterColumns (0), CompressAllowListCheck (0), Etiquetas (1), Inicio de nombres de columna (1), Inicio de datos (2), Líneas para leer ("Todas"), Regla del año ("20xx"))); dtProv: "Todas las causas" n << Establecer nombre ("Todas las causas"); dtFinal: "Todas las causas" n << Establecer nombre ("Todas las causas"); dtAll = dtFinal << Concatenar (dtProv, Crear columna fuente); dtAll << Seleccione dónde (: Jurisdicción del suceso == Filtro de jurisdicción); dtSub = dtAll << Subconjunto (Filas seleccionadas (1), Solo columnas seleccionadas (0)); dtSub << Seleccione dónde (: Fecha de finalización de la semana> fecha dmy (1,3,2020)); // Ocultar nuevos datos dtSub << Ocultar y excluir (1); ts = dtSub << Serie temporal (X (: Fecha de finalización de la semana), Y (: Todas las causas), Retrasos de autocorrelación (52), Períodos de pronóstico (52), ARIMA estacional (0, 0, 0, 1, 1, 1, 52, Intervalo de predicción (0,99)), SendToReport (Despacho ({}, "Diagnóstico básico de la serie temporal", OutlineBox, {Cerrar (1)}), Despacho ({"Comparación de modelos"}, "", ScrollBox, {Color de fondo (117)}))); dtPred = (((ts << XPath ("// OutlineBox [@ helpKey = 'Modelo ARIMA']")) [1]) << Obtener objeto Scriptable) << Guardar fórmula de predicción; dtJoined = dtSub << Unir (Con (dtPred), Combinar columnas del mismo nombre, Por número de fila); dtJoined << Nueva columna ("Exceso v esperado", Numérico, "Continuo", Formato ("Mejor", 12), Fórmula (: Todas las causas -: Todas las causas previstas), Establecer ancho de visualización (100)); dtJoined << Nueva columna ("Acumulado [Exceso v esperado]", Numérico, "Continuo", Formato ("Mejor", 12), Fórmula (Suma acumulada de columnas (: Exceso v esperado)), Establecer ancho de visualización (173)) ; última fila con datos reales = max (dtJoined << Obtener filas donde (! Falta (: Todas las causas))) - 1; última fecha con datos reales = dtJoined: Fecha de finalización de la semana [última fila con datos reales]; dtJoined << Generador de gráficos (Tamaño (1174, 786), Mostrar panel de control (0), Mostrar subtítulo (1), Variables (X (: Fecha de finalización de la semana), Y (: "Acumulativo [Exceso v esperado]" n), Y (: Exceso v esperado), Y (: Todas las causas), Y (: Todas las causas previstas, Posición (3)), Y (: "CL superior (0,99) Todas las causas" n, Posición (3)), Y ( : "CL inferior (0,99) Todas las causas" n, Posición (3))), Tamaños relativos ("Y", [229 204 387]), Elementos (Posición (1, 1), Puntos (X, Y, Leyenda ( 18))), Elementos (Posición (1, 2), Puntos (X, Y, Leyenda (20)), Área (X, Y, Leyenda (22), Estilo de intervalo ("Banda"))), Elementos (Posición (1, 3), Puntos (X, Y (1), Leyenda (13)), Línea (X, Y (2), Y (3), Y (4), Leyenda (15))), Filtro de datos locales (Agregar filtro (columnas (: Fecha de finalización de la semana), Dónde (: Fecha de finalización de la semana> = 3471638400 &: Fecha de finalización de la semana <= última fecha con datos reales))), SendToReport (Despacho ({}, "Fecha de finalización de la semana", ScaleBox , {Min (3462203520), Max (última fecha con datos reales + en días (7)), Intervalo ("Año"), Inc (1), Marcas menores (0 ), Agregar línea de referencia (última fecha con datos reales, "Sólido", "Negro", Formato (última fecha con datos reales, "ddMonyyyy"), 1), Fila de etiquetas (Mostrar cuadrícula principal (1))}), Despacho ({ }, "Acumulativo [Exceso v esperado]", ScaleBox, {Mín (0), // Máx (819328.125380856), Inc (200000), Marcas menores (1), Agregar línea de referencia (0, "Sólido", "Negro" , "", 1), Fila de etiquetas (Mostrar cuadrícula principal (1))}), Despacho ({}, "Exceso v esperado", ScaleBox, {Mín (0), // Máx (26451.6956907074), Inc (10000) , Marcas menores (1), Agregar línea de referencia (0, "Sólido", "Negro", "", 1), Fila de etiquetas ({Mostrar cuadrícula principal (1), Mostrar cuadrícula secundaria (1)})}), Despacho ({}, "Todas las causas", ScaleBox, {Min (0), // Max (90470.7642581624), Inc (10000), Marcas menores (4), Fila de etiquetas ({Mostrar cuadrícula principal (1), Mostrar cuadrícula menor ( 1)})}), Despacho ({}, "400", ScaleBox, {Modelo de leyenda (22, Propiedades (0, {Color de relleno (20), Transparencia (0.5)}, ID de artículo ("Exceso v esperado", 1))), Modelo de leyenda (13, Propiedades (0, {Color de línea (21)}, ID de artículo ("Todas las causas", 1))), Modelo de leyenda (15, Prope rties (0, {Line Color (32)}, ID de artículo ("Predicción por todas las causas", 1)), Propiedades (1, {Line Color (70)}, ID de artículo ("CL superior (0,99) Todas las causas", 1)), Propiedades (2, {Color de línea (70)}, ID de artículo ("CL inferior (0,99) Todas las causas", 1)))}), Despacho ({}, "título de gráfico", TextEditBox, {Set Texto (filtro de jurisdicción || "Muertes reales frente a muertes esperadas (exceso de muertes)")}), Despacho ({}, "título del gráfico 1", TextEditBox, {Establecer texto ("Fuente: https://data.cdc.gov/NCHS/Weekly-Counts -of-Deaths-by-State-and-Select-Causes / 3yf8-kanr y https://data.cdc.gov/NCHS/Weekly-Provisional-Counts-of-Deaths-by-State-and-S/muzy -jte6 ")}), Despacho ({}," Título X ", TextEditBox, {Establecer texto (" ")}), Despacho ({}," Título Y ", TextEditBox, {Establecer texto (" Muertes en exceso acumuladas " )}), Despacho ({}, "Y 1 título", TextEditBox, {Establecer texto ("Exceso v esperado \! NDeaths por mes")}), Despacho ({}, "Y 2 título", TextEditBox, {Establecer Texto ("Muertes por mes")}), Despacho ({}, "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 Tabla de datos (); yMatrix = (dtJoined: All Cause << obtener valores) [54 :: última fila con valores reales] | / Matrix (Reverse (As List ((dtJoined: Predicted All Cause << obtener valores) [54: : última fila con datos reales]))); xMatrix = (dtJoined: Week Ending Date << obtener valores) [54 :: última fila con datos reales] | / Matrix (Reverse (As List ((dtJoined: Week Ending Date << obtener valores) [54 :: última fila con datos reales ]))); Color de relleno ("Verde oscuro medio"); Transparencia (0,5); Polígono (xMatrix, yMatrix); ), Orden de la línea de cuadrícula (1), Orden de la línea de referencia (3)}))); dtFinal << Cerrar ventana; dtProv << Cerrar ventana; dtAll << Cerrar ventana; dtSub << Cerrar ventana; dtPred << Cerrar ventana;

Analítica

Comments

Quick Links