# A6 E4 X( b- t5 ]2 k1 N9 B5 ]+ ~: z2 u6 O/ q( o) {
从上图中我们可以看出,总体人群最终在玩与否的事件上被分成了四个群组。而分组是依据一些特征变量实现的。用来分组的具体指标有很多,比如Gini,information Gain, Chi-square,entropy。 2 }& g) g& m4 S4 y- X4 x, ]3 h4 }. S 3 H { z" L6 `' l, u2 h, U( h: K, v8 N# |" v. J$ M2 s8 D
from sklearn import tree1 Q" ]& {, w1 d" n
+ A W& \/ u% b( S1 J7 w! _; o | 9 G2 _" y4 M; y" t' h) [" m# Create tree object t" A: R; m. D& i# _! bmodel = tree.DecisionTreeClassifier(criterion='gini') # for classification, here you can change the algorithm as gini or entropy (information gain) by default it is gini ! w) ~6 N' N5 Z0 q! |+ t . r4 N+ M$ G$ N1 q+ b# model = tree.DecisionTreeRegressor() for regression/ B V6 _, a( [+ D3 N
& g8 c5 s% }! S
# Train the model using the training sets and check score6 ~3 e$ j- e* q: R
model.fit(X, y) ! G5 O$ o( X& i' ?model.score(X, y)2 a' g8 F/ A: L* V# ]# X
q' k( _' n! Z#Predict Output : l- w _0 E. b2 y" V7 K" {predicted= model.predict(x_test)3 S# `% }3 H* p: q
4. 支持向量机(SVM)+ ~) \0 f& H" s, c3 @4 d& A
这是一个分类算法。在这个算法中我们将每一个数据作为一个点在一个n维空间上作图(n是特征数),每一个特征值就代表对应坐标值的大小。比如说我们有两个特征:一个人的身高和发长。我们可以将这两个变量在一个二维空间上作图,图上的每个点都有两个坐标值(这些坐标轴也叫做支持向量)。0 s' U# d n' t% s6 B& g u% D
& a$ m' \( z8 L9 v0 ^现在我们要在图中找到一条直线能最大程度将不同组的点分开。两组数据中距离这条线最近的点到这条线的距离都应该是最远的。 + o& ~: b* k( ~ 4 E+ L1 r4 R7 r# }! u% l# O * S4 C l. K0 W# C. [6 ]* S5 V" r4 `. K v { e: y, s6 v9 d# ~. U0 {6 @# G9 y
在上图中,黑色的线就是最佳分割线。因为这条线到两组中距它最近的点,点A和B的距离都是最远的。任何其他线必然会使得到其中一个点的距离比这个距离近。这样根据数据点分布在这条线的哪一边,我们就可以将数据归类。" }5 _7 Z) i" |% a' G
, A* h; F( {0 y+ ?+ }5 P" Z" z#Import Library 6 s$ {! L$ h2 `2 Zfrom sklearn import svm) K5 y' F* G4 R8 R) d
#Assumed you have, X (predictor) and Y (target) for training data set and x_test(predictor) of test_dataset % M- x! O) b4 d6 G H, f# Create SVM classification object - H4 Z: T% m( o, N; z# Z/ g' [ 5 S7 @$ _ t( a. i" d! Gmodel = svm.svc() # there is various option associated with it, this is simple for classification. You can refer link, for mo# re detail. ) W( N5 W" i6 Y3 e) x: N2 K5 C- ]+ R& O8 p2 c V, k
# Train the model using the training sets and check score+ I9 g! e6 k3 F7 }$ u, G8 w9 _
model.fit(X, y)# l- m" {( N; D& S
model.score(X, y) # t. M7 Q' Z6 h) C) o4 q8 H" U7 s, t6 U( j ]
#Predict Output % o" Y$ @5 E9 p1 y9 R* _2 q* i' O2 |9 spredicted= model.predict(x_test). y1 ~- {% z% \! t
5. 朴素贝叶斯 V/ G0 O7 m: H这个算法是建立在贝叶斯理论上的分类方法。它的假设条件是自变量之间相互独立。简言之,朴素贝叶斯假定某一特征的出现与其它特征无关。比如说,如果一个水果它是红色的,圆状的,直径大概7cm左右,我们可能猜测它为苹果。即使这些特征之间存在一定关系,在朴素贝叶斯算法中我们都认为红色,圆状和直径在判断一个水果是苹果的可能性上是相互独立的。 . w3 ^9 o- G6 J5 [ + p* H7 g9 h1 g( {朴素贝叶斯的模型易于建造,并且在分析大量数据问题时效率很高。虽然模型简单,但很多情况下工作得比非常复杂的分类方法还要好。 : S& w: b; l1 G9 n$ V0 o/ n ' D* g+ D: U# e5 W! ]0 V贝叶斯理论告诉我们如何从先验概率P(c),P(x)和条件概率P(x|c)中计算后验概率P(c|x)。算法如下: , C3 v. e! M$ C* }! }8 P/ b- ^8 @. H: j* \2 l
4 t- B |6 Q5 iP(c|x)是已知特征x而分类为c的后验概率。4 O/ ~5 v0 Y! `2 L: R8 D
$ D# p% K5 C* Q
P(c)是种类c的先验概率。" q4 `1 |* G. E& }/ M! j
1 ]* \- ^# _" T
P(x|c)是种类c具有特征x的可能性。 6 ~1 U3 B8 l& T4 G T8 M 1 w9 G) q1 u# S" s0 f2 gP(x)是特征x的先验概率。 : e/ i q5 q* A) S7 S: [6 {# Z% [3 Q, T4 o
* \0 c& J" a X6 `
例子: 以下这组训练集包括了天气变量和目标变量“是否出去玩”。我们现在需要根据天气情况将人们分为两组:玩或不玩。整个过程按照如下步骤进行:- P5 B0 ?# I! t1 E% a/ O5 ^0 w
+ \; K) j z8 W
步骤1:根据已知数据做频率表' m6 A% Y3 t# p
$ ~5 v9 O3 Z- x0 D- D2 L9 w步骤2:计算各个情况的概率制作概率表。比如阴天(Overcast)的概率为0.29,此时玩的概率为0.64. 0 p3 V/ j) j3 U* {$ S/ o7 S# q8 p0 W* B# l
+ ^- C [9 b+ h; W" k
步骤3:用朴素贝叶斯计算每种天气情况下玩和不玩的后验概率。概率大的结果为预测值。 c7 K! N2 Q+ ?+ j. j* w
提问: 天气晴朗的情况下(sunny),人们会玩。这句陈述是否正确?' l: p) _/ ~1 i. @# ]
* _3 `: q$ U' ]
我们可以用上述方法回答这个问题。P(Yes | Sunny)=P(Sunny | Yes) * P(Yes) / P(Sunny)。 - A( R$ |% n: [% i8 s. t: w* M
这里,P(Sunny |Yes) = 3/9 = 0.33, P(Sunny) = 5/14 = 0.36, P(Yes)= 9/14 = 0.64。 ( o, m4 U/ Z2 w. K# p5 ]! W7 V& i& O8 \* s, |, x! D% ?: g
那么,P (Yes | Sunny) = 0.33 * 0.64 / 0.36 = 0.60>0.5,说明这个概率值更大。. x/ _9 l+ h2 a# n
X' ?% j& b) L1 k3 w& K
当有多种类别和多种特征时,预测的方法相似。朴素贝叶斯通常用于文本分类和多类别分类问题。 % p& n* o' p+ j 5 K! K' Y; N$ f. F! Y" A& L#Import Library2 c8 |5 i4 n: X% X6 y/ x; ?; C
from sklearn.naive_bayes import GaussianNB " V, C3 t7 v' L; f; g8 q4 _( T#Assumed you have, X (predictor) and Y (target) for training data set and x_test(predictor) of test_dataset ) X- L) n# h6 s7 ]* P0 K . O s8 X" f- O+ s# Create SVM classification object model = GaussianNB() # there is other distribution for multinomial classes like Bernoulli Naive Bayes, Refer link# o m. S% w: [# S9 f I7 c
- b" Z5 ?# a2 u6 @# Train the model using the training sets and check score9 L3 z3 O! }% m" Z. w: y& e- P
model.fit(X, y) ' R& ?' b; M4 {0 l! b+ m* ?4 N: H. B2 q6 b- a7 B
#Predict Output7 _( G( G/ \* {. v0 j2 S- U! p. l
predicted= model.predict(x_test) : ]! s' \, l) b# J2 U6.KNN(K-邻近算法), f8 y( ]; G8 @6 O* l* c; H
这个算法既可以解决分类问题,也可以用于回归问题,但工业上用于分类的情况更多。 KNN先记录所有已知数据,再利用一个距离函数,找出已知数据中距离未知事件最近的K组数据,最后按照这K组数据里最常见的类别预测该事件。1 W$ @3 q2 F" Q4 o \' K
; \. ~/ \8 d. i I [) k9 e& H距离函数可以是欧式距离,曼哈顿距离,闵氏距离 (Minkowski Distance), 和汉明距离(Hamming Distance)。前三种用于连续变量,汉明距离用于分类变量。如果K=1,那问题就简化为根据最近的数据分类。K值的选取时常是KNN建模里的关键。 ; q. ?% n! t6 p! ^6 ?% W" e* V z: o6 z7 U
9 x. F; L0 y- V( r( f! T所有特征应该标准化数量级,否则数量级大的特征在计算距离上会有偏移。8 q4 m3 v- W9 [- `7 S
; P& W: t$ q! C& a6 S. I- o+ U
在进行KNN前预处理数据,例如去除异常值,噪音等。- _6 H* T1 b2 E6 `( `8 E7 t
$ l* p X* c7 f) K# Z5 N
#Import Library ! e4 ]/ I8 |+ O# ?from sklearn.neighbors import KNeighborsClassifier$ G h( m! C# G) e% Y
! j! L0 r6 H) |: P#Assumed you have, X (predictor) and Y (target) for training data set and x_test(predictor) of test_dataset- q' v# O6 ]: B& z' P4 d* _4 K% l4 K
# Create KNeighbors classifier object model " I% r. |' j. y2 E
( D, s) k0 f1 f% cKNeighborsClassifier(n_neighbors=6) # default value for n_neighbors is 5 # A, ?9 P0 L+ e9 j' W) ~8 U: ]5 B9 q, E' K- _0 V
# Train the model using the training sets and check score' D6 I0 n* d9 \
model.fit(X, y)& A$ C( f4 P9 |3 m. M" Q
4 [: i! C8 b0 j" B$ z
#Predict Output 1 \' O& n% s5 r, w1 _1 p+ B8 Epredicted= model.predict(x_test)- C, e- l5 g( l* L% n* `% n7 c
7. K均值算法(K-Means); e& k# ^& g. T6 O* J
这是一种解决聚类问题的非监督式学习算法。这个方法简单地利用了一定数量的集群(假设K个集群)对给定数据进行分类。同一集群内的数据点是同类的,不同集群的数据点不同类。* I! K, G- [8 m. Q& q
' h2 B( P; ?+ D/ {) Q9 k还记得你是怎样从墨水渍中辨认形状的么?K均值算法的过程类似,你也要通过观察集群形状和分布来判断集群数量! ; K0 y( r+ p$ K" t2 x+ f3 s( ^) a 8 @9 H. E# w; K3 j! B" C, p; L r/ G% K( H# F7 U4 V+ k
K均值算法如何划分集群: 5 F+ ?5 p1 @- _# i) r8 ~: A. m0 f( e6 @; K8 U
0 l2 E N! Y9 ^1 }& D/ a1 K, }
; s1 ~4 L! V& g6 s0 [
从每个集群中选取K个数据点作为质心(centroids)。& t2 k# U" K: I) y) A! K) R" L: U$ P: m
9 r: z% `) P' x
将每一个数据点与距离自己最近的质心划分在同一集群,即生成K个新集群。4 s! \7 l2 X4 M; B
, Z# s2 k, v4 a找出新集群的质心,这样就有了新的质心。 0 w8 ~6 v# P" o9 @% M' E- V5 B" ~: M: U$ v; F
重复2和3,直到结果收敛,即不再有新的质心出现。 : I: q2 l5 e. B& T! G 1 A# I; T9 z7 u( \; U9 ]; N1 H2 Z! V0 r A: ^
怎样确定K的值: 8 S/ s' k' C% C4 [5 D, ^4 ]& G # v" w: p1 {. X: c0 y如果我们在每个集群中计算集群中所有点到质心的距离平方和,再将不同集群的距离平方和相加,我们就得到了这个集群方案的总平方和。 8 s5 \1 O2 Q+ J/ g% ^9 V. ^, Z% `4 d" E, J+ S+ F3 s9 u
我们知道,随着集群数量的增加,总平方和会减少。但是如果用总平方和对K作图,你会发现在某个K值之前总平方和急速减少,但在这个K值之后减少的幅度大大降低,这个值就是最佳的集群数。- a! v, S. {6 X4 V6 j2 k8 B
! O K: u8 U" U8 B$ A
+ E# t! c+ B. V. f
#Import Library! U6 L. Q+ O0 X$ A/ k8 y8 x) y
from sklearn.cluster import KMeans& b- ~& ~- H L$ B8 o* x5 y
; k+ u1 l5 ~! `% {. R#Assumed you have, X (attributes) for training data set and x_test(attributes) of test_dataset 2 E0 a" v( h+ F6 g# Create KNeighbors classifier object model " ~0 y$ C& c+ a# z6 t# I# x& I
k_means = KMeans(n_clusters=3, random_state=0)* R; R8 [6 e8 _5 I" e* q5 c! f
6 k4 B1 B3 r# j9 i* y+ N# Train the model using the training sets and check score" x( c. i& B0 C& j; s6 f' N
model.fit(X) E* ~& h- F+ |, ~- a
3 j% h5 D% C: [5 p#Predict Output! n0 f/ ^2 w3 j
predicted= model.predict(x_test) 4 }2 ?' g4 {9 h& i0 R8.随机森林% p* T/ I0 S B8 n5 H& P
随机森林是对决策树集合的特有名称。随机森林里我们有多个决策树(所以叫“森林”)。为了给一个新的观察值分类,根据它的特征,每一个决策树都会给出一个分类。随机森林算法选出投票最多的分类作为分类结果。 ! U: Q( b+ n& y. C+ w: o, ?8 M / p/ e' f; W, Q; j- y+ ~怎样生成决策树: ; k P5 e5 o" Q g% \. e' H3 L' s I8 B5 {# I
如果训练集中有N种类别,则有重复地随机选取N个样本。这些样本将组成培养决策树的训练集。 0 _, x8 ? V+ V0 z5 E# J* @1 t% Y3 i8 ]$ Z( ?
如果有M个特征变量,那么选取数m << M,从而在每个节点上随机选取m个特征变量来分割该节点。m在整个森林养成中保持不变。& t& V" v2 @8 O# V H8 x
' H( d0 l" {; o8 @* J! n每个决策树都最大程度上进行分割,没有剪枝。 * H' V) k0 _: v, [9 V7 [ . t4 O7 Y' C/ B3 v+ k5 N#Import Library" n: L# f% t* e
from sklearn.ensemble import RandomForestClassifier2 y9 N: K9 P# ]9 j( P$ O1 d
#Assumed you have, X (predictor) and Y (target) for training data set and x_test(predictor) of test_dataset }% Y2 `7 t- S; F1 B 1 B8 Q6 G2 h+ y: T. C! |2 ~+ F# Create Random Forest object8 v1 i; \0 a% u9 `; g
model= RandomForestClassifier(): C3 l5 K; \# h- Y* u1 h
; }7 V. r0 d: w* T# b
# Train the model using the training sets and check score+ ]5 N$ d! ]2 @; W+ q8 l
model.fit(X, y)) U5 O+ s/ {/ P' ~- h
9 j6 a. N9 `8 }2 G/ C; c1 u#Predict Output9 J( f4 ^! u* Q8 S3 q y
predicted= model.predict(x_test)9 Z$ e m( A; K; ~3 r$ R7 |
9.降维算法(Dimensionality Reduction Algorithms). H- ?) H1 s% `, H1 j- A _
在过去的4-5年里,可获取的数据几乎以指数形式增长。公司/政府机构/研究组织不仅有了更多的数据来源,也获得了更多维度的数据信息。 + s+ E$ [: W; f3 m& @# ^8 o8 ]- x, C4 _3 b! v( V
例如:电子商务公司有了顾客更多的细节信息,像个人信息,网络浏览历史,个人喜恶,购买记录,反馈信息等,他们关注你的私人特征,比你天天去的超市里的店员更了解你。 % a3 D1 w) }8 N& Q0 S' m- H7 n1 [ : H; A' v3 T8 w2 q作为一名数据科学家,我们手上的数据有非常多的特征。虽然这听起来有利于建立更强大精准的模型,但它们有时候反倒也是建模中的一大难题。怎样才能从1000或2000个变量里找到最重要的变量呢?这种情况下降维算法及其他算法,如决策树,随机森林,PCA,因子分析,相关矩阵,和缺省值比例等,就能帮我们解决难题。2 q% }1 g M7 g. W4 d% e
/ {' J9 t% F& r0 D; ]3 p
6 j- T R5 U$ _$ p9 h- {#Import Library& @: [: O2 Y( |4 l. ~- Q4 T5 r/ Q
from sklearn import decomposition" M. A3 \* p6 P" y
#Assumed you have training and test data set as train and test) H0 Y) @1 [9 {
# Create PCA obeject pca= decomposition.PCA(n_components=k) #default value of k =min(n_sample, n_features)6 z# o: L4 E+ V p. H
# For Factor analysis 1 P* i1 ?7 r- n- l#fa= decomposition.FactorAnalysis(); p, y/ P# _) J0 S, o! m- d
# Reduced the dimension of training dataset using PCA ( C9 v0 F# Z8 L) _2 U; G2 i1 m u- |& Z
train_reduced = pca.fit_transform(train)5 G: U- r; l/ m) u
* j S% [* o; A6 P, q5 h; ]
#Reduced the dimension of test dataset . Z2 r$ R+ D6 W! S, Ttest_reduced = pca.transform(test)" f* d* c) B* o. `$ ~$ }
10.Gradient Boosing 和 AdaBoost * R- z& \ Z, P7 A+ LGBM和AdaBoost都是在有大量数据时提高预测准确度的boosting算法。Boosting是一种集成学习方法。它通过有序结合多个较弱的分类器/估测器的估计结果来提高预测准确度。这些boosting算法在Kaggle,AV Hackthon, CrowdAnalytix等数据科学竞赛中有出色发挥。 2 q* H8 R7 n2 q+ X+ B- \6 O+ l! [, A A" U' F/ G0 r- _$ V* P/ D
#Import Library ) f1 c& B) ?% `% k2 r8 ofrom sklearn.ensemble import GradientBoostingClassifier 8 Y' H! o: N. q9 q( T8 G#Assumed you have, X (predictor) and Y (target) for training data set and x_test(predictor) of test_dataset & Z9 L/ L: D, `1 ]7 X# Create Gradient Boosting Classifier object# Z4 t) x+ N& e7 h* z
model= GradientBoostingClassifier(n_estimators=100, learning_rate=1.0, max_depth=1, random_state=0) 1 p* R3 T6 G* ?8 u7 u- t) S. E
# Train the model using the training sets and check score' Z# r# L# M" e. c9 ^4 w5 z' ]) ]8 x
model.fit(X, y) ' l, [6 x/ f8 O#Predict Output- N- h8 @( o/ ?& A/ l6 `
predicted= model.predict(x_test) ( M' V$ s" }# u2 j7 n- qGradientBoostingClassifier 和随机森林是两种不同的boosting分类树。人们经常提问 这两个算法有什么不同。. o) @" o2 d# L# a
8 R4 C2 g: }" F8 |+ G: u( _8 r
原文链接:http://blog.csdn.net/han_xiaoyang/article/details/51191386 ; v; b1 g: p/ u& t% {1 Q———————————————— s: [9 `" S# m& a2 |; Z
版权声明:本文为CSDN博主「_小羊」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。 # N$ z0 E' G$ j- T0 L6 v原文链接:https://blog.csdn.net/qq_39303465/article/details/79176075 5 A* x$ K$ B9 [) E& Q8 l- q8 q" X+ O5 }( R0 v' j) l0 t$ F
0 i( o7 j5 p7 J