Articles of statsmodels

Pandas con efectos fijos

Estoy usando Pandas en Python 2.7. Tengo datos con las siguientes columnas: estado, año, tasa de desempate, salario Estoy enseñando un curso sobre cómo usar Python para la investigación. Como la culminación de nuestro proyecto, quiero realizar una regresión de UnempRate en el control de salarios para efectos fijos estatales y anuales. Puedo hacer esto […]

Predecir los valores futuros mediante la regresión OLS (Python, StatsModels, Pandas)

Actualmente estoy intentando implementar un MLR en Python y no estoy seguro de cómo aplico los coeficientes que he encontrado a los valores futuros. import pandas as pd import statsmodels.formula.api as sm import statsmodels.api as sm2 TV = [230.1, 44.5, 17.2, 151.5, 180.8] Radio = [37.8,39.3,45.9,41.3,10.8] Newspaper = [69.2,45.1,69.3,58.5,58.4] Sales = [22.1, 10.4, 9.3, 18.5,12.9] […]

Python pandas no tiene atributos ols – Error (Rolling OLS)

Para mi evaluación, quería ejecutar una OLS regression estimation de OLS regression estimation 1000 ventanas rodantes del conjunto de datos que se encuentra en esta URL: https://drive.google.com/open?id=0B2Iv8dfU4fTUa3dPYW5tejA0bzg usando el siguiente script Python . # /usr/bin/python -tt import numpy as np import matplotlib.pyplot as plt import pandas as pd from statsmodels.formula.api import ols df = pd.read_csv(‘estimated.csv’, […]

Errores estándar inesperados con mínimos cuadrados ponderados en Python Pandas

En el código de la clase principal de OLS en Python Pandas , busco ayuda para aclarar qué convenciones se usan para el error estándar y las estadísticas t que se informan cuando se realiza el OLS ponderado. Aquí está mi ejemplo de conjunto de datos, con algunas importaciones para usar Pandas y usar scikits.statsmodels […]

¿Cómo tener varios grupos en el modelo de efectos mixtos lineales de Python statsmodels?

Estoy tratando de usar el modelo de efectos mixtos lineales de Python statsmodels para ajustar un modelo que tiene dos intercepciones aleatorias, por ejemplo, dos grupos. No puedo averiguar cómo inicializar el modelo para poder hacer esto. Aquí está el ejemplo. Tengo datos que se parecen a los siguientes (tomados de aquí ): subject gender […]

¿Qué usar para hacer correlación múltiple?

Estoy tratando de usar python para calcular la regresión lineal múltiple y la correlación múltiple entre una matriz de respuesta y un conjunto de matrices de predictores. Vi el ejemplo muy simple para calcular la regresión lineal múltiple, que es fácil. Pero, ¿cómo calcular la correlación múltiple con statsmodels? O con cualquier otra cosa, como […]

¿Cómo calcular el intervalo de confianza del 99% para la pendiente en un modelo de regresión lineal en python?

Tenemos la siguiente regresión lineal: y ~ b0 + b1 * x1 + b2 * x2. Sé que la función de regresión en Matlab sí la calcula, pero linalg.lstsq de numpy no ( https://docs.scipy.org/doc/numpy-dev/user/numpy-for-matlab-users.html ).

Python ARIMA variable exógena fuera de muestra

Estoy tratando de predecir una serie de tiempo en el paquete ARIMA de statsmodels de python con la inclusión de una variable exógena, pero no puedo encontrar la manera correcta de insertar la variable exógena en el paso de predicción. Vea aquí para los documentos. import numpy as np from scipy import stats import pandas […]

Python Statsmodel ARIMA start

Acabo de empezar a trabajar en el análisis de series de tiempo usando statsmodels. Tengo un conjunto de datos con fechas y valores (durante aproximadamente 3 meses). Estoy enfrentando algunos problemas para proporcionar el orden correcto al modelo ARIMA. Estoy buscando ajustar las tendencias y la estacionalidad y luego calcular los valores atípicos. Mis ‘valores’ […]

Pandas / Statsmodel OLS prediciendo valores futuros

He estado tratando de obtener una predicción de valores futuros en un modelo que he creado. He probado tanto OLS en pandas como statsmodels. Esto es lo que tengo en statsmodels: import statsmodels.api as sm endog = pd.DataFrame(dframe[‘monthly_data_smoothed8’]) smresults = sm.OLS(dframe[‘monthly_data_smoothed8’], dframe[‘date_delta’]).fit() sm_pred = smresults.predict(endog) sm_pred La longitud de la matriz devuelta es igual a […]