# 计算传播与机器学习

wangchengjun@nju.edu.cn

## 1、 监督式学习

• 这个算法由一个目标变量或结果变量（或因变量）组成。
• 这些变量由已知的一系列预示变量（自变量）预测而来。
• 利用这一系列变量，我们生成一个将输入值映射到期望输出值的函数。
• 这个训练过程会一直持续，直到模型在训练数据上获得期望的精确度。
• 监督式学习的例子有：回归、决策树、随机森林、K – 近邻算法、逻辑回归等。

## 2、非监督式学习

• 在这个算法中，没有任何目标变量或结果变量要预测或估计。
• 这个算法用在不同的组内聚类分析。
• 这种分析方式被广泛地用来细分客户，根据干预的方式分为不同的用户组。
• 非监督式学习的例子有：关联算法和 K–均值算法。

## 3、强化学习

• 这个算法训练机器进行决策。
• 它是这样工作的：机器被放在一个能让它通过反复试错来训练自己的环境中。
• 机器从过去的经验中进行学习，并且尝试利用了解最透彻的知识作出精确的商业判断。
• 强化学习的例子有马尔可夫决策过程。alphago

Chess. Here, the agent decides upon a series of moves depending on the state of the board (the environment), and the reward can be defined as win or lose at the end of the game:

• 线性回归
• 逻辑回归
• 决策树
• SVM
• 朴素贝叶斯

• K最近邻算法
• K均值算法
• 随机森林算法
• 降维算法

# 使用sklearn做线性回归

wangchengjun@nju.edu.cn

# 线性回归

• 通常用于估计连续性变量的实际数值（房价、呼叫次数、总销售额等）。
• 通过拟合最佳直线来建立自变量X和因变量Y的关系。
• 这条最佳直线叫做回归线，并且用 $Y= \beta *X + C$ 这条线性等式来表示。
• 系数 a 和 b 可以通过最小二乘法获得


In :

%matplotlib inline

from sklearn import datasets
from sklearn import linear_model
import matplotlib.pyplot as plt
from sklearn.metrics import classification_report
from sklearn.preprocessing import scale

import sklearn
print sklearn.__version__




0.17.1




In :

# boston data
y = boston.target
X = boston.data




In :

' '.join(dir(boston))




Out:

'__class__ __cmp__ __contains__ __delattr__ __delitem__ __dict__ __doc__ __eq__ __format__ __ge__ __getattr__ __getattribute__ __getitem__ __gt__ __hash__ __init__ __iter__ __le__ __len__ __lt__ __module__ __ne__ __new__ __reduce__ __reduce_ex__ __repr__ __setattr__ __setitem__ __setstate__ __sizeof__ __str__ __subclasshook__ __weakref__ clear copy fromkeys get has_key items iteritems iterkeys itervalues keys pop popitem setdefault update values viewitems viewkeys viewvalues'




In :

boston['feature_names']




Out:

array(['CRIM', 'ZN', 'INDUS', 'CHAS', 'NOX', 'RM', 'AGE', 'DIS', 'RAD',
'TAX', 'PTRATIO', 'B', 'LSTAT'],
dtype='|S7')




In :

import numpy as np
import statsmodels.api as sm
import statsmodels.formula.api as smf

# Fit regression model (using the natural log of one of the regressors)
results = smf.ols('boston.target ~ boston.data', data=boston).fit()

print results.summary()




