http://wiki.math.uwaterloo.ca/statwiki/api.php?action=feedcontributions&user=Z66song&feedformat=atomstatwiki - User contributions [US]2022-01-22T14:33:59ZUser contributionsMediaWiki 1.28.3http://wiki.math.uwaterloo.ca/statwiki/index.php?title=summary&diff=34517summary2018-03-17T15:37:31Z<p>Z66song: /* Split Finding Algorithms */</p>
<hr />
<div>'''''XGBoost: A Scalable Tree Boosting System'''''<br />
<br />
== Presented by ==<br />
<br />
Jiang, Cong<br />
<br />
Song, Ziwei<br />
<br />
Ye, Zhaoshan<br />
<br />
Zhang, Wenling<br />
<br />
= Introduction =<br />
<br />
The Extreme Gradiant Boosting, also known as XGBoost, has been the leading model that dominates Kaggle, as well as other machine learning and data mining challenges. The main reason to its success is that this tree boosting system is highly scalable , which means it could be used to solve various problems with significantly less time and fewer resources. <br />
<br />
This paper is written by the father of XGBoost, Tianqi Chen. It gives us an overview on how he used algorithmic optimizations as well as some important systems to develop XGBoost. He explained it in the following manner:<br />
<br />
1. Based on gradient tree boosting algorithm, Chen makes some modifications to the regularized objective function and introduces an end-to-end tree boosting system. He also applies shrinkage and feature subsampling to prevent overfitting, <br />
<br />
2. The system supports existing split finding algorithms including exact greedy and approximate algorithms. For weighted data, Chen proposes a distributed weighted quantile sketch algorithm which is provable; while for sparse data, he develops a novel sparsity-aware algorithm which is much faster than the naive model.<br />
<br />
3. Storing the data in ''blocks'' is Chen's method to reduce sorting cost. The cache-aware block structure is highly effective for out-of-core tree learning.<br />
<br />
4. XGBoost is evaluated by designed experiments with 4 datasets. Chen compares it with adjusted setting as well as other systems to test its performance.<br />
<br />
= Tree Boosting=<br />
<br />
<br />
=== Regularized Learning Objective ===<br />
Suppose the given data set has m features and n examples. A tree ensemble model uses K additive functions to predict the output.<br />
[[File:eq1.png|center|eq1.png]]<br />
where F is the space of regression trees(CART).<br />
The author suggests the following regularized objective function to make improvements comparing to gradient boosting<br />
[[File:eq2.png|center|eq2.png]]<br />
Regression tree contains continuous scores on each leaf, which is different with decision trees we have learned in class. The following 2 figures represent how to use the decision rules in the trees to classify it into the leaves and calculate the final prediction by summing up the score in the corresponding leaves.<br />
<br />
[[File:Picture2.1.png|500px| thumb| center| Figure 0: Prediction Score of Tree Structure.]]<br />
<br />
<br />
[[File:Picture2.2.png|500px|thumb|center|Figure 1: Tree Ensemble Model. The final prediction for a given example is the sum of predictions from each tree. ]]<br />
<br />
=== Gradient Tree Boosting ===<br />
The tree ensemble model (Equation 2.) is trained in an additive manner. We need to add a new <math>f</math> to our objective.<br />
[[File:eq2.2.png|center|eq2.2.png]]<br />
After applying the second-order Taylor expansion to the above formula and removing constant terms, we can finally have<br />
[[File:eq3.png|center|eq3.png]]<br />
Define <math>I_{j} = \{i| q(X_{i}) \}</math><br />
As the instance set of leaf j. We can rewrite Equation (3) as<br />
[[File:eq4.png|center|eq4.png]]<br />
By computing the optimal weight and corresponding optimal value, we can use following equation as the scoring function which can measure the quality of a tree structure. Lower score means better tree structure.<br />
[[File:eq5.png|center|eq5.png]]<br />
The following figure shows an example how this score is calculated.<br />
[[File:Picture2.3.png|500px| thumb|center|Figure 2: Structure Score Calculation. We only need to sum up the gradient and second order gradient statistics on each leaf, then apply the scoring formula to get the quality score. ]]<br />
<br />
<br />
<br />
Normally it is impossible to enumerate all the possible tree structures. Assume <math>I_{L}</math> and <math>I_{R}</math> are the instance sets of left and right nodes after the split, and <math> I = I_{L} + I_{R}</math>. Then the loss reduction after the split is given by<br />
[[File:eq6.png|center]]<br />
<br />
=== Shrinkage and Column Subsampling ===<br />
In addition to the regularized objective, the paper also introduces two techniques to prevent over-fitting. <br />
<br />
Method 1: Shrinkage. After each step of tree boosting, we scale newly added weights by a shrinkage factor η, which reduces the influence of every individual tree and leaves the space for future trees. <br />
<br />
Method 2: Column (feature) subsampling. We try to build each individual tree with only a portion of predictors chosen randomly. This procedure encourages the variance between the trees (i.e. no same tree occurs). Implemented in an opensource packages, this method prevents over-fitting, compared with the traditional row subsampling, and speeds up the computations of parallel algorithm as well.<br />
<br />
= Split Finding Algorithms =<br />
The ways to find the best split can be realized by two main methods. <br />
<br />
Basic exact greedy algorithm: enumerate over all the possible splits on all the features, which is computationally demanding. Sort the data first, and visit the data to accumulate the gradient statistics greedily. <br />
<br />
[[File:exact_greedy.png|500px| thumb|center ]]<br />
[[File:appro_greedy.png|500px| thumb|center]]<br />
<br />
Approximate algorithm: To achieve an effective gradient tree boosting model, The author introduces an approximate algorithm. The algorithm first proposes candidate splitting points by the percentiles of feature distribution, maps the continuous features into buckets split and then finds the best solution among the proposals by the aggregated statistics. <br />
The algorithm can be varied by the time the proposal is given. <br />
<br />
The global variant: all candidate splits are proposed at the initial phase of tree construction and are used for all split finding. It requires fewer proposal steps. <br />
<br />
The local variant: the candidate splits are refined after each split. It requires fewer candidate points, so it is potentially more appropriate for deeper trees. <br />
Methods can be chosen according to users’ needs. <br />
<br />
[[File:auc.png|500px| thumb|center]]<br />
<br />
=== Weighted Quantile Sketch ===<br />
<br />
=== Sparsity--aware Split Finding ===<br />
<br />
In real world problem, it is quite common for the input <math> X </math> to be sparse. The reasons that causes for sparsity are following: 1. presence of missing values in the data; 2.frequent zero entries in the statistics; and, 3.artifacts of feature engineering such as one-hot encoding.<br />
<br />
[[File:Figure 4.png|500px| thumb|center | Figure 4. Tree structure with default directions. An example will be classified into the default direction when the feature needed for the split is missing]]<br />
This paper propose to add a default direction in each tree node, which is shown in Figure 4. When a value is missing in the sparse matrix <math> X </math>, the instance is classified into the default direction. There are two choices of default direction in each branch. The optimal default direction is learn from data. Algorithm 3 shows Sparsity aware Split Finding. <br />
<br />
[[File:Algorithm 3.png|500px| thumb|center ]]<br />
The important improvement is to only visit the non-missing entries <math> I_k</math>. The presented algorithm treats the non-presence as a missing value and learns the best direction to handle missing values.<br />
<br />
=Evaluations=<br />
===System Implementation===<br />
Implemented as an open source package, XGBoost is flexible, portable and reusable. It also supports various languages and ecosystems.<br />
<br />
===Dataset and Setup===<br />
Table 2 summarizes the 4 datasets used in the experiments.<br />
[[File:t2.png|500px| thumb|center ]]<br />
<br />
<br />
The first three datasets are used for the single machine parallel setting, which is run on a Dell PowerEdge R420 with two 8-core Intel Xeon and memory 64GB. Meanwhile, AWS c3.8xlarge machine is used for experiment in distributed and the out-of-core setting, which uses the criteo dataset. All available cores in the machines are devoted into the experiments, unless otherwise specified.<br />
<br />
=== Experiments===<br />
The classification experiment uses Higgs Boson dataset to classify whether an event corresponds to the Higgs boson. The goal is to evaluate XGBoost using Exact Greedy Methods and compare it with column sampling XGBoost and two other well-known systems. Table 3 below summarizes the results.<br />
[[File:t3.png|500px| thumb|center ]]<br />
Overall, both XGBoost and scikit-learn performs better than R.gbm by 20% of Test AUC. In terms of time, however, XGBoost runs significantly faster than both scikit-learn and R.gbm. Note that XGBoost with column sampling performs slightly worse because the number of key features is limited. We could conclude that XGBoost using the exact greedy algorithm without column sampling works the best under this scenario.<br />
<br />
The learning to rank experiment compares XGBoost(using exact greedy algorithm) with pGBRT(using approximate algorithm) in order to evaluate the system's performance on rank problem.<br />
<br />
[[File:t4.png|500px| thumb|center ]]<br />
[[File:f10.png|500px| thumb|center ]]<br />
According to Table 4, XGBoost with column sampling performs almost as well as pGBRT, while using only 1/5 of the time. Note that pGBRT is the best previously published system to solve this rank problem. <br />
<br />
The out-of-score and distributed experiments demonstrate XGBoost's ability to handle large data with a fast speed.<br />
<br />
= Conclusion=<br />
The writer, Tianqi Chen, proposed a scalable end-to-end tree boosting system with cache access patterns, block structure and other essential elements. He also introduced two novel split finding algorithms, one is sparsity-aware algorithm and the other is weighted quantile sketch with proofs. Compared to other popular systems, XGBoost gives state-of-the-art results in different experiment settings, which explains why it it considered one of the best and fastest algorithm for structured or tabular data.<br />
<br />
= Source =</div>Z66songhttp://wiki.math.uwaterloo.ca/statwiki/index.php?title=summary&diff=34504summary2018-03-17T14:08:09Z<p>Z66song: /* Split Finding Algorithms */</p>
<hr />
<div>'''''XGBoost: A Scalable Tree Boosting System'''''<br />
<br />
== Presented by ==<br />
<br />
Jiang, Cong<br />
<br />
Song, Ziwei<br />
<br />
Ye, Zhaoshan<br />
<br />
Zhang, Wenling<br />
<br />
= Introduction =<br />
<br />
The Extreme Gradiant Boosting, also known as XGBoost, has been the leading model that dominates Kaggle, as well as other machine learning and data mining challenges. The main reason to its success is that this tree boosting system is highly scalable , which means it could be used to solve various problems with significantly less time and fewer resources. <br />
<br />
This paper is written by the father of XGBoost, Tianqi Chen. It gives us an overview on how he used algorithmic optimizations as well as some important systems to develop XGBoost. He explained it in the following manner:<br />
<br />
1. Based on gradient tree boosting algorithm, Chen makes some modifications to the regularized objective function and introduces an end-to-end tree boosting system. He also applies shrinkage and feature subsampling to prevent overfitting, <br />
<br />
2. The system supports existing split finding algorithms including exact greedy and approximate algorithms. For weighted data, Chen proposes a distributed weighted quantile sketch algorithm which is provable; while for sparse data, he develops a novel sparsity-aware algorithm which is much faster than the naive model.<br />
<br />
3. Storing the data in ''blocks'' is Chen's method to reduce sorting cost. The cache-aware block structure is highly effective for out-of-core tree learning.<br />
<br />
4. XGBoost is evaluated by designed experiments with 4 datasets. Chen compares it with adjusted setting as well as other systems to test its performance.<br />
<br />
= Tree Boosting=<br />
<br />
<br />
=== Regularized Learning Objective ===<br />
Suppose the given data set has m features and n examples. A tree ensemble model uses K additive functions to predict the output.<br />
[[File:eq1.png|center|eq1.png]]<br />
where F is the space of regression trees(CART).<br />
The author suggests the following regularized objective function to make improvements comparing to gradient boosting<br />
[[File:eq2.png|center|eq2.png]]<br />
Regression tree contains continuous scores on each leaf, which is different with decision trees we have learned in class. The following 2 figures represent how to use the decision rules in the trees to classify it into the leaves and calculate the final prediction by summing up the score in the corresponding leaves.<br />
<br />
[[File:Picture2.1.png|500px| thumb| center| Figure 0: Prediction Score of Tree Structure.]]<br />
<br />
<br />
[[File:Picture2.2.png|500px|thumb|center|Figure 1: Tree Ensemble Model. The final prediction for a given example is the sum of predictions from each tree. ]]<br />
<br />
=== Gradient Tree Boosting ===<br />
The tree ensemble model (Equation 2.) is trained in an additive manner. We need to add a new <math>f</math> to our objective.<br />
[[File:eq2.2.png|center|eq2.2.png]]<br />
After applying the second-order Taylor expansion to the above formula and removing constant terms, we can finally have<br />
[[File:eq3.png|center|eq3.png]]<br />
Define <math>I_{j} = \{i| q(X_{i}) \}</math><br />
As the instance set of leaf j. We can rewrite Equation (3) as<br />
[[File:eq4.png|center|eq4.png]]<br />
By computing the optimal weight and corresponding optimal value, we can use following equation as the scoring function which can measure the quality of a tree structure. Lower score means better tree structure.<br />
[[File:eq5.png|center|eq5.png]]<br />
The following figure shows an example how this score is calculated.<br />
[[File:Picture2.3.png|500px| thumb|center|Figure 2: Structure Score Calculation. We only need to sum up the gradient and second order gradient statistics on each leaf, then apply the scoring formula to get the quality score. ]]<br />
<br />
<br />
<br />
Normally it is impossible to enumerate all the possible tree structures. Assume <math>I_{L}</math> and <math>I_{R}</math> are the instance sets of left and right nodes after the split, and <math> I = I_{L} + I_{R}</math>. Then the loss reduction after the split is given by<br />
[[File:eq6.png|center]]<br />
<br />
=== Shrinkage and Column Subsampling ===<br />
In addition to the regularized objective, the paper also introduces two techniques to prevent over-fitting. <br />
<br />
Method 1: Shrinkage. After each step of tree boosting, we scale newly added weights by a shrinkage factor η, which reduces the influence of every individual tree and leaves the space for future trees. <br />
<br />
Method 2: Column (feature) subsampling. We try to build each individual tree with only a portion of predictors chosen randomly. This procedure encourages the variance between the trees (i.e. no same tree occurs). Implemented in an opensource packages, this method prevents over-fitting, compared with the traditional row subsampling, and speeds up the computations of parallel algorithm as well.<br />
<br />
= Split Finding Algorithms =<br />
The ways to find the best split can be realized by two main methods. <br />
<br />
Basic exact greedy algorithm: enumerate over all the possible splits on all the features, which is computationally demanding. Sort the data first, and visit the data to accumulate the gradient statistics greedily. <br />
<br />
[[File:exact_greedy.png|500px| thumb|center ]]<br />
[[File:appro_greedy.png|500px| thumb|center]]<br />
<br />
Approximate algorithm: To achieve an effective gradient tree boosting model, The author introduces an approximate algorithm. The algorithm first proposes candidate splitting points by the percentiles of feature distribution, maps the continuous features into buckets split and then finds the best solution among the proposals by the aggregated statistics. <br />
The algorithm can be varied by the time the proposal is given. <br />
The global variant: all candidate splits are proposed at the initial phase of tree construction and are used for all split finding. It requires fewer proposal steps. <br />
The local variant: the candidate splits are refined after each split. It requires fewer candidate points, so it is potentially more appropriate for deeper trees. <br />
Methods can be chosen according to users’ needs. <br />
<br />
[[File:auc.png|500px| thumb|center]]<br />
<br />
=== Weighted Quantile Sketch ===<br />
<br />
=== Sparsity--aware Split Finding ===<br />
<br />
In real world problem, it is quite common for the input X to be sparse. The reasons that causes for sparsity are following: 1. presence of missing values in the data; 2.frequent zero entries in the statistics; and, 3.artifacts of feature engineering such as one-hot encoding.<br />
<br />
=Evaluations=<br />
===System Implementation===<br />
Implemented as an open source package, XGBoost is flexible, portable and reusable. It also supports various languages and ecosystems.<br />
<br />
===Dataset and Setup===<br />
Table 2 summarizes the 4 datasets used in the experiments.<br />
[[File:t2.png|500px| thumb|center ]]<br />
<br />
<br />
The first three datasets are used for the single machine parallel setting, which is run on a Dell PowerEdge R420 with two 8-core Intel Xeon and memory 64GB. Meanwhile, AWS c3.8xlarge machine is used for experiment in distributed and the out-of-core setting, which uses the criteo dataset. All available cores in the machines are devoted into the experiments, unless otherwise specified.<br />
<br />
=== Experiments===<br />
The classification experiment uses Higgs Boson dataset to classify whether an event corresponds to the Higgs boson. The goal is to evaluate XGBoost using Exact Greedy Methods and compare it with column sampling XGBoost and two other well-known systems. Table 3 below summarizes the results.<br />
[[File:t3.png|500px| thumb|center ]]<br />
Overall, both XGBoost and scikit-learn performs better than R.gbm by 20% of Test AUC. In terms of time, however, XGBoost runs significantly faster than both scikit-learn and R.gbm. Note that XGBoost with column sampling performs slightly worse because the number of key features is limited. We could conclude that XGBoost using the exact greedy algorithm without column sampling works the best under this scenario.<br />
<br />
The learning to rank experiment compares XGBoost(using exact greedy algorithm) with pGBRT(using approximate algorithm) in order to evaluate the system's performance on rank problem.<br />
<br />
[[File:t4.png|500px| thumb|center ]]<br />
[[File:f10.png|500px| thumb|center ]]<br />
According to Table 4, XGBoost with column sampling performs almost as well as pGBRT, while using only 1/5 of the time. Note that pGBRT is the best previously published system to solve this rank problem. <br />
<br />
The out-of-score and distributed experiments demonstrate XGBoost's ability to handle large data with a fast speed.<br />
<br />
= Conclusion=<br />
The writer, Tianqi Chen, proposed a scalable end-to-end tree boosting system with cache access patterns, block structure and other essential elements. He also introduced two novel split finding algorithms, one is sparsity-aware algorithm and the other is weighted quantile sketch with proofs. Compared to other popular systems, XGBoost gives outstanding performance in different experiment settings, which explains why it is considered the best and fastest algorithm for structured or tabular data.<br />
<br />
= Source =</div>Z66songhttp://wiki.math.uwaterloo.ca/statwiki/index.php?title=File:auc.png&diff=34503File:auc.png2018-03-17T14:06:50Z<p>Z66song: </p>
<hr />
<div></div>Z66songhttp://wiki.math.uwaterloo.ca/statwiki/index.php?title=summary&diff=34502summary2018-03-17T14:03:55Z<p>Z66song: /* Split Finding Algorithms */</p>
<hr />
<div>'''''XGBoost: A Scalable Tree Boosting System'''''<br />
<br />
== Presented by ==<br />
<br />
Jiang, Cong<br />
<br />
Song, Ziwei<br />
<br />
Ye, Zhaoshan<br />
<br />
Zhang, Wenling<br />
<br />
= Introduction =<br />
<br />
The Extreme Gradiant Boosting, also known as XGBoost, has been the leading model that dominates Kaggle, as well as other machine learning and data mining challenges. The main reason to its success is that this tree boosting system is highly scalable , which means it could be used to solve various problems with significantly less time and fewer resources. <br />
<br />
This paper is written by the father of XGBoost, Tianqi Chen. It gives us an overview on how he used algorithmic optimizations as well as some important systems to develop XGBoost. He explained it in the following manner:<br />
<br />
1. Based on gradient tree boosting algorithm, Chen makes some modifications to the regularized objective function and introduces an end-to-end tree boosting system. He also applies shrinkage and feature subsampling to prevent overfitting, <br />
<br />
2. The system supports existing split finding algorithms including exact greedy and approximate algorithms. For weighted data, Chen proposes a distributed weighted quantile sketch algorithm which is provable; while for sparse data, he develops a novel sparsity-aware algorithm which is much faster than the naive model.<br />
<br />
3. Storing the data in ''blocks'' is Chen's method to reduce sorting cost. The cache-aware block structure is highly effective for out-of-core tree learning.<br />
<br />
4. XGBoost is evaluated by designed experiments with 4 datasets. Chen compares it with adjusted setting as well as other systems to test its performance.<br />
<br />
= Tree Boosting=<br />
<br />
<br />
=== Regularized Learning Objective ===<br />
Suppose the given data set has m features and n examples. A tree ensemble model uses K additive functions to predict the output.<br />
[[File:eq1.png|center|eq1.png]]<br />
where F is the space of regression trees(CART).<br />
The author suggests the following regularized objective function to make improvements comparing to gradient boosting<br />
[[File:eq2.png|center|eq2.png]]<br />
Regression tree contains continuous scores on each leaf, which is different with decision trees we have learned in class. The following 2 figures represent how to use the decision rules in the trees to classify it into the leaves and calculate the final prediction by summing up the score in the corresponding leaves.<br />
<br />
[[File:Picture2.1.png|500px| thumb| center| Figure 0: Prediction Score of Tree Structure.]]<br />
<br />
<br />
[[File:Picture2.2.png|500px|thumb|center|Figure 1: Tree Ensemble Model. The final prediction for a given example is the sum of predictions from each tree. ]]<br />
<br />
=== Gradient Tree Boosting ===<br />
The tree ensemble model (Equation 2.) is trained in an additive manner. We need to add a new <math>f</math> to our objective.<br />
[[File:eq2.2.png|center|eq2.2.png]]<br />
After applying the second-order Taylor expansion to the above formula and removing constant terms, we can finally have<br />
[[File:eq3.png|center|eq3.png]]<br />
Define <math>I_{j} = \{i| q(X_{i}) \}</math><br />
As the instance set of leaf j. We can rewrite Equation (3) as<br />
[[File:eq4.png|center|eq4.png]]<br />
By computing the optimal weight and corresponding optimal value, we can use following equation as the scoring function which can measure the quality of a tree structure. Lower score means better tree structure.<br />
[[File:eq5.png|center|eq5.png]]<br />
The following figure shows an example how this score is calculated.<br />
[[File:Picture2.3.png|500px| thumb|center|Figure 2: Structure Score Calculation. We only need to sum up the gradient and second order gradient statistics on each leaf, then apply the scoring formula to get the quality score. ]]<br />
<br />
<br />
<br />
Normally it is impossible to enumerate all the possible tree structures. Assume <math>I_{L}</math> and <math>I_{R}</math> are the instance sets of left and right nodes after the split, and <math> I = I_{L} + I_{R}</math>. Then the loss reduction after the split is given by<br />
[[File:eq6.png|center]]<br />
<br />
=== Shrinkage and Column Subsampling ===<br />
In addition to the regularized objective, the paper also introduces two techniques to prevent over-fitting. <br />
<br />
Method 1: Shrinkage. After each step of tree boosting, we scale newly added weights by a shrinkage factor η, which reduces the influence of every individual tree and leaves the space for future trees. <br />
<br />
Method 2: Column (feature) subsampling. We try to build each individual tree with only a portion of predictors chosen randomly. This procedure encourages the variance between the trees (i.e. no same tree occurs). Implemented in an opensource packages, this method prevents over-fitting, compared with the traditional row subsampling, and speeds up the computations of parallel algorithm as well.<br />
<br />
= Split Finding Algorithms =<br />
The ways to find the best split can be realized by two main methods. <br />
<br />
Basic exact greedy algorithm: enumerate over all the possible splits on all the features, which is computationally demanding. Sort the data first, and visit the data to accumulate the gradient statistics greedily. <br />
<br />
Approximate algorithm: To achieve an effective gradient tree boosting model, The author introduces an approximate algorithm. The algorithm first proposes candidate splitting points by the percentiles of feature distribution, maps the continuous features into buckets split and then finds the best solution among the proposals by the aggregated statistics. <br />
The algorithm can be varied by the time the proposal is given. <br />
The global variant: all candidate splits are proposed at the initial phase of tree construction and are used for all split finding. It requires fewer proposal steps. <br />
The local variant: the candidate splits are refined after each split. It requires fewer candidate points, so it is potentially more appropriate for deeper trees. <br />
Methods can be chosen according to users’ needs. <br />
<br />
[[File:exact_greedy.png|500px| thumb|center ]]<br />
[[File:appro_greedy.png|500px| thumb|center]]<br />
<br />
<br />
=== Weighted Quantile Sketch ===<br />
<br />
=== Sparsity--aware Split Finding ===<br />
<br />
In real world problem, it is quite common for the input X to be sparse. The reasons that causes for sparsity are following: 1. presence of missing values in the data; 2.frequent zero entries in the statistics; and, 3.artifacts of feature engineering such as one-hot encoding.<br />
<br />
=Evaluations=<br />
===System Implementation===<br />
Implemented as an open source package, XGBoost is flexible, portable and reusable. It also supports various languages and ecosystems.<br />
<br />
===Dataset and Setup===<br />
Table 2 summarizes the 4 datasets used in the experiments.<br />
[[File:t2.png|500px| thumb|center ]]<br />
<br />
<br />
The first three datasets are used for the single machine parallel setting, which is run on a Dell PowerEdge R420 with two 8-core Intel Xeon and memory 64GB. Meanwhile, AWS c3.8xlarge machine is used for experiment in distributed and the out-of-core setting, which uses the criteo dataset. All available cores in the machines are devoted into the experiments, unless otherwise specified.<br />
<br />
=== Experiments===<br />
The classification experiment uses Higgs Boson dataset to classify whether an event corresponds to the Higgs boson. The goal is to evaluate XGBoost using Exact Greedy Methods and compare it with column sampling XGBoost and two other well-known systems. Table 3 below summarizes the results.<br />
[[File:t3.png|500px| thumb|center ]]<br />
Overall, both XGBoost and scikit-learn performs better than R.gbm by 20% of Test AUC. In terms of time, however, XGBoost runs significantly faster than both scikit-learn and R.gbm. Note that XGBoost with column sampling performs slightly worse because the number of key features is limited. We could conclude that XGBoost using the exact greedy algorithm without column sampling works the best under this scenario.<br />
<br />
The learning to rank experiment compares XGBoost(using exact greedy algorithm) with pGBRT(using approximate algorithm) in order to evaluate the system's performance on rank problem.<br />
<br />
[[File:t4.png|500px| thumb|center ]]<br />
[[File:f10.png|500px| thumb|center ]]<br />
According to Table 4, XGBoost with column sampling performs almost as well as pGBRT, while using only 1/5 of the time. Note that pGBRT is the best previously published system to solve this rank problem. <br />
<br />
The out-of-score and distributed experiments demonstrate XGBoost's ability to handle large data with a fast speed.<br />
<br />
= Conclusion=<br />
The writer, Tianqi Chen, proposed a scalable end-to-end tree boosting system with cache access patterns, block structure and other essential elements. He also introduced two novel split finding algorithms, one is sparsity-aware algorithm and the other is weighted quantile sketch with proofs. Compared to other popular systems, XGBoost gives outstanding performance in different experiment settings, which explains why it is considered the best and fastest algorithm for structured or tabular data.<br />
<br />
= Source =</div>Z66songhttp://wiki.math.uwaterloo.ca/statwiki/index.php?title=summary&diff=34500summary2018-03-17T14:02:55Z<p>Z66song: /* Shrinkage and Column Subsampling */</p>
<hr />
<div>'''''XGBoost: A Scalable Tree Boosting System'''''<br />
<br />
== Presented by ==<br />
<br />
Jiang, Cong<br />
<br />
Song, Ziwei<br />
<br />
Ye, Zhaoshan<br />
<br />
Zhang, Wenling<br />
<br />
= Introduction =<br />
<br />
The Extreme Gradiant Boosting, also known as XGBoost, has been the leading model that dominates Kaggle, as well as other machine learning and data mining challenges. The main reason to its success is that this tree boosting system is highly scalable , which means it could be used to solve various problems with significantly less time and fewer resources. <br />
<br />
This paper is written by the father of XGBoost, Tianqi Chen. It gives us an overview on how he used algorithmic optimizations as well as some important systems to develop XGBoost. He explained it in the following manner:<br />
<br />
1. Based on gradient tree boosting algorithm, Chen makes some modifications to the regularized objective function and introduces an end-to-end tree boosting system. He also applies shrinkage and feature subsampling to prevent overfitting, <br />
<br />
2. The system supports existing split finding algorithms including exact greedy and approximate algorithms. For weighted data, Chen proposes a distributed weighted quantile sketch algorithm which is provable; while for sparse data, he develops a novel sparsity-aware algorithm which is much faster than the naive model.<br />
<br />
3. Storing the data in ''blocks'' is Chen's method to reduce sorting cost. The cache-aware block structure is highly effective for out-of-core tree learning.<br />
<br />
4. XGBoost is evaluated by designed experiments with 4 datasets. Chen compares it with adjusted setting as well as other systems to test its performance.<br />
<br />
= Tree Boosting=<br />
<br />
<br />
=== Regularized Learning Objective ===<br />
Suppose the given data set has m features and n examples. A tree ensemble model uses K additive functions to predict the output.<br />
[[File:eq1.png|center|eq1.png]]<br />
where F is the space of regression trees(CART).<br />
The author suggests the following regularized objective function to make improvements comparing to gradient boosting<br />
[[File:eq2.png|center|eq2.png]]<br />
Regression tree contains continuous scores on each leaf, which is different with decision trees we have learned in class. The following 2 figures represent how to use the decision rules in the trees to classify it into the leaves and calculate the final prediction by summing up the score in the corresponding leaves.<br />
<br />
[[File:Picture2.1.png|500px| thumb| center| Figure 0: Prediction Score of Tree Structure.]]<br />
<br />
<br />
[[File:Picture2.2.png|500px|thumb|center|Figure 1: Tree Ensemble Model. The final prediction for a given example is the sum of predictions from each tree. ]]<br />
<br />
=== Gradient Tree Boosting ===<br />
The tree ensemble model (Equation 2.) is trained in an additive manner. We need to add a new <math>f</math> to our objective.<br />
[[File:eq2.2.png|center|eq2.2.png]]<br />
After applying the second-order Taylor expansion to the above formula and removing constant terms, we can finally have<br />
[[File:eq3.png|center|eq3.png]]<br />
Define <math>I_{j} = \{i| q(X_{i}) \}</math><br />
As the instance set of leaf j. We can rewrite Equation (3) as<br />
[[File:eq4.png|center|eq4.png]]<br />
By computing the optimal weight and corresponding optimal value, we can use following equation as the scoring function which can measure the quality of a tree structure. Lower score means better tree structure.<br />
[[File:eq5.png|center|eq5.png]]<br />
The following figure shows an example how this score is calculated.<br />
[[File:Picture2.3.png|500px| thumb|center|Figure 2: Structure Score Calculation. We only need to sum up the gradient and second order gradient statistics on each leaf, then apply the scoring formula to get the quality score. ]]<br />
<br />
<br />
<br />
Normally it is impossible to enumerate all the possible tree structures. Assume <math>I_{L}</math> and <math>I_{R}</math> are the instance sets of left and right nodes after the split, and <math> I = I_{L} + I_{R}</math>. Then the loss reduction after the split is given by<br />
[[File:eq6.png|center]]<br />
<br />
=== Shrinkage and Column Subsampling ===<br />
In addition to the regularized objective, the paper also introduces two techniques to prevent over-fitting. <br />
<br />
Method 1: Shrinkage. After each step of tree boosting, we scale newly added weights by a shrinkage factor η, which reduces the influence of every individual tree and leaves the space for future trees. <br />
<br />
Method 2: Column (feature) subsampling. We try to build each individual tree with only a portion of predictors chosen randomly. This procedure encourages the variance between the trees (i.e. no same tree occurs). Implemented in an opensource packages, this method prevents over-fitting, compared with the traditional row subsampling, and speeds up the computations of parallel algorithm as well.<br />
<br />
= Split Finding Algorithms =<br />
The ways to find the best split can be realized by two main methods. <br />
<br />
Basic exact greedy algorithm: enumerate over all the possible splits on all the features, which is computationally demanding. Sort the data first, and visit the data to accumulate the gradient statistics. <br />
<br />
Approximate algorithm: To achieve an effective gradient tree boosting model, we first propose candidate splitting points by percentiles of feature distribution, map the continuous features into buckets split and then finds the best solution among the proposals by the aggregated statistics. <br />
The algorithm can be varied by the time the proposal is given. <br />
The global variant: all candidate splits are proposed at the initial phase of tree construction and used at all levels. <br />
The local variant: the candidate splits are refined after each splits. It is potential for deeper trees and it requires fewer candidates. <br />
Methods can be chosen according to users’ needs.<br />
[[File:exact_greedy.png|500px| thumb|center ]]<br />
[[File:appro_greedy.png|500px| thumb|center]]<br />
<br />
<br />
=== Weighted Quantile Sketch ===<br />
<br />
=== Sparsity--aware Split Finding ===<br />
<br />
In real world problem, it is quite common for the input X to be sparse. The reasons that causes for sparsity are following: 1. presence of missing values in the data; 2.frequent zero entries in the statistics; and, 3.artifacts of feature engineering such as one-hot encoding.<br />
<br />
=Evaluations=<br />
===System Implementation===<br />
Implemented as an open source package, XGBoost is flexible, portable and reusable. It also supports various languages and ecosystems.<br />
<br />
===Dataset and Setup===<br />
Table 2 summarizes the 4 datasets used in the experiments.<br />
[[File:t2.png|500px| thumb|center ]]<br />
<br />
<br />
The first three datasets are used for the single machine parallel setting, which is run on a Dell PowerEdge R420 with two 8-core Intel Xeon and memory 64GB. Meanwhile, AWS c3.8xlarge machine is used for experiment in distributed and the out-of-core setting, which uses the criteo dataset. All available cores in the machines are devoted into the experiments, unless otherwise specified.<br />
<br />
=== Experiments===<br />
The classification experiment uses Higgs Boson dataset to classify whether an event corresponds to the Higgs boson. The goal is to evaluate XGBoost using Exact Greedy Methods and compare it with column sampling XGBoost and two other well-known systems. Table 3 below summarizes the results.<br />
[[File:t3.png|500px| thumb|center ]]<br />
Overall, both XGBoost and scikit-learn performs better than R.gbm by 20% of Test AUC. In terms of time, however, XGBoost runs significantly faster than both scikit-learn and R.gbm. Note that XGBoost with column sampling performs slightly worse because the number of key features is limited. We could conclude that XGBoost using the exact greedy algorithm without column sampling works the best under this scenario.<br />
<br />
The learning to rank experiment compares XGBoost(using exact greedy algorithm) with pGBRT(using approximate algorithm) in order to evaluate the system's performance on rank problem.<br />
<br />
[[File:t4.png|500px| thumb|center ]]<br />
[[File:f10.png|500px| thumb|center ]]<br />
According to Table 4, XGBoost with column sampling performs almost as well as pGBRT, while using only 1/5 of the time. Note that pGBRT is the best previously published system to solve this rank problem. <br />
<br />
The out-of-score and distributed experiments demonstrate XGBoost's ability to handle large data with a fast speed.<br />
<br />
= Conclusion=<br />
The writer, Tianqi Chen, proposed a scalable end-to-end tree boosting system with cache access patterns, block structure and other essential elements. He also introduced two novel split finding algorithms, one is sparsity-aware algorithm and the other is weighted quantile sketch with proofs. Compared to other popular systems, XGBoost gives outstanding performance in different experiment settings, which explains why it is considered the best and fastest algorithm for structured or tabular data.<br />
<br />
= Source =</div>Z66songhttp://wiki.math.uwaterloo.ca/statwiki/index.php?title=summary&diff=34450summary2018-03-16T20:24:58Z<p>Z66song: </p>
<hr />
<div>'''''XGBoost: A Scalable Tree Boosting System'''''<br />
<br />
== Presented by ==<br />
<br />
Jiang, Cong<br />
<br />
Song, Ziwei<br />
<br />
Ye, Zhaoshan<br />
<br />
Zhang, Wenling<br />
<br />
= Introduction =<br />
<br />
In recent years, XGBoost has been the most popular tool in Kaggle, as well as other machine learning and data mining challenges. The main reason to its success is that this tree boosting system is highly scalable , which means it could be used to solve a variety of different problems. <br />
<br />
This paper is written by XGBoost's father, Tianqi Chen. It gives us an overview on how he used algorithmic optimizations as well as some important systems to develop XGBoost. He explained it in the following manner:<br />
<br />
1. Based on gradient tree boosting algorithm as well as using shrinkage and feature subsampling to prevent overfitting, Chen introduced an end-to-end tree boosting system.<br />
<br />
2.<br />
<br />
= Tree Boosting=<br />
<br />
<br />
=== Regularized Learning Objective ===<br />
Suppose the given data set has m features and n examples. A tree ensemble model uses K additive functions to predict the output.<br />
[[File:eq1.png|center|eq1.png]]<br />
where F is the space of regression trees(CART).<br />
The author suggests the following regularized objective function to make improvements comparing to gradient boosting<br />
[[File:eq2.png|center|eq2.png]]<br />
Regression tree contains continuous scores on each leaf, which is different with decision trees we have learned in class. The following 2 figures represent how to use the decision rules in the trees to classify it into the leaves and calculate the final prediction by summing up the score in the corresponding leaves.<br />
<br />
[[File:Picture2.1.png|500px| thumb| center| Figure 0: Prediction Score of Tree Structure.]]<br />
<br />
<br />
[[File:Picture2.2.png|500px|thumb|center|Figure 1: Tree Ensemble Model. The final prediction for a given example is the sum of predictions from each tree. ]]<br />
<br />
=== Gradient Tree Boosting ===<br />
The tree ensemble model (Equation 2.) is trained in an additive manner. We need to add a new <math>f</math> to our objective.<br />
[[File:eq2.2.png|center|eq2.2.png]]<br />
After applying the second-order Taylor expansion to the above formula and removing constant terms, we can finally have<br />
[[File:eq3.png|center|eq3.png]]<br />
Define <math>I_{j} = \{i| q(X_{i}) \}</math><br />
As the instance set of leaf j. We can rewrite Equation (3) as<br />
[[File:eq4.png|center|eq4.png]]<br />
By computing the optimal weight and corresponding optimal value, we can use following equation as the scoring function which can measure the quality of a tree structure. Lower score means better tree structure.<br />
[[File:eq5.png|center|eq5.png]]<br />
The following figure shows an example how this score is calculated.<br />
[[File:Picture2.3.png|500px| thumb|center|Figure 2: Structure Score Calculation. We only need to sum up the gradient and second order gradient statistics on each leaf, then apply the scoring formula to get the quality score. ]]<br />
<br />
<br />
<br />
Normally it is impossible to enumerate all the possible tree structures. Assume <math>I_{L}</math> and <math>I_{R}</math> are the instance sets of left and right nodes after the split, and <math> I = I_{L} + I_{R}</math>. Then the loss reduction after the split is given by<br />
[[File:eq6.png|center]]<br />
<br />
=== Shrinkage and Column Subsampling ===<br />
In addition to the regularized objective, the paper also introduces two techniques to prevent over-fitting. <br />
<br />
Method 1: Shrinkage. After each step of tree boosting, we scale newly added weights by a shrinkage factor η, which reduces the influence of every individual tree and leaves the space for future trees. <br />
<br />
Method 2: Column (feature) subsampling. We try to build each individual tree with only a portion of predictors chosen randomly. This procedure encourages the variance between the trees (i.e. no same tree occurs). Implemented in an opensource packages, this method prevents over-fitting and speeds up the computations of parallel algorithm as well.<br />
<br />
= Split Finding Algorithms =<br />
The ways to find the best split can be realized by two main methods. <br />
<br />
Basic exact greedy algorithm: enumerate over all the possible splits on all the features, which is computationally demanding. Sort the data first, and visit the data to accumulate the gradient statistics. <br />
<br />
Approximate algorithm: To achieve an effective gradient tree boosting model, we first propose candidate splitting points by percentiles of feature distribution, map the continuous features into buckets split and then finds the best solution among the proposals by the aggregated statistics. <br />
The algorithm can be varied by the time the proposal is given. <br />
The global variant: all candidate splits are proposed at the initial phase of tree construction and used at all levels. <br />
The local variant: the candidate splits are refined after each splits. It is potential for deeper trees and it requires fewer candidates. <br />
Methods can be chosen according to users’ needs.<br />
[[File:exact_greedy.png|500px| thumb|center ]]<br />
[[File:appro_greedy.png|500px| thumb|center]]<br />
<br />
<br />
= Conclusion=<br />
<br />
<br />
<br />
= Source =<br />
<br />
===**Sample format===<br />
[http://www.wildml.com/2015/09/recurrent-neural-networks-tutorial-part-1-introduction-to-rnns/ Recurrent neural networks] are a variation of deep neural networks that are capable of storing information about previous hidden states in special memory layers.<ref name=lstm><br />
Hochreiter, Sepp, and Jürgen Schmidhuber. [http://deeplearning.cs.cmu.edu/pdfs/Hochreiter97_lstm.pdf "Long short-term memory."] Neural computation 9.8 (1997): 1735-1780.<br />
</ref> Unlike feed forward neural networks that take in a single fixed length vector input and output a fixed length vector output, recurrent neural networks can take in a sequence of fixed length vectors as input, because of their ability to store information and maintain a connection between inputs through this memory layer. By comparison, previous inputs would have no impact on current output for feed forward neural networks, whereas they can impact current input in a recurrent neural network. (This paper used the LSTM formulation from Graves<ref name=grave><br />
Graves, Alex. [http://arxiv.org/pdf/1308.0850.pdf "Generating sequences with recurrent neural networks."] arXiv preprint arXiv:1308.0850 (2013).<br />
</ref>)<br />
<br />
<br />
Where <math>\,S</math> is the base/source sentence, <math>\,T</math> is the paired translated sentence and <math>\,T_r</math> is the total training set. This objective function is to maximize the log probability of a correct translation <math>\,T</math> given the base/source sentence <math>\,S</math> over the entire training set. Once the training is complete, translations are produced by finding the most likely translation according to LSTM:<br />
<br />
:<math>\hat{T} = \underset{T}{\operatorname{arg\ max}}\ p(T|S)</math><br />
<br />
<br />It has been showed that Long Short-Term Memory recurrent neural networks have the ability to generate both discrete and real-valued sequences with complex, long-range structure using next-step prediction <ref name=grave><br />
Reference<br />
</ref>.</div>Z66songhttp://wiki.math.uwaterloo.ca/statwiki/index.php?title=summary&diff=34449summary2018-03-16T20:22:42Z<p>Z66song: /* Split Finding Algorithms */</p>
<hr />
<div>'''''XGBoost: A Scalable Tree Boosting System'''''<br />
<br />
== Presented by ==<br />
<br />
Jiang, Cong<br />
<br />
Song, Ziwei<br />
<br />
Ye, Zhaoshan<br />
<br />
Zhang, Wenling<br />
<br />
= Introduction =<br />
<br />
In recent years, XGBoost has been the most popular tool in Kaggle, as well as other machine learning and data mining challenges. The main reason to its success is that this tree boosting system is highly scalable , which means it could be used to solve a variety of different problems. <br />
<br />
This paper is written by XGBoost's father, Tianqi Chen. It gives us an overview on how he used algorithmic optimizations as well as some important systems to develop XGBoost. He explained it in the following manner:<br />
<br />
1. Based on gradient tree boosting algorithm as well as using shrinkage and feature subsampling to prevent overfitting, Chen introduced an end-to-end tree boosting system.<br />
<br />
2.<br />
<br />
= Tree Boosting=<br />
<br />
<br />
=== Regularized Learning Objective ===<br />
Suppose the given data set has m features and n examples. A tree ensemble model uses K additive functions to predict the output.<br />
[[File:eq1.png|center|eq1.png]]<br />
where F is the space of regression trees(CART).<br />
The author suggests the following regularized objective function to make improvements comparing to gradient boosting<br />
[[File:eq2.png|center|eq2.png]]<br />
Regression tree contains continuous scores on each leaf, which is different with decision trees we have learned in class. The following 2 figures represent how to use the decision rules in the trees to classify it into the leaves and calculate the final prediction by summing up the score in the corresponding leaves.<br />
<br />
[[File:Picture2.1.png|500px| thumb| center| Figure 0: Prediction Score of Tree Structure.]]<br />
<br />
<br />
[[File:Picture2.2.png|500px|thumb|center|Figure 1: Tree Ensemble Model. The final prediction for a given example is the sum of predictions from each tree. ]]<br />
<br />
=== Gradient Tree Boosting ===<br />
The tree ensemble model (Equation 2.) is trained in an additive manner. We need to add a new <math>f</math> to our objective.<br />
[[File:eq2.2.png|center|eq2.2.png]]<br />
After applying the second-order Taylor expansion to the above formula and removing constant terms, we can finally have<br />
[[File:eq3.png|center|eq3.png]]<br />
Define <math>I_{j} = \{i| q(X_{i}) \}</math><br />
As the instance set of leaf j. We can rewrite Equation (3) as<br />
[[File:eq4.png|center|eq4.png]]<br />
By computing the optimal weight and corresponding optimal value, we can use following equation as the scoring function which can measure the quality of a tree structure. Lower score means better tree structure.<br />
[[File:eq5.png|center|eq5.png]]<br />
The following figure shows an example how this score is calculated.<br />
[[File:Picture2.3.png|500px| thumb|center|Figure 2: Structure Score Calculation. We only need to sum up the gradient and second order gradient statistics on each leaf, then apply the scoring formula to get the quality score. ]]<br />
<br />
<br />
<br />
Normally it is impossible to enumerate all the possible tree structures. Assume <math>I_{L}</math> and <math>I_{R}</math> are the instance sets of left and right nodes after the split, and <math> I = I_{L} + I_{R}</math>. Then the loss reduction after the split is given by<br />
[[File:eq6.png|center]]<br />
<br />
=== Shrinkage and Column Subsampling ===<br />
In addition to the regularized objective, the paper also introduces two techniques to prevent over-fitting. <br />
<br />
Method 1: Shrinkage. After each step of tree boosting, we scale newly added weights by a shrinkage factor η, which reduces the influence of every individual tree and leaves the space for future trees. <br />
<br />
Method 2: Column (feature) subsampling. We try to build each individual tree with only a portion of predictors chosen randomly. This procedure encourages the variance between the trees (i.e. no same tree occurs). Implemented in an opensource packages, this method prevents over-fitting and speeds up the computations of parallel algorithm as well.<br />
<br />
= Split Finding Algorithms =<br />
The ways to find the best split can be realized by two main methods. <br />
<br />
Basic exact greedy algorithm: enumerate over all the possible splits on all the features, which is computationally demanding. Sort the data first, and visit the data to accumulate the gradient statistics. <br />
<br />
Approximate algorithm: To achieve an effective gradient tree boosting model, we first propose candidate splitting points by percentiles of feature distribution, map the continuous features into buckets split and then finds the best solution among the proposals by the aggregated statistics. <br />
The algorithm can be varied by the time the proposal is given. <br />
The global variant: all candidate splits are proposed at the initial phase of tree construction and used at all levels. <br />
The local variant: the candidate splits are refined after each splits. It is potential for deeper trees and it requires fewer candidates. <br />
Methods can be chosen according to users’ needs.<br />
[[File:exact_greedy.png|500px| thumb|center ]]<br />
[[File:appro_greedy.png|500px| thumb|center]]<br />
<br />
= System Design =<br />
<br />
=Related Works=<br />
<br />
= End to End Evaluations =<br />
<br />
<br />
=== Results ===<br />
<br />
<br />
<br />
<br />
= Conclusion=<br />
<br />
<br />
<br />
<br />
<br />
<br />
= Criticisms =<br />
<br />
= Source =<br />
<br />
===**Sample format===<br />
[http://www.wildml.com/2015/09/recurrent-neural-networks-tutorial-part-1-introduction-to-rnns/ Recurrent neural networks] are a variation of deep neural networks that are capable of storing information about previous hidden states in special memory layers.<ref name=lstm><br />
Hochreiter, Sepp, and Jürgen Schmidhuber. [http://deeplearning.cs.cmu.edu/pdfs/Hochreiter97_lstm.pdf "Long short-term memory."] Neural computation 9.8 (1997): 1735-1780.<br />
</ref> Unlike feed forward neural networks that take in a single fixed length vector input and output a fixed length vector output, recurrent neural networks can take in a sequence of fixed length vectors as input, because of their ability to store information and maintain a connection between inputs through this memory layer. By comparison, previous inputs would have no impact on current output for feed forward neural networks, whereas they can impact current input in a recurrent neural network. (This paper used the LSTM formulation from Graves<ref name=grave><br />
Graves, Alex. [http://arxiv.org/pdf/1308.0850.pdf "Generating sequences with recurrent neural networks."] arXiv preprint arXiv:1308.0850 (2013).<br />
</ref>)<br />
<br />
<br />
Where <math>\,S</math> is the base/source sentence, <math>\,T</math> is the paired translated sentence and <math>\,T_r</math> is the total training set. This objective function is to maximize the log probability of a correct translation <math>\,T</math> given the base/source sentence <math>\,S</math> over the entire training set. Once the training is complete, translations are produced by finding the most likely translation according to LSTM:<br />
<br />
:<math>\hat{T} = \underset{T}{\operatorname{arg\ max}}\ p(T|S)</math><br />
<br />
<br />It has been showed that Long Short-Term Memory recurrent neural networks have the ability to generate both discrete and real-valued sequences with complex, long-range structure using next-step prediction <ref name=grave><br />
Reference<br />
</ref>.</div>Z66songhttp://wiki.math.uwaterloo.ca/statwiki/index.php?title=summary&diff=34448summary2018-03-16T20:22:10Z<p>Z66song: /* Split Finding Algorithms */</p>
<hr />
<div>'''''XGBoost: A Scalable Tree Boosting System'''''<br />
<br />
== Presented by ==<br />
<br />
Jiang, Cong<br />
<br />
Song, Ziwei<br />
<br />
Ye, Zhaoshan<br />
<br />
Zhang, Wenling<br />
<br />
= Introduction =<br />
<br />
In recent years, XGBoost has been the most popular tool in Kaggle, as well as other machine learning and data mining challenges. The main reason to its success is that this tree boosting system is highly scalable , which means it could be used to solve a variety of different problems. <br />
<br />
This paper is written by XGBoost's father, Tianqi Chen. It gives us an overview on how he used algorithmic optimizations as well as some important systems to develop XGBoost. He explained it in the following manner:<br />
<br />
1. Based on gradient tree boosting algorithm as well as using shrinkage and feature subsampling to prevent overfitting, Chen introduced an end-to-end tree boosting system.<br />
<br />
2.<br />
<br />
= Tree Boosting=<br />
<br />
<br />
=== Regularized Learning Objective ===<br />
Suppose the given data set has m features and n examples. A tree ensemble model uses K additive functions to predict the output.<br />
[[File:eq1.png|center|eq1.png]]<br />
where F is the space of regression trees(CART).<br />
The author suggests the following regularized objective function to make improvements comparing to gradient boosting<br />
[[File:eq2.png|center|eq2.png]]<br />
Regression tree contains continuous scores on each leaf, which is different with decision trees we have learned in class. The following 2 figures represent how to use the decision rules in the trees to classify it into the leaves and calculate the final prediction by summing up the score in the corresponding leaves.<br />
<br />
[[File:Picture2.1.png|500px| thumb| center| Figure 0: Prediction Score of Tree Structure.]]<br />
<br />
<br />
[[File:Picture2.2.png|500px|thumb|center|Figure 1: Tree Ensemble Model. The final prediction for a given example is the sum of predictions from each tree. ]]<br />
<br />
=== Gradient Tree Boosting ===<br />
The tree ensemble model (Equation 2.) is trained in an additive manner. We need to add a new <math>f</math> to our objective.<br />
[[File:eq2.2.png|center|eq2.2.png]]<br />
After applying the second-order Taylor expansion to the above formula and removing constant terms, we can finally have<br />
[[File:eq3.png|center|eq3.png]]<br />
Define <math>I_{j} = \{i| q(X_{i}) \}</math><br />
As the instance set of leaf j. We can rewrite Equation (3) as<br />
[[File:eq4.png|center|eq4.png]]<br />
By computing the optimal weight and corresponding optimal value, we can use following equation as the scoring function which can measure the quality of a tree structure. Lower score means better tree structure.<br />
[[File:eq5.png|center|eq5.png]]<br />
The following figure shows an example how this score is calculated.<br />
[[File:Picture2.3.png|500px| thumb|center|Figure 2: Structure Score Calculation. We only need to sum up the gradient and second order gradient statistics on each leaf, then apply the scoring formula to get the quality score. ]]<br />
<br />
<br />
<br />
Normally it is impossible to enumerate all the possible tree structures. Assume <math>I_{L}</math> and <math>I_{R}</math> are the instance sets of left and right nodes after the split, and <math> I = I_{L} + I_{R}</math>. Then the loss reduction after the split is given by<br />
[[File:eq6.png|center]]<br />
<br />
=== Shrinkage and Column Subsampling ===<br />
In addition to the regularized objective, the paper also introduces two techniques to prevent over-fitting. <br />
<br />
Method 1: Shrinkage. After each step of tree boosting, we scale newly added weights by a shrinkage factor η, which reduces the influence of every individual tree and leaves the space for future trees. <br />
<br />
Method 2: Column (feature) subsampling. We try to build each individual tree with only a portion of predictors chosen randomly. This procedure encourages the variance between the trees (i.e. no same tree occurs). Implemented in an opensource packages, this method prevents over-fitting and speeds up the computations of parallel algorithm as well.<br />
<br />
= Split Finding Algorithms =<br />
The ways to find the best split can be realized by two main methods. <br />
<br />
Basic exact greedy algorithm: enumerate over all the possible splits on all the features, which is computationally demanding. Sort the data first, and visit the data to accumulate the gradient statistics. <br />
<br />
Approximate algorithm: To achieve an effective gradient tree boosting model, we first propose candidate splitting points by percentiles of feature distribution, map the continuous features into buckets split and then finds the best solution among the proposals by the aggregated statistics. <br />
The algorithm can be varied by the time the proposal is given. <br />
The global variant: all candidate splits are proposed at the initial phase of tree construction and used at all levels. <br />
The local variant: the candidate splits are refined after each splits. It is potential for deeper trees and it requires fewer candidates. <br />
Methods can be chosen according to users’ needs.<br />
[[File:exact_greedy.png|500px| thumb|left ]]<br />
[[File:appro_greedy.png|500px| thumb|left ]]<br />
<br />
= System Design =<br />
<br />
=Related Works=<br />
<br />
= End to End Evaluations =<br />
<br />
<br />
=== Results ===<br />
<br />
<br />
<br />
<br />
= Conclusion=<br />
<br />
<br />
<br />
<br />
<br />
<br />
= Criticisms =<br />
<br />
= Source =<br />
<br />
===**Sample format===<br />
[http://www.wildml.com/2015/09/recurrent-neural-networks-tutorial-part-1-introduction-to-rnns/ Recurrent neural networks] are a variation of deep neural networks that are capable of storing information about previous hidden states in special memory layers.<ref name=lstm><br />
Hochreiter, Sepp, and Jürgen Schmidhuber. [http://deeplearning.cs.cmu.edu/pdfs/Hochreiter97_lstm.pdf "Long short-term memory."] Neural computation 9.8 (1997): 1735-1780.<br />
</ref> Unlike feed forward neural networks that take in a single fixed length vector input and output a fixed length vector output, recurrent neural networks can take in a sequence of fixed length vectors as input, because of their ability to store information and maintain a connection between inputs through this memory layer. By comparison, previous inputs would have no impact on current output for feed forward neural networks, whereas they can impact current input in a recurrent neural network. (This paper used the LSTM formulation from Graves<ref name=grave><br />
Graves, Alex. [http://arxiv.org/pdf/1308.0850.pdf "Generating sequences with recurrent neural networks."] arXiv preprint arXiv:1308.0850 (2013).<br />
</ref>)<br />
<br />
<br />
Where <math>\,S</math> is the base/source sentence, <math>\,T</math> is the paired translated sentence and <math>\,T_r</math> is the total training set. This objective function is to maximize the log probability of a correct translation <math>\,T</math> given the base/source sentence <math>\,S</math> over the entire training set. Once the training is complete, translations are produced by finding the most likely translation according to LSTM:<br />
<br />
:<math>\hat{T} = \underset{T}{\operatorname{arg\ max}}\ p(T|S)</math><br />
<br />
<br />It has been showed that Long Short-Term Memory recurrent neural networks have the ability to generate both discrete and real-valued sequences with complex, long-range structure using next-step prediction <ref name=grave><br />
Reference<br />
</ref>.</div>Z66songhttp://wiki.math.uwaterloo.ca/statwiki/index.php?title=File:appro_greedy.png&diff=34447File:appro greedy.png2018-03-16T20:19:22Z<p>Z66song: </p>
<hr />
<div></div>Z66songhttp://wiki.math.uwaterloo.ca/statwiki/index.php?title=File:exact_greedy.png&diff=34446File:exact greedy.png2018-03-16T20:18:56Z<p>Z66song: </p>
<hr />
<div></div>Z66songhttp://wiki.math.uwaterloo.ca/statwiki/index.php?title=summary&diff=34445summary2018-03-16T20:17:27Z<p>Z66song: /* Split Finding Algorithms */</p>
<hr />
<div>'''''XGBoost: A Scalable Tree Boosting System'''''<br />
<br />
== Presented by ==<br />
<br />
Jiang, Cong<br />
<br />
Song, Ziwei<br />
<br />
Ye, Zhaoshan<br />
<br />
Zhang, Wenling<br />
<br />
= Introduction =<br />
<br />
In recent years, XGBoost has been the most popular tool in Kaggle, as well as other machine learning and data mining challenges. The main reason to its success is that this tree boosting system is highly scalable , which means it could be used to solve a variety of different problems. <br />
<br />
This paper is written by XGBoost's father, Tianqi Chen. It gives us an overview on how he used algorithmic optimizations as well as some important systems to develop XGBoost. He explained it in the following manner:<br />
<br />
1. Based on gradient tree boosting algorithm as well as using shrinkage and feature subsampling to prevent overfitting, Chen introduced an end-to-end tree boosting system.<br />
<br />
2.<br />
<br />
= Tree Boosting=<br />
<br />
<br />
=== Regularized Learning Objective ===<br />
Suppose the given data set has m features and n examples. A tree ensemble model uses K additive functions to predict the output.<br />
[[File:eq1.png|center|eq1.png]]<br />
where F is the space of regression trees(CART).<br />
The author suggests the following regularized objective function to make improvements comparing to gradient boosting<br />
[[File:eq2.png|center|eq2.png]]<br />
Regression tree contains continuous scores on each leaf, which is different with decision trees we have learned in class. The following 2 figures represent how to use the decision rules in the trees to classify it into the leaves and calculate the final prediction by summing up the score in the corresponding leaves.<br />
<br />
[[File:Picture2.1.png|500px| thumb| center| Figure 0: Prediction Score of Tree Structure.]]<br />
<br />
<br />
[[File:Picture2.2.png|500px|thumb|center|Figure 1: Tree Ensemble Model. The final prediction for a given example is the sum of predictions from each tree. ]]<br />
<br />
=== Gradient Tree Boosting ===<br />
The tree ensemble model (Equation 2.) is trained in an additive manner. We need to add a new <math>f</math> to our objective.<br />
[[File:eq2.2.png|center|eq2.2.png]]<br />
After applying the second-order Taylor expansion to the above formula and removing constant terms, we can finally have<br />
[[File:eq3.png|center|eq3.png]]<br />
Define <math>I_{j} = \{i| q(X_{i}) \}</math><br />
As the instance set of leaf j. We can rewrite Equation (3) as<br />
[[File:eq4.png|center|eq4.png]]<br />
By computing the optimal weight and corresponding optimal value, we can use following equation as the scoring function which can measure the quality of a tree structure. Lower score means better tree structure.<br />
[[File:eq5.png|center|eq5.png]]<br />
The following figure shows an example how this score is calculated.<br />
[[File:Picture2.3.png|500px| thumb|center|Figure 2: Structure Score Calculation. We only need to sum up the gradient and second order gradient statistics on each leaf, then apply the scoring formula to get the quality score. ]]<br />
<br />
<br />
<br />
Normally it is impossible to enumerate all the possible tree structures. Assume <math>I_{L}</math> and <math>I_{R}</math> are the instance sets of left and right nodes after the split, and <math> I = I_{L} + I_{R}</math>. Then the loss reduction after the split is given by<br />
[[File:eq6.png|center]]<br />
<br />
=== Shrinkage and Column Subsampling ===<br />
In addition to the regularized objective, the paper also introduces two techniques to prevent over-fitting. <br />
<br />
Method 1: Shrinkage. After each step of tree boosting, we scale newly added weights by a shrinkage factor η, which reduces the influence of every individual tree and leaves the space for future trees. <br />
<br />
Method 2: Column (feature) subsampling. We try to build each individual tree with only a portion of predictors chosen randomly. This procedure encourages the variance between the trees (i.e. no same tree occurs). Implemented in an opensource packages, this method prevents over-fitting and speeds up the computations of parallel algorithm as well.<br />
<br />
= Split Finding Algorithms =<br />
The ways to find the best split can be realized by two main methods. <br />
<br />
Basic exact greedy algorithm: enumerate over all the possible splits on all the features, which is computationally demanding. Sort the data first, and visit the data to accumulate the gradient statistics. <br />
<br />
Approximate algorithm: To achieve an effective gradient tree boosting model, we first propose candidate splitting points by percentiles of feature distribution, map the continuous features into buckets split and then finds the best solution among the proposals by the aggregated statistics. <br />
The algorithm can be varied by the time the proposal is given. <br />
The global variant: all candidate splits are proposed at the initial phase of tree construction and used at all levels. <br />
The local variant: the candidate splits are refined after each splits. It is potential for deeper trees and it requires fewer candidates. <br />
Methods can be chosen according to users’ needs.<br />
<br />
= System Design =<br />
<br />
=Related Works=<br />
<br />
= End to End Evaluations =<br />
<br />
<br />
=== Results ===<br />
<br />
<br />
<br />
<br />
= Conclusion=<br />
<br />
<br />
<br />
<br />
<br />
<br />
= Criticisms =<br />
<br />
= Source =<br />
<br />
===**Sample format===<br />
[http://www.wildml.com/2015/09/recurrent-neural-networks-tutorial-part-1-introduction-to-rnns/ Recurrent neural networks] are a variation of deep neural networks that are capable of storing information about previous hidden states in special memory layers.<ref name=lstm><br />
Hochreiter, Sepp, and Jürgen Schmidhuber. [http://deeplearning.cs.cmu.edu/pdfs/Hochreiter97_lstm.pdf "Long short-term memory."] Neural computation 9.8 (1997): 1735-1780.<br />
</ref> Unlike feed forward neural networks that take in a single fixed length vector input and output a fixed length vector output, recurrent neural networks can take in a sequence of fixed length vectors as input, because of their ability to store information and maintain a connection between inputs through this memory layer. By comparison, previous inputs would have no impact on current output for feed forward neural networks, whereas they can impact current input in a recurrent neural network. (This paper used the LSTM formulation from Graves<ref name=grave><br />
Graves, Alex. [http://arxiv.org/pdf/1308.0850.pdf "Generating sequences with recurrent neural networks."] arXiv preprint arXiv:1308.0850 (2013).<br />
</ref>)<br />
<br />
<br />
Where <math>\,S</math> is the base/source sentence, <math>\,T</math> is the paired translated sentence and <math>\,T_r</math> is the total training set. This objective function is to maximize the log probability of a correct translation <math>\,T</math> given the base/source sentence <math>\,S</math> over the entire training set. Once the training is complete, translations are produced by finding the most likely translation according to LSTM:<br />
<br />
:<math>\hat{T} = \underset{T}{\operatorname{arg\ max}}\ p(T|S)</math><br />
<br />
<br />It has been showed that Long Short-Term Memory recurrent neural networks have the ability to generate both discrete and real-valued sequences with complex, long-range structure using next-step prediction <ref name=grave><br />
Reference<br />
</ref>.</div>Z66songhttp://wiki.math.uwaterloo.ca/statwiki/index.php?title=summary&diff=34444summary2018-03-16T20:16:29Z<p>Z66song: /* Tree Boosting */</p>
<hr />
<div>'''''XGBoost: A Scalable Tree Boosting System'''''<br />
<br />
== Presented by ==<br />
<br />
Jiang, Cong<br />
<br />
Song, Ziwei<br />
<br />
Ye, Zhaoshan<br />
<br />
Zhang, Wenling<br />
<br />
= Introduction =<br />
<br />
In recent years, XGBoost has been the most popular tool in Kaggle, as well as other machine learning and data mining challenges. The main reason to its success is that this tree boosting system is highly scalable , which means it could be used to solve a variety of different problems. <br />
<br />
This paper is written by XGBoost's father, Tianqi Chen. It gives us an overview on how he used algorithmic optimizations as well as some important systems to develop XGBoost. He explained it in the following manner:<br />
<br />
1. Based on gradient tree boosting algorithm as well as using shrinkage and feature subsampling to prevent overfitting, Chen introduced an end-to-end tree boosting system.<br />
<br />
2.<br />
<br />
= Tree Boosting=<br />
<br />
<br />
=== Regularized Learning Objective ===<br />
Suppose the given data set has m features and n examples. A tree ensemble model uses K additive functions to predict the output.<br />
[[File:eq1.png|center|eq1.png]]<br />
where F is the space of regression trees(CART).<br />
The author suggests the following regularized objective function to make improvements comparing to gradient boosting<br />
[[File:eq2.png|center|eq2.png]]<br />
Regression tree contains continuous scores on each leaf, which is different with decision trees we have learned in class. The following 2 figures represent how to use the decision rules in the trees to classify it into the leaves and calculate the final prediction by summing up the score in the corresponding leaves.<br />
<br />
[[File:Picture2.1.png|500px| thumb| center| Figure 0: Prediction Score of Tree Structure.]]<br />
<br />
<br />
[[File:Picture2.2.png|500px|thumb|center|Figure 1: Tree Ensemble Model. The final prediction for a given example is the sum of predictions from each tree. ]]<br />
<br />
=== Gradient Tree Boosting ===<br />
The tree ensemble model (Equation 2.) is trained in an additive manner. We need to add a new <math>f</math> to our objective.<br />
[[File:eq2.2.png|center|eq2.2.png]]<br />
After applying the second-order Taylor expansion to the above formula and removing constant terms, we can finally have<br />
[[File:eq3.png|center|eq3.png]]<br />
Define <math>I_{j} = \{i| q(X_{i}) \}</math><br />
As the instance set of leaf j. We can rewrite Equation (3) as<br />
[[File:eq4.png|center|eq4.png]]<br />
By computing the optimal weight and corresponding optimal value, we can use following equation as the scoring function which can measure the quality of a tree structure. Lower score means better tree structure.<br />
[[File:eq5.png|center|eq5.png]]<br />
The following figure shows an example how this score is calculated.<br />
[[File:Picture2.3.png|500px| thumb|center|Figure 2: Structure Score Calculation. We only need to sum up the gradient and second order gradient statistics on each leaf, then apply the scoring formula to get the quality score. ]]<br />
<br />
<br />
<br />
Normally it is impossible to enumerate all the possible tree structures. Assume <math>I_{L}</math> and <math>I_{R}</math> are the instance sets of left and right nodes after the split, and <math> I = I_{L} + I_{R}</math>. Then the loss reduction after the split is given by<br />
[[File:eq6.png|center]]<br />
<br />
=== Shrinkage and Column Subsampling ===<br />
In addition to the regularized objective, the paper also introduces two techniques to prevent over-fitting. <br />
<br />
Method 1: Shrinkage. After each step of tree boosting, we scale newly added weights by a shrinkage factor η, which reduces the influence of every individual tree and leaves the space for future trees. <br />
<br />
Method 2: Column (feature) subsampling. We try to build each individual tree with only a portion of predictors chosen randomly. This procedure encourages the variance between the trees (i.e. no same tree occurs). Implemented in an opensource packages, this method prevents over-fitting and speeds up the computations of parallel algorithm as well.<br />
<br />
= Split Finding Algorithms =<br />
<br />
= System Design =<br />
<br />
=Related Works=<br />
<br />
= End to End Evaluations =<br />
<br />
<br />
=== Results ===<br />
<br />
<br />
<br />
<br />
= Conclusion=<br />
<br />
<br />
<br />
<br />
<br />
<br />
= Criticisms =<br />
<br />
= Source =<br />
<br />
===**Sample format===<br />
[http://www.wildml.com/2015/09/recurrent-neural-networks-tutorial-part-1-introduction-to-rnns/ Recurrent neural networks] are a variation of deep neural networks that are capable of storing information about previous hidden states in special memory layers.<ref name=lstm><br />
Hochreiter, Sepp, and Jürgen Schmidhuber. [http://deeplearning.cs.cmu.edu/pdfs/Hochreiter97_lstm.pdf "Long short-term memory."] Neural computation 9.8 (1997): 1735-1780.<br />
</ref> Unlike feed forward neural networks that take in a single fixed length vector input and output a fixed length vector output, recurrent neural networks can take in a sequence of fixed length vectors as input, because of their ability to store information and maintain a connection between inputs through this memory layer. By comparison, previous inputs would have no impact on current output for feed forward neural networks, whereas they can impact current input in a recurrent neural network. (This paper used the LSTM formulation from Graves<ref name=grave><br />
Graves, Alex. [http://arxiv.org/pdf/1308.0850.pdf "Generating sequences with recurrent neural networks."] arXiv preprint arXiv:1308.0850 (2013).<br />
</ref>)<br />
<br />
<br />
Where <math>\,S</math> is the base/source sentence, <math>\,T</math> is the paired translated sentence and <math>\,T_r</math> is the total training set. This objective function is to maximize the log probability of a correct translation <math>\,T</math> given the base/source sentence <math>\,S</math> over the entire training set. Once the training is complete, translations are produced by finding the most likely translation according to LSTM:<br />
<br />
:<math>\hat{T} = \underset{T}{\operatorname{arg\ max}}\ p(T|S)</math><br />
<br />
<br />It has been showed that Long Short-Term Memory recurrent neural networks have the ability to generate both discrete and real-valued sequences with complex, long-range structure using next-step prediction <ref name=grave><br />
Reference<br />
</ref>.</div>Z66song