Переписать линейную модель R на Python

Помогите переписать линейную модель R на Python. R-код:

x <- rnorm(10)
y <- 1+x+rnorm(10)
model <- lm(y~x)
res = summary(model)$r.squared
print(res)

Код Python выдает ошибку — «установка элемента массива с помощью последовательности». Кажется чего-то не хватает, не могу понять чего

x = np.random.normal(0, 1, 10)
y = [1 + np.random.normal() + v for v in x]
new_list = [x, y]
array = np.array(new_list)
df = pd.DataFrame({'x': [x], 'y': [y]})
model = LinearRegression()
X, y = df[['x', 'y']], df
model.fit(X, y)

Какой у Вас вопрос?

Martin Gal 09.04.2022 11:14

мой код не работает, выдает ошибку: «установка элемента массива с последовательностью». Я вообще не уверен, что это правильный подход

Aella 09.04.2022 11:26
Анализ настроения постов в Twitter с помощью Python, Tweepy и Flair
Анализ настроения постов в Twitter с помощью Python, Tweepy и Flair
Анализ настроения текстовых сообщений может быть настолько сложным или простым, насколько вы его сделаете. Как и в любом ML-проекте, вы можете выбрать...
7 лайфхаков для начинающих Python-программистов
7 лайфхаков для начинающих Python-программистов
В этой статье мы расскажем о хитростях и советах по Python, которые должны быть известны разработчику Python.
Установка Apache Cassandra на Mac OS
Установка Apache Cassandra на Mac OS
Это краткое руководство по установке Apache Cassandra.
Сертификатная программа "Кванты Python": Бэктестер ансамблевых методов на основе ООП
Сертификатная программа "Кванты Python": Бэктестер ансамблевых методов на основе ООП
В одном из недавних постов я рассказал о том, как я использую навыки количественных исследований, которые я совершенствую в рамках программы TPQ...
Создание персонального файлового хранилища
Создание персонального файлового хранилища
Вы когда-нибудь хотели поделиться с кем-то файлом, но он содержал конфиденциальную информацию? Многие думают, что электронная почта безопасна, но это...
Создание приборной панели для анализа данных на GCP - часть I
Создание приборной панели для анализа данных на GCP - часть I
Недавно я столкнулся с интересной бизнес-задачей - визуализацией сбоев в цепочке поставок лекарств, которую могут просматривать врачи и...
0
2
32
1
Перейти к ответу Данный вопрос помечен как решенный

Ответы 1

Ответ принят как подходящий

Библиотека statsmodels обеспечивает простую реализацию линейной регрессии со сводной таблицей, аналогичной R. Вы можете найти документацию здесь.

Питон:

import numpy as np
import statsmodels.api as sm

x = np.random.normal(0, 1, 10)
y = [1 + np.random.normal() + v for v in x]

#add intercept to x
x = sm.add_constant(x)

#statsmodels ordinary linear regression
model = sm.OLS(y, x)
results = model.fit()

print(results.summary())

Другие вопросы по теме