Machine Learning - Logistic Regression

Logistic regression aims to solve classification problems. It does this by predicting categorical outcomes, unlike linear regression that predicts a continuous outcome. In the simplest case there are two outcomes, which is called binomial, an example of which is predicting if a tumor is malignant or benign. Other cases have more than two outcomes to classify, in this case it is called multinomial. A common example for multinomial logistic regression would be predicting the class of an iris flower between 3 different species. Here we will be using basic logistic regression to predict a binomial variable. This means it has only two possible outcomes.

In Python we have modules that will do the work for us. Start by importing the NumPy module. import numpy Store the independent variables in X. Store the dependent variable in y.

#X represents the size of a tumor in centimeters.

#Note: X has to be reshaped into a column from a row for the LogisticRegression() function to work. #y represents whether or not the tumor is cancerous (0 for "No", 1 for "Yes").

We will use a method from the sklearn module, so we will have to import that module as well: from sklearn import linear_model From the sklearn module we will use the LogisticRegression() method to create a logistic regression object.

that takes the independent and dependent values as parameters and fills the regression object with data that describes the relationship:

logr.fit(X,y) Now we have a logistic regression object that is ready to whether a tumor is cancerous based on the tumor size: #predict if tumor is cancerous where the size is 3.46mm:

See the whole example in action: import numpy from sklearn import linear_model #Reshaped for Logistic function.

logr.fit(X,y) #predict if tumor is cancerous where the size is 3.46mm:

[0] We have predicted that a tumor with a size of 3.46mm will not be cancerous.

See the whole example in action: import numpy from sklearn import linear_model #Reshaped for Logistic function.

logr.fit(X,y)

[4.03541657] This tells us that as the size of a tumor increases by 1mm the odds of it being a cancerous tumor increases by 4x.

The coefficient and intercept values can be used to find the probability that each tumor is cancerous. Create a function that uses the model's coefficient and intercept values to return a new value. This new value represents probability that the given observation is a tumor:

To find the log-odds for each observation, we must first create a formula that looks similar to the one from linear regression, extracting the coefficient and the intercept.

Now that we have the odds, we can convert it to probability by dividing it by 1 plus the odds.

Let us now use the function with what we have learned to find out the probability that each tumor is cancerous.

See the whole example in action: import numpy from sklearn import linear_model

logr.fit(X,y)

[[0.60749955] [0.19268876] [0.12775886] [0.00955221] [0.08038616] [0.07345637] [0.88362743] [0.77901378] [0.88924409] [0.81293497] [0.57719129] [0.96664243]]

3.78 0.61 The probability that a tumor with the size 3.78cm is cancerous is 61%. 2.44 0.19 The probability that a tumor with the size 2.44cm is cancerous is 19%. 2.09 0.13 The probability that a tumor with the size 2.09cm is cancerous is 13%.