Linear Regression is a very commonly used statistical method that allows us to determine and study the relationship between two continuous variables. The various properties of linear regression and its Python implementation have been covered in this article previously. Now, we shall find out how to implement this in PyTorch, a very popular deep learning library that is being developed by Facebook.
Firstly, you will need to install PyTorch into your Python environment. The easiest way to do this is to use the pip or conda tool. Visit pytorch.org and install the version of your Python interpreter and the package manager that you would like to use.
Python3
# We can run this Python code on a Jupyter notebook# to automatically install the correct version of# PyTorch.# http://pytorch.org / from os import pathfrom wheel.pep425tags import get_abbr_impl, get_impl_ver, get_abi_tagplatform = '{}{}-{}'.format(get_abbr_impl(), get_impl_ver(), get_abi_tag())accelerator = 'cu80' if path.exists('/opt / bin / nvidia-smi') else 'cpu'! pip install -q http://download.pytorch.org / whl/{accelerator}/torch-1.3.1.post4-{platform}-linux_x86_64.whl torchvision |
With PyTorch installed, let us now have a look at the code.
Write the two lines given below to import the necessary library functions and objects.
Python3
import torchfrom torch.autograd import Variable |
We also define some data and assign them to variables x_data and y_data as given below:
Python3
x_data = Variable(torch.Tensor([[1.0], [2.0], [3.0]]))y_data = Variable(torch.Tensor([[2.0], [4.0], [6.0]])) |
Here, x_data is our independent variable and y_data is our dependent variable. This will be our dataset for now. Next, we need to define our model. There are two main steps associated with defining our model. They are:
- Initializing our model.
- Declaring the forward pass.
We use the class given below:
Python3
class LinearRegressionModel(torch.nn.Module): def __init__(self): super(LinearRegressionModel, self).__init__() self.linear = torch.nn.Linear(1, 1) # One in and one out def forward(self, x): y_pred = self.linear(x) return y_pred |
As you can see, our Model class is a subclass of torch.nn.module. Also, since here we have only one input and one output, we use a Linear model with both the input and output dimension as 1.
Next, we create an object of this model.
Python3
# our modelour_model = LinearRegressionModel() |
After this, we select the optimizer and the loss criteria. Here, we will use the mean squared error (MSE) as our loss function and stochastic gradient descent (SGD) as our optimizer. Also, we arbitrarily fix a learning rate of 0.01.
Python3
criterion = torch.nn.MSELoss(size_average = False)optimizer = torch.optim.SGD(our_model.parameters(), lr = 0.01) |
We now arrive at our training step. We perform the following tasks 500 times during training:
- Perform a forward pass bypassing our data and finding out the predicted value of y.
- Compute the loss using MSE.
- Reset all the gradients to 0, perform a backpropagation and then, update the weights.
Python3
for epoch in range(500): # Forward pass: Compute predicted y by passing # x to the model pred_y = our_model(x_data) # Compute and print loss loss = criterion(pred_y, y_data) # Zero gradients, perform a backward pass, # and update the weights. optimizer.zero_grad() loss.backward() optimizer.step() print('epoch {}, loss {}'.format(epoch, loss.item())) |
Once the training is completed, we test if we are getting correct results using the model that we defined. So, we test it for an unknown value of x_data, in this case, 4.0.
Python3
new_var = Variable(torch.Tensor([[4.0]]))pred_y = our_model(new_var)print("predict (after training)", 4, our_model(new_var).item()) |
If you performed all steps correctly, you will see that for input 4.0, you are getting a value that is very close to 8.0 as below. So, our model inherently learns the relationship between the input data and the output data without being programmed explicitly.
predict (after training) 4 7.966438293457031
For your reference, you can find the entire code of this article given below:
Python3
import torchfrom torch.autograd import Variablex_data = Variable(torch.Tensor([[1.0], [2.0], [3.0]]))y_data = Variable(torch.Tensor([[2.0], [4.0], [6.0]]))class LinearRegressionModel(torch.nn.Module): def __init__(self): super(LinearRegressionModel, self).__init__() self.linear = torch.nn.Linear(1, 1) # One in and one out def forward(self, x): y_pred = self.linear(x) return y_pred# our modelour_model = LinearRegressionModel()criterion = torch.nn.MSELoss(size_average = False)optimizer = torch.optim.SGD(our_model.parameters(), lr = 0.01)for epoch in range(500): # Forward pass: Compute predicted y by passing # x to the model pred_y = our_model(x_data) # Compute and print loss loss = criterion(pred_y, y_data) # Zero gradients, perform a backward pass, # and update the weights. optimizer.zero_grad() loss.backward() optimizer.step() print('epoch {}, loss {}'.format(epoch, loss.item()))new_var = Variable(torch.Tensor([[4.0]]))pred_y = our_model(new_var)print("predict (after training)", 4, our_model(new_var).item()) |