OLS Regression Results
==============================================================================
Dep. Variable:          boston.target   R-squared:                       0.741
Method:                 Least Squares   F-statistic:                     108.1
Date:                Sun, 08 May 2016   Prob (F-statistic):          6.95e-135
Time:                        20:14:27   Log-Likelihood:                -1498.8
No. Observations:                 506   AIC:                             3026.
Df Residuals:                     492   BIC:                             3085.
Df Model:                          13
Covariance Type:            nonrobust
===================================================================================
coef    std err          t      P>|t|      [95.0% Conf. Int.]
-----------------------------------------------------------------------------------
Intercept          36.4911      5.104      7.149      0.000        26.462    46.520
boston.data     -0.1072      0.033     -3.276      0.001        -0.171    -0.043
boston.data      0.0464      0.014      3.380      0.001         0.019     0.073
boston.data      0.0209      0.061      0.339      0.735        -0.100     0.142
boston.data      2.6886      0.862      3.120      0.002         0.996     4.381
boston.data    -17.7958      3.821     -4.658      0.000       -25.302   -10.289
boston.data      3.8048      0.418      9.102      0.000         2.983     4.626
boston.data      0.0008      0.013      0.057      0.955        -0.025     0.027
boston.data     -1.4758      0.199     -7.398      0.000        -1.868    -1.084
boston.data      0.3057      0.066      4.608      0.000         0.175     0.436
boston.data     -0.0123      0.004     -3.278      0.001        -0.020    -0.005
boston.data    -0.9535      0.131     -7.287      0.000        -1.211    -0.696
boston.data     0.0094      0.003      3.500      0.001         0.004     0.015
boston.data    -0.5255      0.051    -10.366      0.000        -0.625    -0.426
==============================================================================
Omnibus:                      178.029   Durbin-Watson:                   1.078
Prob(Omnibus):                  0.000   Jarque-Bera (JB):              782.015
Skew:                           1.521   Prob(JB):                    1.54e-170
Kurtosis:                       8.276   Cond. No.                     1.51e+04
==============================================================================

Warnings:
 Standard Errors assume that the covariance matrix of the errors is correctly specified.
 The condition number is large, 1.51e+04. This might indicate that there are
strong multicollinearity or other numerical problems.




In :

regr = linear_model.LinearRegression()
lm = regr.fit(boston.data, y)




In :

lm.intercept_, lm.coef_, lm.score(boston.data, y)




Out:

(36.491103280361614,
array([ -1.07170557e-01,   4.63952195e-02,   2.08602395e-02,
2.68856140e+00,  -1.77957587e+01,   3.80475246e+00,
7.51061703e-04,  -1.47575880e+00,   3.05655038e-01,
-1.23293463e-02,  -9.53463555e-01,   9.39251272e-03,
-5.25466633e-01]),
0.7406077428649428)




In :

predicted = regr.predict(boston.data)




In :

fig, ax = plt.subplots()
ax.scatter(y, predicted)
ax.plot([y.min(), y.max()], [y.min(), y.max()], 'k--', lw=4)
ax.set_xlabel('$Measured$', fontsize = 20)
ax.set_ylabel('$Predicted$', fontsize = 20)
plt.show()






## 训练集和测试集



In :

boston.data




Out:

array([[  6.32000000e-03,   1.80000000e+01,   2.31000000e+00, ...,
1.53000000e+01,   3.96900000e+02,   4.98000000e+00],
[  2.73100000e-02,   0.00000000e+00,   7.07000000e+00, ...,
1.78000000e+01,   3.96900000e+02,   9.14000000e+00],
[  2.72900000e-02,   0.00000000e+00,   7.07000000e+00, ...,
1.78000000e+01,   3.92830000e+02,   4.03000000e+00],
...,
[  6.07600000e-02,   0.00000000e+00,   1.19300000e+01, ...,
2.10000000e+01,   3.96900000e+02,   5.64000000e+00],
[  1.09590000e-01,   0.00000000e+00,   1.19300000e+01, ...,
2.10000000e+01,   3.93450000e+02,   6.48000000e+00],
[  4.74100000e-02,   0.00000000e+00,   1.19300000e+01, ...,
2.10000000e+01,   3.96900000e+02,   7.88000000e+00]])




In :

from sklearn.cross_validation import train_test_split
Xs_train, Xs_test, y_train, y_test = train_test_split(boston.data, boston.target, test_size=0.2, random_state=42)




In :

regr = linear_model.LinearRegression()
lm = regr.fit(Xs_train, y_train)




In :

lm.intercept_, lm.coef_, lm.score(Xs_train, y_train)




Out:

