Tyler Caraza-Harter
In this semester, we've learned two ways to run Python code: in interactive mode, where code runs as we type it, and script mode, where we write our code in a .py file, which we then run.
Now, we're going to learn a third way: with Notebooks. Notebooks are a popular tool of data scientists for interleaving code with analysis results (e.g., tables and plots). You're going to need to install Anaconda for this part of the course (if you haven't already).
To get started, run the following in the terminal (make sure to use the "cd" command to move to a directory where you want to do your work, such as your documents directory):
ipython notebook
If things are working properly, you'll see something like the following in the terminal:
ty-mac:lec-26$ ipython notebook
[TerminalIPythonApp] WARNING | Subcommand `ipython notebook` is deprecated and will be removed in future versions.
[TerminalIPythonApp] WARNING | You likely want to use `jupyter notebook` in the future
[I 10:36:31.799 NotebookApp] JupyterLab beta preview extension loaded from /anaconda3/lib/python3.6/site-packages/jupyterlab
[I 10:36:31.799 NotebookApp] JupyterLab application directory is /anaconda3/share/jupyter/lab
[I 10:36:31.808 NotebookApp] Serving notebooks from local directory: /Users/trh/git_co/caraza-harter-com/tyler/cs301/fall18/materials/readings/lec-26
[I 10:36:31.808 NotebookApp] 0 active kernels
[I 10:36:31.808 NotebookApp] The Jupyter Notebook is running at:
[I 10:36:31.808 NotebookApp] http://localhost:8888/?token=1254607604289711443154541
Also, your web browser should pop up a new page that says Jupyter in the top left. Jupyter Notebooks are an interactive tool for doing data science in Julia, Python, and R (we'll only be working with Python).
Click the "New" button, then choose "Python 3". A new document should open, with a box labeled "In [ ]" in which you can start typing. From here, you can write documentation and Python code, summarize data, and display plots. In fact, this document you're reading now is itself a Python Notebook, and it will include a mix of documentation and code.
We can type Python expressions in these boxes and see the result in an "Out[N]" box. For example:
1 + 2
This should feel similar to interactive Python, but one difference is that we only see the output from the last line. For example, in the following example, we do two computations, but only see the second result.
3*4 # expression 1
10*10 # expression 2
Of course, we can still always print things:
print(3*4)
print(10*10)
One difference between an ipython notebook and interactive mode is that you can go back and edit earlier cells (each "In [N]" box is called a cell). While this is useful, there's one gotcha you need to be aware of.
In interactive mode, this code would fail:
>>> print(x)
>>> x = 3
Similarly, if you try to assign x in a later cell and print it in an earlier cell, it should (and usually does) fail. One catch is that notebooks don't force you to execute each cell in order. You could do something like this:
EARLIER CELL: print("hi")
LATER CELL: x = 3
...Run everything...
Modify EARLIER cell:
EARLIER CELL: print(x)
Run just the EARLIER cell. Surprisingly, (and deceptively), this works, because you ran the cells out of order. However, when you share your notebook with somebody else, it will fail for them because they might try rerunning everything from scratch, in order.
Important: to make sure you don't make this common blunder, be sure to occasionally re-run all your code from the start by clicking "Kernel" from the notebook menu and selecting "Restart & Run All".
A common module you'll use when working in a notebook is called "pandas". Pandas will allow you to do many of the things you might do in Excel (or another spreadsheet program) programmatically with Python. Let's import it.
import pandas as pd
The "as pd" expression may be new for you. This just gives the pandas module a new name in our code, so we can type things like pd.some_function()
to call a function named some_function
rather than type out pandas.some_function
. You could also have just used import pandas
or even given it another name with an import like import pandas as silly_bear
, but we recommend you import as "pd", because most pandas users do so by convention.
We'll also be using two new data types (Series and DataFrame) from Pandas very often, so let's import these directly so we don't even need to prefix their use with "pd.".
from pandas import Series, DataFrame
A Series is very similar to a Python list, and we can convert back and forth between them:
num_list = [100, 200, 300]
print(type(num_list))
num_series = Series(num_list) # create Series from list
print(type(num_series))
# displaying a list:
num_list
# displaying a Series:
num_series
Notice that both the list and the Series contain the same values. However, there are some differences:
"dtype" stands for data type. In this case, it means that the series contains integers, each of which require 64 bits of memory (this detail is not important for us). Although you could create a Series containing different types of data (as with lists), we'll avoid doing so because working with Series of one type will often be more convenient.
Going from a Series back to a list is just as easy as going from a list to a Series:
list(num_series)
With Series, it is easy to apply the same operation to every value in the Series with a single line of code (instead of a loop).
For example, suppose we wanted to add 1 to every item in a list. We would need to write something like this:
orig_nums = [100, 200, 300]
new_nums = []
for x in orig_nums:
new_nums.append(x+1)
new_nums
With a Series, we can do the same like this:
nums = Series([100, 200, 300])
nums + 1
This means multiplication means something very different for lists than for Series.
[1,2,3] * 3
Series([1,2,3]) * 3
Whereas a "+" means concatenate for lists, it means element-wise addition for Series:
[10, 20] + [3, 4]
Series([10, 20]) + Series([3, 4])
One implication of this is that you might not get what you expect if you add Series of different sizes:
Series([10,20,30]) + Series([1,2])
The 10 gets added with the 1, and the 20 gets added with the 2, but there's nothing in the second series to add with 30. 30 plus nothing doesn't make sense, so Pandas gives "NaN". This stands for "Not a Number".
Consider the following:
nums = Series([1,9,8,2])
nums
nums > 5
This example shows that you can do element-wise comparisons as well. The result is a Series of booleans. If the value in the original Series is greater than 5, we see True at the same position in the output Series. Otherwise, the value at the same position in the output Series is False.
We can also chain these operations together:
nums = Series([7,5,8,2,3])
nums
mod_2 = nums % 2
mod_2
odd = mod_2 == 1
odd
As you can see, we first obtained an integer Series (mod_2
) by computing the value of every number modulo 2 (mod_2
) will of course contain only 1's and 0's).
We then create a Boolean series (odd
) by comparing the mod_2
series to 1.
If a number in the nums
Series is odd, then the value at the same position in the odd
series will be True.
Notice what happens when we create a series from a list:
Series([100,200,300])
We see the following:
One interesting difference between lists and Series is that with Series, the index does not always need to correspond so closely with the position; that's just a default that can be overridden.
For example:
nums1 = Series([100, 200, 300], index=[2,1,0])
nums1
Now we see indexes are assigned based on the argument we passed for index (not the position):
When we do element-wise operations between two Sersies, Pandas lines up the data based on index, not position. As a concrete example, consider three Series:
X = Series([100, 200, 300])
Y = Series([10, 20, 30])
Z = Series([10, 20, 30], index=[2,1,0])
X
Y
Z
Note: Y and Z are nearly the same (numbers 10, 20, and 30, in that order), except for the index. Let's see the difference between X+Y
and Y+Z
:
X+Y
X+Z
For X+Y
, Pandas adds the number at index 0 in X (100) with the value at index 0 in Y (10), such that the value in the output at index 0 is 110.
For X+Z
, Pandas adds the number at index 0 in X (100) with the value at index 0 in Y (30), such that the value in the output at index 0 is 130. It doesn't matter that the first number in Z is 10, because Pandas does element-wise operations based on index, not position.
We will often want to extract one or more values from a Series. Like lists, Series support indexing and slicing.
letter_list = ["A", "B", "C", "D"]
letter_series = Series(letter_list)
letter_series
letter_list[0]
letter_series[0]
letter_list[3]
letter_series[3]
letter_list[-1]
# but be careful! Series don't support negative indexes like lists do
try:
print(letter_series[-1])
except Exception as e:
print(type(e))
Series slicing works much like list slicing:
print("list slice:")
print(letter_list[:2])
print("\nseries slice:")
print(letter_series[:2])
print("list slice:")
print(letter_list[2:])
print("\nseries slice:")
print(letter_series[2:])
# although we CANNOT do negative indexing with a Series
# we CAN use negative numbers in a Series slice
print("list slice:")
print(letter_list[:-1])
print("\nseries slice:")
print(letter_series[:-1])
In addition to supporting indexing and slicing, Series allow you to pull out items using a Boolean series. For example:
letters = Series(["A", "B", "C", "D"])
letters
bool_series = Series([True, True, False, False])
bool_series
# we can used the bool_series almost like an index
# to pull values out of letters:
letters[bool_series]
# We could also create the Boolean Series on the fly:
letters[Series([True, True, False, False])]
# Let's grab the last two letterrs:
letters[Series([False, False, True, True])]
# Let's grab the first and last (can't do this with a slice):
letters[Series([True, False, False, True])]
As we just saw, we can use a Boolean series (let's call it B) to select values from another Series (let's call it S).
A common pattern is to create B by performing operation on S, then using B to select from S. Let's try doing this to pull all the numbers greater than 5 from a Series.
Example 1
# we want to pull out 9 and 8
S = Series([1,9,2,3,8])
S
B = S > 5
B
# this will pull out values from S at index 1 and 4,
# because the values in B at index 1 and 4 are True
S[B]
Example 2
Let's try to pull out all the upper case strings from a series:
words = Series(["APPLE", "boy", "CAT", "dog"])
words
# we can use .str.upper() to get upper case version of words
upper_words = words.str.upper()
upper_words
# B will be True where the original word equals the upper-case version
B = words == upper_words
B
# pull out the just words that were orginally uppercase
words[B]
We have done this example in several steps to illustrate what is happening, but it could have been simplified. Recall that B is words == upper_words
. Thus we could have done this without ever storing a Boolean series in B:
words[words == upper_words]
Let's simplify one step further (instead of using upper_words, let's paste the expression we used to compute it earlier):
words[words == words.str.upper()]
Example 3
Let's try to pull out all the odd numbers from this Series:
nums = Series([11,12,19,18,15,17])
nums
nums % 2
well produce a Series of 1's (for odd numbers) and 0's (for even numbers). Thus nums % 2 == 1
produces a Boolean Series of True's (for odd numbers) and False's (for even numbers). Let's use that Boolean Series to pull out the odd numbers:
nums[nums % 2 == 1]
Pandas will often be used to deal with tabular data (much as in Excel).
In many tables, all the data in the same column is similar, so Pandas represents each column in a table as a Series object. A table is represented as a DataFrame, which is just a collection of named Series (one for each column).
We can use a dictionary of aligned Series objects to create a dictionary. For example:
name_column = Series(["Alice", "Bob", "Cindy", "Dan"])
score_column = Series([100, 150, 160, 120])
table = DataFrame({'name': name_column, 'score': score_column})
table
Or, if we want, we can create a DataFrame table from a dictionary of lists, and Pandas will implicitly create the Series for each column for us:
data = {"name": ["Alice", "Bob", "Cindy", "Dan"],
"score": [100, 150, 160, 120]}
df = DataFrame(data)
df
There are a few things we might want to do:
# we'll use the DataFrame of scores defined
# in the previous section
df
# let's grab the name cell using DataFrame["COL NAME"]
df["name"]
# or we could extract the score column:
df["score"]
# if we want to generate some simple stats over a column,
# we can use .describe()
df["score"].describe()
# lookup is done for columns by default (df[x] looks up column named x)
# we can also lookup a row, but we need to use df.loc[y]. ("loc" stands for location)
# for example, let's get Bob's row:
df.loc[1]
# if we want a particular cell, we can use df.loc[row,col].
# for example, this is Bob's score:
df.loc[1, "score"]
# we can also use this to modify cells:
df.loc[1, "score"] += 5
df
Most of the time, we'll let Pandas directly load a CSV file to a DataFrame (instead of creating a dictionary of lists ourselves). We can easily do this with pd.read_csv(path)
(recall that we imported pandas as import pandas as pd
):
# movies is a DataFrame
movies = pd.read_csv('IMDB-Movie-Data.csv')
# how many are there?
print("Number of movies:", len(movies))
# it's large, but we can preview the first few with DataFrame.head()
movies.head()
# we can pull out Runtime minutes if we like
runtime = movies["Runtime"]
# it's still long (same length as movies), but let's preview the first 10 runtime minutes
runtime.head(10)
# what is the mean runtime, in hours?
runtime.mean() / 60
# what if we want stats about movies from 2016?
# use .head() on results to make it shorter
(movies["Year"] == 2016).head()
Observe:
Let's pull out the movies from 2016 using this Boolean Series:
movies_2016 = movies[movies["Year"] == 2016]
print("there are " + str(len(movies_2016)) + " movies in 2016")
movies_2016.head(10)
# let's get some general stats about movies from 2016
movies_2016.describe()
We see (amonge other things) that the average Runtime is 107.34 minutes.
Writing Python in an interactive Jupyter notebook is a powerful way to write code and see the data we're working with in the same place. Data comes in many different forms, but tabular data is especially common. The Pandas module helps us work with tabular data and integrates with ipython, making it fast and easy to compute simple statistics over columns within our dataset. In this lesson, we learned to do the following: