Нахождение среднеквадратичной ошибки для линейной регрессии в python (с помощью scikit learn) - PullRequest
0 голосов
/ 27 июня 2018

Я пытаюсь сделать простую линейную регрессию в Python с переменной x, являющейся словом количество описания проекта и значение y, являющееся скоростью финансирования в днях.

Я немного запутался, поскольку среднеквадратическая ошибка (RMSE) для теста составляет 13,77. и 13,88 для данных обучения. Во-первых, не должно ли RMSE быть между 0 и 1? И, во-вторых, разве RMSE для тестовых данных должна быть выше, чем для обучающих данных? Так что, я думаю, я сделал что-то не так, но не уверен, где ошибка.

Кроме того, мне нужно знать весовой коэффициент для регрессии, но, к сожалению, не знаю, как его распечатать, так как он скрыт внутри методов sklearn. Кто-нибудь может помочь здесь?

Это то, что я имею до сих пор:

import numpy as np
import matplotlib.pyplot as plt
import sqlite3
from sklearn.model_selection import train_test_split
from sklearn import linear_model

con = sqlite3.connect('database.db')
cur = con.cursor()

# y-variable in regression is funding speed ("DAYS_NEEDED")    
cur.execute("SELECT DAYS_NEEDED FROM success")
y = cur.fetchall()                  # list of tuples
y = np.array([i[0] for i in y])     # list of int   # y.shape = (1324476,)

# x-variable in regression is the project description length ("WORD_COUNT")
cur.execute("SELECT WORD_COUNT FROM success")
x = cur.fetchall()
x = np.array([i[0] for i in x])     # list of int   # x.shape = (1324476,)

# Get the train and test data split
x_train, x_test, y_train, y_test = train_test_split(x, y, test_size=0.2, random_state=42)

# Fit a model
lm = linear_model.LinearRegression()
x_train = x_train.reshape(-1, 1)    # new shape: (1059580, 1)
y_train = y_train.reshape(-1, 1)    # new shape: (1059580, 1)
model = lm.fit(x_train, y_train)
x_test = x_test.reshape(-1, 1)      # new shape: (264896, 1)
predictions_test = lm.predict(x_test)
predictions_train = lm.predict(x_train)

print("y_test[5]: ", y_test[5])     # 14
print("predictions[5]: ", predictions_test[5]) # [ 12.6254537]

# Calculate the root mean square error (RMSE) for test and training data
N = len(y_test)
rmse_test = np.sqrt(np.sum((np.array(y_test).flatten() - np.array(predictions_test).flatten())**2)/N)
print("RMSE TEST: ", rmse_test)     # 13.770731326

N = len(y_train)
rmse_train = np.sqrt(np.sum((np.array(y_train).flatten() - np.array(predictions_train).flatten())**2)/N)
print("RMSE train: ", rmse_train)   # 13.8817814595

Любая помощь очень ценится! Спасибо!

1 Ответ

0 голосов
/ 27 июня 2018
  1. RMSE имеет ту же единицу, что и зависимая переменная. Это означает, что если переменная, которую вы пытаетесь предсказать, варьируется от 0 до 100, среднеквадратическое среднеквадратичное значение 99 ужасно! Если, скажем, у вас RMSE 5 для данных в диапазоне от 0 до 100, RMSE 5 впечатляет. НО, если RMSE равно 5 для данных в диапазоне от 1 до 10, то у вас проблема! Я надеюсь, что это способно привести точку в чувство.

  2. Поскольку RMSE вашего поезда и теста похожи, похлопайте себя по спине! Вы на самом деле сделали хорошую работу! Если RMSE теста> тренируется, вы слегка одеваетесь.

В соответствии с тем, что Уманг сказал в комментариях, вы используете model.coef_ и model.intercept_, чтобы напечатать веса, которые ваша модель рассчитала, чтобы быть оптимальными.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...