(30.288948339368776,
array([ -1.12463481e-01,   3.00810168e-02,   4.07309919e-02,
2.78676719e+00,  -1.72406347e+01,   4.43248784e+00,
-6.23998173e-03,  -1.44848504e+00,   2.62113793e-01,
-1.06390978e-02,  -9.16398679e-01,   1.24516469e-02,
-5.09349120e-01]),
0.75088377867329148)




In :

predicted = regr.predict(Xs_test)




In :

fig, ax = plt.subplots()
ax.scatter(y_test, predicted)
ax.plot([y.min(), y.max()], [y.min(), y.max()], 'k--', lw=4)
ax.set_xlabel('$Measured$', fontsize = 20)
ax.set_ylabel('$Predicted$', fontsize = 20)
plt.show()






# 交叉验证



In :

from sklearn.cross_validation import cross_val_score

regr = linear_model.LinearRegression()
scores = cross_val_score(regr, boston.data , boston.target, cv = 3)
scores.mean()




Out:

-1.5787701857179994




In :

scores = [cross_val_score(regr, data_X_scale, boston.target, cv = int(i)).mean() for i in range(3, 50)]
plt.plot(range(3, 50), scores,'r-o')
plt.show()







In :

data_X_scale = scale(boston.data)
scores = cross_val_score(regr, boston.data, boston.target, cv = 7)
scores.mean()




Out:

0.453848713596956



# 使用天涯bbs数据



In :

import pandas as pd

df[:2]




Out:

title
author
author_page
click
time

0
【民间语文第161期】宁波px启示:船进港湾人应上岸
/post-free-2849477-1.shtml
贾也
http://www.tianya.cn/50499450
194675
2703
2012-10-29 07:59

1
宁波镇海PX项目引发群体上访 当地政府发布说明(转载)
/post-free-2839539-1.shtml
无上卫士ABC
http://www.tianya.cn/74341835
88244
1041
2012-10-24 12:41




In :

def randomSplit(dataX, dataY, num):
dataX_train = []
dataX_test = []
dataY_train = []
dataY_test = []
import random
test_index = random.sample(range(len(df)), num)
for k in range(len(dataX)):
if k in test_index:
dataX_test.append([dataX[k]])
dataY_test.append(dataY[k])
else:
dataX_train.append([dataX[k]])
dataY_train.append(dataY[k])
return dataX_train, dataX_test, dataY_train, dataY_test,




In :

import numpy as np

# Use only one feature
# Split the data into training/testing sets
data_X_train, data_X_test, data_y_train, data_y_test = randomSplit(np.log(df.click+1), np.log(df.reply+1), 20)
# Create linear regression object
regr = linear_model.LinearRegression()
# Train the model using the training sets
regr.fit(data_X_train, data_y_train)
# Explained variance score: 1 is perfect prediction
print'Variance score: %.2f' % regr.score(data_X_test, data_y_test)




Variance score: 0.78




In :

y_true, y_pred = data_y_test, regr.predict(data_X_test)




In :

plt.scatter(y_pred, y_true,  color='black')
plt.show()







In :

# Plot outputs
plt.scatter(data_X_test, data_y_test,  color='black')
plt.plot(data_X_test, regr.predict(data_X_test), color='blue', linewidth=3)
plt.show()







In :

# The coefficients
print 'Coefficients: \n', regr.coef_




Coefficients:
[ 0.68250313]




In :

# The mean square error
print "Residual sum of squares: %.2f" % np.mean((regr.predict(data_X_test) - data_y_test) ** 2)




Residual sum of squares: 0.42




In :

df.click_log = [[df.click[i]] for i in range(len(df))]




In :

from sklearn.cross_validation import train_test_split
Xs_train, Xs_test, y_train, y_test = train_test_split(df.click_log, df.reply_log,test_size=0.2, random_state=0)

# Create linear regression object
regr = linear_model.LinearRegression()
# Train the model using the training sets
regr.fit(Xs_train, y_train)
# Explained variance score: 1 is perfect prediction
print'Variance score: %.2f' % regr.score(Xs_test, y_test)




Variance score: 0.95




In :

