2011-11-10から1日間の記事一覧
#-*- coding:utf-8 -*- import numpy as np from numpy import linalg as LA def linreg(input,target): nData = np.shape(input)[0] input = np.concatenate((-np.ones((nData,1)),input),axis=1) a = LA.inv(np.dot(input.T,input)) b = np.dot(a,input.T)…