Numpy에서 다항식 회귀 통계 구하기

Question

Numpy에서 다항식 회귀 (polyfit)를 사용하여 계수와 요격을 얻으 려하지만 다항식 함수를 얻기 위해 스크립트를 작성하는 방법을 모르겠습니다. 나는 아래 첨부 내가 코드 선형 회귀에 대한 코드를 만든

는 :

import matplotlib.pyplot as plt 
import sys 
from numpy import * 
import numpy as np 
import numpy.polynomial.polynomial as poly 
import pylab 
from scipy import stats 
from scipy.interpolate import * 
from datetime import datetime, timedelta 

#Open dataset1,dataset2 data 

data1 = np.loadtxt('/home/script/2_columns', delimiter=',', skiprows=0) 

data2 = np.loadtxt('/home/script/2_columns_a', delimiter=',', skiprows=0) 

#Define first column as dataset1 

#Define second column as dataset2 

x = data1[:,0] 
y = data1[:,1] 

#The stuff... 
slope, intercept, r_value, p_value, std_err = stats.linregress(x,y) 

Answer 1

당신은 통계가 필요하십니까? 방금 계수를 필요로하지 않을 경우, 그것은 numpy.polyfit을 사용하여 정말 아주 간단합니다 : 당신이 다른 값을 필요로 할 경우

import numpy as np 

# from your code 
data1 = np.loadtxt('/home/script/2_columns', delimiter=',', skiprows=0) 

x = data1[:,0] 
y = data1[:,1] 

degree = 3 
coeffs = np.polyfit(x, y, degree) 

# now, coeffs is an array which contains the polynomial coefficients 
# in ascending order, i.e. x^0, x^1, x^2 
intercept, linear, quadratic, cubic = coeffs 

는, 예를 들어 r_value이 x 사이의 상관 계수이기 때문에 당신이 필요한 지정하세요 y으로, 데이터가 선형이 아닌 것으로 예상 할 때 그다지 유용하지 않습니다.