# Plot outputs
plt.scatter(Xs_test, y_test,  color='black')
plt.plot(Xs_test, regr.predict(Xs_test), color='blue', linewidth=3)
plt.show()







In :

from sklearn.cross_validation import cross_val_score

regr = linear_model.LinearRegression()
scores = cross_val_score(regr, df.click_log, df.reply_log, cv = 3)
scores.mean()




Out:

0.21630869764168115




In :

regr = linear_model.LinearRegression()
scores = cross_val_score(regr, df.click_log, df.reply_log, cv = 4)
scores.mean()




Out:

0.17723080221577134



# 使用sklearn做logistic回归

wangchengjun@nju.edu.cn

• logistic回归是一个分类算法而不是一个回归算法。
• 可根据已知的一系列因变量估计离散数值（比方说二进制数值 0 或 1 ，是或否，真或假）。
• 简单来说，它通过将数据拟合进一个逻辑函数(logistic function)来预估一个事件出现的概率。
• 因此，它也被叫做逻辑回归。因为它预估的是概率，所以它的输出值大小在 0 和 1 之间（正如所预计的一样）。
$$odds= \frac{p}{1-p} = \frac{probability\: of\: event\: occurrence} {probability \:of \:not\: event\: occurrence}$$
$$ln(odds)= ln(\frac{p}{1-p})$$
$$logit(x) = ln(\frac{p}{1-p}) = b_0+b_1X_1+b_2X_2+b_3X_3....+b_kX_k$$



In :

repost = []
for i in df.title:
if u'转载' in i.decode('utf8'):
repost.append(1)
else:
repost.append(0)




In :

data_X = [[df.click[i], df.reply[i]] for i in range(len(df))]
data_X[:3]




Out:

[[194675, 2703], [88244, 1041], [82779, 625]]




In :

from sklearn.linear_model import LogisticRegression
df['repost'] = repost
model = LogisticRegression()
model.fit(data_X,df.repost)
model.score(data_X,df.repost)




Out:

0.61241970021413272




In :

def randomSplitLogistic(dataX, dataY, num):
dataX_train = []
dataX_test = []
dataY_train = []
dataY_test = []
import random
test_index = random.sample(range(len(df)), num)
for k in range(len(dataX)):
if k in test_index:
dataX_test.append(dataX[k])
dataY_test.append(dataY[k])
else:
dataX_train.append(dataX[k])
dataY_train.append(dataY[k])
return dataX_train, dataX_test, dataY_train, dataY_test,




In :

# Split the data into training/testing sets
data_X_train, data_X_test, data_y_train, data_y_test = randomSplitLogistic(data_X, df.repost, 20)
# Create logistic regression object
log_regr = LogisticRegression()
# Train the model using the training sets
log_regr.fit(data_X_train, data_y_train)
# Explained variance score: 1 is perfect prediction
print'Variance score: %.2f' % log_regr.score(data_X_test, data_y_test)




Variance score: 0.40




In :

y_true, y_pred = data_y_test, log_regr.predict(data_X_test)




In :

y_true, y_pred




Out:

([0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1],
array([1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]))




In :

print(classification_report(y_true, y_pred))




precision    recall  f1-score   support

0       0.00      0.00      0.00        10
1       0.50      1.00      0.67        10

avg / total       0.25      0.50      0.33        20




In :

from sklearn.cross_validation import train_test_split
Xs_train, Xs_test, y_train, y_test = train_test_split(data_X, df.repost, test_size=0.2, random_state=42)




In :

# Create logistic regression object
log_regr = LogisticRegression()
# Train the model using the training sets
log_regr.fit(Xs_train, y_train)
# Explained variance score: 1 is perfect prediction
print'Variance score: %.2f' % log_regr.score(Xs_test, y_test)




Variance score: 0.60




In :

print('Logistic score for test set: %f' % log_regr.score(Xs_test, y_test))
print('Logistic score for training set: %f' % log_regr.score(Xs_train, y_train))
y_true, y_pred = y_test, log_regr.predict(Xs_test)
print(classification_report(y_true, y_pred))




Logistic score for test set: 0.595745
Logistic score for training set: 0.613941
precision    recall  f1-score   support

