Логистическая регрессия

@Key27 · Регистрация: 02.04.2017

Студворк — интернет-сервис помощи студентам

Нашел готовую реализацию алгоритма логистической регрессии, но работает он довольно странно. Может быть кто-нибудь поможет найти проблему?

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
class LogisticRegression(object):
    def __init__(self, eta=0.1, n_iter=50):
        self.eta = eta
        self.n_iter = n_iter
 
    def fit(self, X, y):
        X = np.insert(X, 0, 1, axis=1)
        self.w = np.ones(X.shape[1])
        m = X.shape[0]
 
        for _ in range(self.n_iter):
            output = X.dot(self.w)
            errors = y - self._sigmoid(output)
            self.w += self.eta / m * errors.dot(X)
        return self
 
    def predict(self, X):
        output = np.insert(X, 0, 1, axis=1).dot(self.w)
        return (np.floor(self._sigmoid(output) + .5)).astype(int)
 
    def score(self, X, y):
        return sum(self.predict(X) == y) / len(y)
 
    def _sigmoid(self, x):
        return 1 / (1 + np.exp(-x))

Если делать с помощью этого класса:

Python
1
2
3
log1 = LogisticRegression().fit(X, y)
print(log1.predict(X))
print(collections.Counter(log1.predict(X)))

То результат будет такой: [0 0 0 ..., 0 0 0] Counter({0: 4899}) - одни нули.

При этом, если использовать библиотеку sklearn,

Python
1
2
3
4
5
from sklearn.linear_model import LogisticRegression
 
log2 = LogisticRegression()
log2.fit(X,y)
print(log2.predict(X))

то результат будет:
[ 1. 1. 1. ..., 0. 1. 0.] Counter({1.0: 2653, 0.0: 2246})

PS

Кликните здесь для просмотра всего текста

На всякий случай приложу кусок данных

Bash
1
2
X_list=[[10.0, 1.0, 1999.0, 14.0, 14.25, 14.0, 14.125, 9.904241, 3162400.0], [10.0, 4.0, 1999.0, 14.125, 14.125, 14.0, 14.0, 9.816596, 533400.0], [10.0, 5.0, 1999.0, 14.0, 14.125, 14.0, 14.125, 9.904241, 270600.0], [10.0, 6.0, 1999.0, 14.125, 14.875, 14.125, 14.875, 10.430133, 205400.0], [10.0, 7.0, 1999.0, 14.875, 14.875, 14.0625, 14.125, 9.904241, 214200.0], [10.0, 8.0, 1999.0, 14.1875, 14.1875, 14.0625, 14.0625, 9.860417, 65600.0], [10.0, 11.0, 1999.0, 14.0625, 14.125, 14.0, 14.0625, 9.860417, 77500.0], [10.0, 12.0, 1999.0, 14.0625, 14.125, 14.0, 14.125, 9.904241, 49200.0], [10.0, 13.0, 1999.0, 14.125, 14.375, 14.0, 14.0625, 9.860417, 49900.0], [10.0, 14.0, 1999.0, 14.0, 14.0625, 13.25, 13.625, 9.55365, 430900.0]]
y_list = [1.0, 0.0, 1.0, 1.0, 0.0, 0.0, 1.0, 1.0, 0.0, 0.0]

Далее делаю

Python
1
2
X = np.array(X_list)
y = np.array(y_list2)

и передаю X и y в функцию

@Qoffee · 17.04.2019, 02:08

Она не странно, она неверно работает. Взгляните на формулы и сверьте)

Новые блоги и статьи Все статьи Все блоги /
Доступность команды формы по условию Maks 07.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: сделать доступной кнопку (команда формы "ЗавершитьСписание") при. . .	Уведомление о неверно выбранном значении справочника Maks 06.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "НарядПутевка", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если в документе выбран неверный склад. . .	Установка Qt Creator для C и C++: ставим среду, CMake и MinGW без фреймворка Qt 8Observer8 05.04.2026 Среду разработки Qt Creator можно установить без фреймворка Qt. Есть отдельный репозиторий для этой среды: https:/ / github. com/ qt-creator/ qt-creator, где можно скачать установщик, на вкладке Releases:. . .	AkelPad-скрипты, структуры, и немного лирики.. testuser2 05.04.2026 Такая программа, как AkelPad существует уже давно, и также давно существуют скрипты под нее. Тем не менее, прога живет, периодически что-то не спеша дополняется, улучшается. Что меня в первую очередь. . .
Отображение реквизитов в документе по условию и контроль их заполнения Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеСпецтехники", разработанного в конфигурации КА2. Данный документ берёт данные из другого нетипового документа. . .	Фото всей Земли с борта корабля Orion миссии Artemis II kumehtar 04.04.2026 Это первое подобное фото сделанное человеком за 50 лет. Снимок называют новым вариантом легендарной фотографии «The Blue Marble» 1972 года, сделанной с борта корабля «Аполлон-17». Новое фото. . .	Вывод диалогового окна перед закрытием, если документ не проведён Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать программный контроль на предмет проведения документа. . .	Программный контроль заполнения реквизитов табличной части документа Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: 1. Реализовать контроль заполнения реквизита. . .

@Qoffee 0 / 0 / 0 Регистрация: 20.09.2016 Сообщений: 3
	17.04.2019, 02:08
	Она не странно, она неверно работает. Взгляните на формулы и сверьте) 0