Python R 线性回归 高斯回归 比较

使用的数据是公路一氧化碳数据，相应细节可参见下面链接：
数据下载链接：http://www.statsci.org/data/general/cofreewy.html

R设定工作目录指令setwd
下面先使用R 的逐步回归选取AIC最小的普通线性模型实行最小二乘估计：

w = read.table("COfreewy.txt", header = T)
a = lm(CO~.,w)
summary(a)
b = step(a, direction = 'backward')
summary(b)
shapiro.test(b$res)

并后面使用斯皮尔曼正态检验，不能拒绝正态性假设。（0.05601）

下面给出上述过程的python简化版本：

from __future__ import division

import numpy as np 
import scipy as sci
import pandas as pd
from math import log
from sklearn import linear_model 
from sklearn.gaussian_process import GaussianProcess

from sklearn.cross_validation import LeaveOneOut

import itertools