0       1.00      0.03      0.05        39
1       0.59      1.00      0.74        55

avg / total       0.76      0.60      0.46        94




In :

logre = LogisticRegression()
scores = cross_val_score(logre, data_X, df.repost, cv = 3)
scores.mean()




Out:

0.53333333333333333




In :

logre = LogisticRegression()
data_X_scale = scale(data_X)
# The importance of preprocessing in data science and the machine learning pipeline I:
scores = cross_val_score(logre, data_X_scale, df.repost, cv = 3)
scores.mean()




Out:

0.62948717948717947



# 使用sklearn实现贝叶斯预测

wangchengjun@nju.edu.cn

# Naive Bayes algorithm

It is a classification technique based on Bayes’ Theorem with an assumption of independence among predictors.

In simple terms, a Naive Bayes classifier assumes that the presence of a particular feature in a class is unrelated to the presence of any other feature.

why it is known as ‘Naive’? For example, a fruit may be considered to be an apple if it is red, round, and about 3 inches in diameter. Even if these features depend on each other or upon the existence of the other features, all of these properties independently contribute to the probability that this fruit is an apple.

$$p(c|x) = \frac{p(x|c) p(c)}{p(x)}$$
• P(c|x) is the posterior probability of class (c, target) given predictor (x, attributes).
• P(c) is the prior probability of class.
• P(x|c) is the likelihood which is the probability of predictor given class.
• P(x) is the prior probability of predictor.

Step 1: Convert the data set into a frequency table

Step 2: Create Likelihood table by finding the probabilities like:

• p(Overcast) = 0.29, p(rainy) = 0.36, p(sunny) = 0.36
• p(playing) = 0.64, p(rest) = 0.36

Step 3: Now, use Naive Bayesian equation to calculate the posterior probability for each class. The class with the highest posterior probability is the outcome of prediction.

## Problem: Players will play if weather is sunny. Is this statement is correct?

We can solve it using above discussed method of posterior probability.

$P(Yes | Sunny) = \frac{P( Sunny | Yes) * P(Yes) } {P (Sunny)}$

Here we have P (Sunny |Yes) = 3/9 = 0.33, P(Sunny) = 5/14 = 0.36, P( Yes)= 9/14 = 0.64

Now, $P (Yes | Sunny) = \frac{0.33 * 0.64}{0.36} = 0.60$, which has higher probability.



In :

from sklearn import naive_bayes
'  '.join(dir(naive_bayes))




Out:

'ABCMeta  BaseDiscreteNB  BaseEstimator  BaseNB  BernoulliNB  ClassifierMixin  GaussianNB  LabelBinarizer  MultinomialNB  __all__  __builtins__  __doc__  __file__  __name__  __package__  _check_partial_fit_first_call  abstractmethod  binarize  check_X_y  check_array  check_is_fitted  in1d  issparse  label_binarize  logsumexp  np  safe_sparse_dot  six'


• naive_bayes.GaussianNB Gaussian Naive Bayes (GaussianNB)
• naive_bayes.MultinomialNB([alpha, ...]) Naive Bayes classifier for multinomial models
• naive_bayes.BernoulliNB([alpha, binarize, ...]) Naive Bayes classifier for multivariate Bernoulli models.


In :

#Import Library of Gaussian Naive Bayes model
from sklearn.naive_bayes import GaussianNB
import numpy as np

#assigning predictor and target variables
x= np.array([[-3,7],[1,5], [1,2], [-2,0], [2,3], [-4,0], [-1,1], [1,1], [-2,2], [2,7], [-4,1], [-2,7]])
Y = np.array([3, 3, 3, 3, 4, 3, 3, 4, 3, 4, 4, 4])




In :

#Create a Gaussian Classifier
model = GaussianNB()

# Train the model using the training sets
model.fit(x[:8], Y[:8])

#Predict Output
predicted= model.predict([[1,2],[3,4]])
print predicted




[4 3]




In :

model.score(x[8:], Y[8:])




Out:

0.25



# cross-validation

