Articles of statsmodels

Prueba de cointegración de Johansen en python

No puedo encontrar ninguna referencia sobre la funcionalidad para realizar la prueba de cointegración de Johansen en ningún módulo de Python que se ocupe de las estadísticas y el análisis de series de tiempo (pandas y statsmodel). ¿Alguien sabe si hay algún código alrededor que pueda realizar tal prueba de cointegración entre series de tiempo? […]

¿Cómo funciona exactamente BIC en pruebas de Dickey aumentada-Fuller en Python?

Esta pregunta es sobre la implementación de la prueba de Dickey-Fuller aumentada en statsmodels.tsa.stattools python library – adfuller (). En principio, se supone que AIC y BIC computan el criterio de información para un conjunto de modelos disponibles y seleccionan el mejor (el que tiene la pérdida de información más baja). Pero, ¿cómo operan en […]

¿Cómo calcular la probabilidad de ajuste de la curva en scipy?

Tengo un ajuste de modelo no lineal que se ve así: La línea sólida oscura es el ajuste del modelo, y la parte gris son los datos en bruto. Versión corta de la pregunta: ¿cómo obtengo la probabilidad de que este modelo se ajuste, por lo que puedo realizar la prueba de relación de probabilidad […]

¿Cómo generar muchos términos de interacción en Pandas?

Me gustaría estimar un modelo de regresión IV usando muchas interacciones con maniquíes de año, demográficos, etc. No puedo encontrar un método explícito para hacer esto en Pandas y tengo curiosidad si alguien tiene consejos. Estoy pensando en probar scikit-learn y esta función: http://scikit-learn.org/stable/modules/generated/sklearn.preprocessing.PolynomialFeatures.html

Intervalo de confianza para LOWESS en Python

¿Cómo podría calcular los intervalos de confianza para una regresión LOWESS en Python? Me gustaría agregar estos como una región sombreada a la gráfica LOESS creada con el siguiente código (otros paquetes además de statsmodels también están bien). import numpy as np import pylab as plt import statsmodels.api as sm x = np.linspace(0,2*np.pi,100) y = […]

Cómo obtener la intersección de regresión usando Statsmodels.api

Estoy intentando calcular una salida de regresión utilizando la biblioteca de Python, pero no puedo obtener el valor de intercepción cuando uso la biblioteca: import statsmodels.api as sm Imprime todo el análisis de regresión excepto el intercepto. pero cuando uso: from pandas.stats.api import ols Mi código para pandas: Regression = ols(y= Sorted_Data3[‘net_realization_rate’],x = Sorted_Data3[[‘Cohort_2′,’Cohort_3’]]) print […]

Modelo ARIMA de Python, los valores predichos se desplazan

Soy nuevo en la implementación de Python ARIMA. Tengo unos datos a 15 min de frecuencia durante unos meses. En mi bash de seguir el método de Box-Jenkins para ajustar un modelo de series temporales. Me encontré con un problema hacia el final. Se dan el gráfico ACF-PACF para la serie de tiempo (ts) y […]

Deshacer una serie de diferencias

Tengo una serie de pandas con datos mensuales ( df.sales ). Necesitaba restar los datos 12 meses antes para ajustar una serie de tiempo, así que ejecuté este comando: sales_new = df.sales.diff(periods=12) Entonces me ajusté a un modelo ARMA y predije el futuro: model = ARMA(sales_new, order=(2,0)).fit() model.predict(‘2015-01-01’, ‘2017-01-01’) Debido a que había difuminado los […]

Cómo configurar correctamente los parámetros de inicio / finalización de statsmodels.tsa.ar_model.AR.predict

Tengo una estructura de datos de los costos del proyecto de una serie de tiempo espaciada irregularmente en la que me gustaría intentar aplicar el modelo de statsmodel AR . Esta es una muestra de los datos en su dataframe: cost date 2015-07-16 35.98 2015-08-11 25.00 2015-08-11 43.94 2015-08-13 26.25 2015-08-18 15.38 2015-08-24 77.72 2015-09-09 […]

Statsmodels Poisson glm diferente de R

Estoy tratando de ajustar algunos modelos (modelos de interacción espacial) de acuerdo con el código que se proporciona en R. He podido hacer que algunos de los códigos funcionen utilizando statsmodels en un marco de python, pero algunos de ellos no coinciden en absoluto. Creo que el código que tengo para R y Python debería […]