Revision as of 03:39, 18 October 2017

Modeling

VOC Classification

Overview

The VOC device is designed to judge whether the tobacco is heathy or gets infected. Since this is an inquiry experiment, algorithms in data analysis are widely use in our modeling. We do data preprocessing, data analysis, and algorithm optimization on the data collected by VOC device. Finally, we use Logistic regression and detect the infected tobacco with 91% confidence.

Data preprocessing

First we defragment the raw input data, and reorganize them into a matrix. 10 VOC factors are served as features, and the status(heathy or infected) is served as tag to be predicted.

Then we analysis the data using box plot and discover that most data are normal, but some records are singular, whose box plot are show as folowing:

We remove those records with singular value, and the data left obey normal distribution:

Data analysis

Our target is to create a model and predict tobacco's status according to 10 input features. This is a classic two classification problem, and there are several algrithm to solve it. The sampling algorithm is cross validation and the scoring policy we apply is ridit test.

Decision Tree

First we use decision tree based on information theory. ID3 decision tree is used to reduce the most information gain, and CART tree is used to reduce the GINI index. The performance of these two algorithm is almost the same.

$R=0.83$

MLP

The second algorithm we apply is Multi-Layer Perception, also called neutral network. In this model, we use more than 100 neurons in each layer and the activation function is relu.

The result of MLP is much better than decision tree.

$R=0.89$

Leaner Model

Although the performance of MLP has been good enough, it's difficult to extract konwledge learn by algorithm, the interpretability is weak. Why don't we try a simple model with high interpretability? First we try LDA algorithm to compress the 10dimensions data into 2 dimensions.

$J=\frac{||w^T\mu_0-w^T\mu_1||^2}{w^T\Sigma_0w+w^T\Sigma_1w}=\frac{w^T(\mu_0-\mu_1)(\mu_0-\mu_1)^Tw}{w^T\Sigma_0w+w^T\Sigma_1w}$

We define $S_w$ as within-class scatter matrix

$S_w=\Sigma_0+\Sigma_1=\sum_{x\in X_0}(x-\mu_0)(x-\mu_0)^T+\sum_{x\in X_1}(x-\mu_1)(x-\mu_1)^T$

We define $S_b$ as between-class scatter matrix

$S_b=(\mu_0-\mu_1)(\mu_0-\mu_1)^T$

So

$J=\frac{w^TS_bw}{w^TS_ww}$ , the target of LDA is maxmize the $J$

The result of LDA algorithm is as following and $R=0.89$ :

This result prove the data are linear separable, then we choose logistics regression algorithm.

We difine $logitP=ln\frac{y}{1-y}\in (-\infty,+\infty)$

$p(y=1|x)=\frac{e^{w^Tx+b}}{1+e^{w^Tx+b}}$

$p(y=1|x)=\frac{1}{1+e^{w^Tx+b}}$

$l(w,b)=\sum_{i=1}^{m}lnp(y_i|x_i;w,b)$

Then we can apply maximum likelihood method algorithm to estimate the paramaters.

The result is as following:

Weight:

[[ 0.1819504 0.38788225 0.01350023 0.39594948 0.17799418

0.42087034

-0.57733395 -0.23876003 -0.00532918 -0.46174515]]

Intercept:

[ 0.00937812]

Effect:

D 35.300735

B 22.596339

F 18.289277

E 10.265025

C 0.393225

I -1.575564

A -10.679026

H -14.398440

G -26.211964

J -39.130542

dtype: float64

Score:

0.894333333333

Algorithm optimization

From the result of logistics regression, factor C and I and etc. are with less important weight, these factors maybe disturb the classifaction. We try to reduce unimportant factors and simplify the model.

Finally, we reserve 4 factors with which we can predict the tobacco in 91% confidence and also reduce the VOC device.

Weight:

[[ 0.53196697 0.3404023 -0.53555988 -0.45588715]]

Intercept:

[-0.01204088]

Effect:

D 33.217011

F 15.492680

G -17.319760

J -33.967849