k-fold CV, the training set is split into k smaller sets (other approaches are described below, but generally follow the same principles). The following procedure is followed for each of the k “folds”:

• A model is trained using k-1 of the folds as training data;
• the resulting model is validated on the remaining part of the data (i.e., it is used as a test set to compute a performance measure such as accuracy).


In :

data_X_train, data_X_test, data_y_train, data_y_test = randomSplit(df.click, df.reply, 20)
# Train the model using the training sets
model.fit(data_X_train, data_y_train)

#Predict Output
predicted= model.predict(data_X_test)
print predicted




[41 16  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0]




In :

model.score(data_X_test, data_y_test)




Out:

0.65000000000000002




In :

from sklearn.cross_validation import cross_val_score

model = GaussianNB()
scores = cross_val_score(model, [[c] for c in df.click], df.reply, cv = 5)
scores.mean()




Out:

0.49403904714780522



# 使用sklearn实现决策树

wangchengjun@nju.edu.cn

# 决策树

• 这个监督式学习算法通常被用于分类问题。
• 它同时适用于分类变量和连续因变量。
• 在这个算法中，我们将总体分成两个或更多的同类群。
• 这是根据最重要的属性或者自变量来分成尽可能不同的组别。

## 在上图中你可以看到，根据多种属性，人群被分成了不同的四个小组，来判断 “他们会不会去玩”。

### 为了把总体分成不同组别，需要用到许多技术，比如说 Gini、Information Gain、Chi-square、entropy。



In :

from sklearn import tree
model = tree.DecisionTreeClassifier(criterion='gini')




In :

data_X_train, data_X_test, data_y_train, data_y_test = randomSplitLogistic(data_X, df.repost, 20)
model.fit(data_X_train,data_y_train)
model.score(data_X_train,data_y_train)




Out:

0.91722595078299773




In :

# Predict
model.predict(data_X_test)




Out:

array([1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1])




In :

# crossvalidation
scores = cross_val_score(model, data_X, df.repost, cv = 3)
scores.mean()




Out:

0.33461538461538459



# 使用sklearn实现SVM支持向量机

wangchengjun@nju.edu.cn

• 将每个数据在N维空间中用点标出（N是你所有的特征总数），每个特征的值是一个坐标的值。
• 举个例子，如果我们只有身高和头发长度两个特征，我们会在二维空间中标出这两个变量，每个点有两个坐标（这些坐标叫做支持向量）。

• 现在，我们会找到将两组不同数据分开的一条直线。
• 两个分组中距离最近的两个点到这条线的距离同时最优化。

## 上面示例中的黑线将数据分类优化成两个小组

• 两组中距离最近的点（图中A、B点）到达黑线的距离满足最优条件。
• 这条直线就是我们的分割线。接下来，测试数据落到直线的哪一边，我们就将它分到哪一类去。


In :

from sklearn import svm
# Create SVM classification object
model=svm.SVC()




In :

' '.join(dir(svm))




Out:

'LinearSVC LinearSVR NuSVC NuSVR OneClassSVM SVC SVR __all__ __builtins__ __doc__ __file__ __name__ __package__ __path__ base bounds classes l1_min_c liblinear libsvm libsvm_sparse'




In :

data_X_train, data_X_test, data_y_train, data_y_test = randomSplitLogistic(data_X, df.repost, 20)
model.fit(data_X_train,data_y_train)
model.score(data_X_train,data_y_train)




Out:

0.90156599552572703




In :

# Predict
model.predict(data_X_test)




Out:

array([0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1])




In :

# crossvalidation
scores = []
cvs = [3, 5, 10, 25, 50, 75, 100]
for i in cvs:
score = cross_val_score(model, data_X, df.repost, cv = i)
scores.append(score.mean() ) # Try to tune cv




In :

plt.plot(cvs, scores, 'b-o')
plt.xlabel('$cv$', fontsize = 20)
plt.ylabel('$Score$', fontsize = 20)
plt.show()






# 泰坦尼克号数据分析

wangchengjun@nju.edu.cn



In :

#Import the Numpy library
import numpy as np
#Import 'tree' from scikit-learn library
from sklearn import tree




In :

import pandas as pd
train = pd.read_csv('/Users/chengjun/github/cjc2016/data/tatanic_train.csv', sep = ",")




In :




Out:

Unnamed: 0
PassengerId
Survived
Pclass
Name
Sex
Age
SibSp
Parch
Ticket
Fare
Cabin
Embarked

0
0
1
0
3
Braund, Mr. Owen Harris
male
22.0
1
0
A/5 21171
7.2500
NaN
S

1
1
2
1
1
Cumings, Mrs. John Bradley (Florence Briggs Th...
female
38.0
1
0
PC 17599
71.2833
C85
C

2
2
3
1
3
Heikkinen, Miss. Laina
female
26.0
0
0
STON/O2. 3101282
7.9250
NaN
S

3
3
4
1
1
Futrelle, Mrs. Jacques Heath (Lily May Peel)
female
35.0
1
0
113803
53.1000
C123
S

4
4
5
0
3
Allen, Mr. William Henry
male
35.0
0
0
373450
8.0500
NaN
S




In :

train["Age"] = train["Age"].fillna(train["Age"].median())
#Convert the male and female groups to integer form
train["Sex"][train["Sex"] == "male"] = 0
train["Sex"][train["Sex"] == "female"] = 1

#Impute the Embarked variable
train["Embarked"] = train["Embarked"].fillna('S')
#Convert the Embarked classes to integer form
train["Embarked"][train["Embarked"] == "S"] = 0
train["Embarked"][train["Embarked"] == "C"] = 1
train["Embarked"][train["Embarked"] == "Q"] = 2




/Users/chengjun/anaconda/lib/python2.7/site-packages/ipykernel/__main__.py:3: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame

See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy
app.launch_new_instance()
/Users/chengjun/anaconda/lib/python2.7/site-packages/ipykernel/__main__.py:4: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame

See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy
/Users/chengjun/anaconda/lib/python2.7/site-packages/ipykernel/__main__.py:9: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame

See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy
/Users/chengjun/anaconda/lib/python2.7/site-packages/ipykernel/__main__.py:10: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame

See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy
/Users/chengjun/anaconda/lib/python2.7/site-packages/ipykernel/__main__.py:11: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame

See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy




In :

#Create the target and features numpy arrays: target, features_one
target = train['Survived'].values
features_one = train[["Pclass", "Sex", "Age", "Fare"]].values

#Fit your first decision tree: my_tree_one
my_tree_one = tree.DecisionTreeClassifier()

my_tree_one = my_tree_one.fit(features_one, target)
#Look at the importance of the included features and print the score
print(my_tree_one.feature_importances_)
print(my_tree_one.score(features_one, target))




[ 0.12629524  0.31274009  0.2387992   0.32216546]
0.977553310887




In :

test = pd.read_csv('/Users/chengjun/github/cjc2016/data/tatanic_test.csv', sep = ",")
# Impute the missing value with the median
test.Fare = test.Fare.median()
test["Age"] = test["Age"].fillna(test["Age"].median())
#Convert the male and female groups to integer form
test["Sex"][test["Sex"] == "male"] = 0
test["Sex"][test["Sex"] == "female"] = 1

#Impute the Embarked variable
test["Embarked"] = test["Embarked"].fillna('S')
#Convert the Embarked classes to integer form
test["Embarked"][test["Embarked"] == "S"] = 0
test["Embarked"][test["Embarked"] == "C"] = 1
test["Embarked"][test["Embarked"] == "Q"] = 2

# Extract the features from the test set: Pclass, Sex, Age, and Fare.
test_features = test[["Pclass","Sex", "Age", "Fare"]].values

# Make your prediction using the test set
my_prediction = my_tree_one.predict(test_features)

# Create a data frame with two columns: PassengerId & Survived. Survived contains your predictions
PassengerId =np.array(test['PassengerId']).astype(int)
my_solution = pd.DataFrame(my_prediction, PassengerId, columns = ["Survived"])




/Users/chengjun/anaconda/lib/python2.7/site-packages/ipykernel/__main__.py:3: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame

See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy
app.launch_new_instance()
/Users/chengjun/anaconda/lib/python2.7/site-packages/ipykernel/__main__.py:6: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame

See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy
/Users/chengjun/anaconda/lib/python2.7/site-packages/ipykernel/__main__.py:7: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame

See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy
/Users/chengjun/anaconda/lib/python2.7/site-packages/ipykernel/__main__.py:12: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame

See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy
/Users/chengjun/anaconda/lib/python2.7/site-packages/ipykernel/__main__.py:13: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame

See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy
/Users/chengjun/anaconda/lib/python2.7/site-packages/ipykernel/__main__.py:14: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame

See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy




In :

print my_solution[:3]




Survived
892         0
893         0
894         1




In :

# Check that your data frame has 418 entries
print my_solution.shape




(418, 1)




In :

# Write your solution to a csv file with the name my_solution.csv
my_solution.to_csv("/Users/chengjun/github/cjc2016/data/tatanic_solution_one.csv", index_label = ["PassengerId"])




In :

# Create a new array with the added features: features_two
features_two = train[["Pclass","Age","Sex","Fare", "SibSp", "Parch", "Embarked"]].values

#Control overfitting by setting "max_depth" to 10 and "min_samples_split" to 5 : my_tree_two
max_depth = 10
min_samples_split = 5
my_tree_two = tree.DecisionTreeClassifier(max_depth = max_depth, min_samples_split = min_samples_split, random_state = 1)
my_tree_two = my_tree_two.fit(features_two, target)

#Print the score of the new decison tree
print(my_tree_two.score(features_two, target))




0.905723905724




In :

# create a new train set with the new variable
train_two = train
train_two['family_size'] = train.SibSp + train.Parch + 1

# Create a new decision tree my_tree_three
features_three = train[["Pclass", "Sex", "Age", "Fare", "SibSp", "Parch", "family_size"]].values

my_tree_three = tree.DecisionTreeClassifier()
my_tree_three = my_tree_three.fit(features_three, target)

# Print the score of this decision tree
print(my_tree_three.score(features_three, target))




0.979797979798




In :

#Import the RandomForestClassifier
from sklearn.ensemble import RandomForestClassifier

#We want the Pclass, Age, Sex, Fare,SibSp, Parch, and Embarked variables
features_forest = train[["Pclass", "Age", "Sex", "Fare", "SibSp", "Parch", "Embarked"]].values

#Building the Forest: my_forest
n_estimators = 100
forest = RandomForestClassifier(max_depth = 10, min_samples_split=2, n_estimators = n_estimators, random_state = 1)
my_forest = forest.fit(features_forest, target)

#Print the score of the random forest
print(my_forest.score(features_forest, target))

#Compute predictions and print the length of the prediction vector:test_features, pred_forest
test_features = test[["Pclass", "Age", "Sex", "Fare", "SibSp", "Parch", "Embarked"]].values
pred_forest = my_forest.predict(test_features)
print(len(test_features))
print(pred_forest[:3])




0.939393939394
418
[0 0 0]




In :

#Request and print the .feature_importances_ attribute
print(my_tree_two.feature_importances_)
print(my_forest.feature_importances_)

#Compute and print the mean accuracy score for both models
print(my_tree_two.score(features_two, target))
print(my_forest.score(features_two, target))




[ 0.14130255  0.17906027  0.41616727  0.17938711  0.05039699  0.01923751
0.0144483 ]
[ 0.10384741  0.20139027  0.31989322  0.24602858  0.05272693  0.04159232
0.03452128]
0.905723905724
0.939393939394



# 阅读材料

The "Python Machine Learning" book code repository and info resource https://github.com/rasbt/python-machine-learning-book

An Introduction to Statistical Learning (James, Witten, Hastie, Tibshirani, 2013) : Python code https://github.com/JWarmenhoven/ISLR-python

BuildingMachineLearningSystemsWithPython https://github.com/luispedro/BuildingMachineLearningSystemsWithPython



In [ ]: