from IPython.core.display import HTML
HTML("<iframe src=http://pandas.pydata.org width=800 height=350></iframe>")

%matplotlib inline
import pandas as pd
import numpy as np

# Set some Pandas options
pd.set_option('html', False)
pd.set_option('max_columns', 30)
pd.set_option('max_rows', 20)

"""
Explanation: Introduction to Pandas
pandas is a Python package providing fast, flexible, and expressive data structures designed to work with relational or labeled data both. It is a fundamental high-level building block for doing practical, real world data analysis in Python. 
pandas is well suited for:

Tabular data with heterogeneously-typed columns, as in an SQL table or Excel spreadsheet
Ordered and unordered (not necessarily fixed-frequency) time series data.
Arbitrary matrix data (homogeneously typed or heterogeneous) with row and column labels
Any other form of observational / statistical data sets. The data actually need not be labeled at all to be placed into a pandas data structure

Key features:

Easy handling of missing data
Size mutability: columns can be inserted and deleted from DataFrame and higher dimensional objects
Automatic and explicit data alignment: objects can be explicitly aligned to a set of labels, or the data can be aligned automatically
Powerful, flexible group by functionality to perform split-apply-combine operations on data sets
Intelligent label-based slicing, fancy indexing, and subsetting of large data sets
Intuitive merging and joining data sets
Flexible reshaping and pivoting of data sets
Hierarchical labeling of axes
Robust IO tools for loading data from flat files, Excel files, databases, and HDF5
Time series functionality: date range generation and frequency conversion, moving window statistics, moving window linear regressions, date shifting and lagging, etc.
End of explanation
"""


counts = pd.Series([632, 1638, 569, 115])
counts

"""
Explanation: Pandas Data Structures
Series
A Series is a single vector of data (like a NumPy array) with an index that labels each element in the vector.
End of explanation
"""


counts.values

counts.index

"""
Explanation: If an index is not specified, a default sequence of integers is assigned as the index. A NumPy array comprises the values of the Series, while the index is a pandas Index object.
End of explanation
"""


bacteria = pd.Series([632, 1638, 569, 115], 
    index=['Firmicutes', 'Proteobacteria', 'Actinobacteria', 'Bacteroidetes'])

bacteria

"""
Explanation: We can assign meaningful labels to the index, if they are available:
End of explanation
"""


bacteria['Actinobacteria']

bacteria[[name.endswith('bacteria') for name in bacteria.index]]

[name.endswith('bacteria') for name in bacteria.index]

"""
Explanation: These labels can be used to refer to the values in the Series.
End of explanation
"""


bacteria[0]

"""
Explanation: Notice that the indexing operation preserved the association between the values and the corresponding indices.
We can still use positional indexing if we wish.
End of explanation
"""


bacteria.name = 'counts'
bacteria.index.name = 'phylum'
bacteria

"""
Explanation: We can give both the array of values and the index meaningful labels themselves:
End of explanation
"""


np.log(bacteria)

"""
Explanation: NumPy's math functions and other operations can be applied to Series without losing the data structure.
End of explanation
"""


bacteria[bacteria>1000]

"""
Explanation: We can also filter according to the values in the Series:
End of explanation
"""


bacteria_dict = {'Firmicutes': 632, 'Proteobacteria': 1638, 'Actinobacteria': 569, 'Bacteroidetes': 115}
pd.Series(bacteria_dict)

"""
Explanation: A Series can be thought of as an ordered key-value store. In fact, we can create one from a dict:
End of explanation
"""


bacteria2 = pd.Series(bacteria_dict, index=['Cyanobacteria','Firmicutes','Proteobacteria','Actinobacteria'])
bacteria2

bacteria2.isnull()

"""
Explanation: Notice that the Series is created in key-sorted order.
If we pass a custom index to Series, it will select the corresponding values from the dict, and treat indices without corrsponding values as missing. Pandas uses the NaN (not a number) type for missing values.
End of explanation
"""


bacteria + bacteria2

"""
Explanation: Critically, the labels are used to align data when used in operations with other Series objects:
End of explanation
"""


data = pd.DataFrame({'value':[632, 1638, 569, 115, 433, 1130, 754, 555],
                     'patient':[1, 1, 1, 1, 2, 2, 2, 2],
                     'phylum':['Firmicutes', 'Proteobacteria', 'Actinobacteria', 
    'Bacteroidetes', 'Firmicutes', 'Proteobacteria', 'Actinobacteria', 'Bacteroidetes']})
data

"""
Explanation: Contrast this with NumPy arrays, where arrays of the same length will combine values element-wise; adding Series combined values with the same label in the resulting series. Notice also that the missing values were propogated by addition.
DataFrame
Inevitably, we want to be able to store, view and manipulate data that is multivariate, where for every index there are multiple fields or columns of data, often of varying data type.
A DataFrame is a tabular data structure, encapsulating multiple series like columns in a spreadsheet. Data are stored internally as a 2-dimensional object, but the DataFrame allows us to represent and manipulate higher-dimensional data.
End of explanation
"""


data[['phylum','value','patient']]

"""
Explanation: Notice the DataFrame is sorted by column name. We can change the order by indexing them in the order we desire:
End of explanation
"""


data.columns

"""
Explanation: A DataFrame has a second index, representing the columns:
End of explanation
"""


data['value']

data.value

type(data.value)

type(data[['value']])

"""
Explanation: If we wish to access columns, we can do so either by dict-like indexing or by attribute:
End of explanation
"""


data.ix[3]

"""
Explanation: Notice this is different than with Series, where dict-like indexing retrieved a particular element (row). If we want access to a row in a DataFrame, we index its ix attribute.
End of explanation
"""


data = pd.DataFrame({0: {'patient': 1, 'phylum': 'Firmicutes', 'value': 632},
                    1: {'patient': 1, 'phylum': 'Proteobacteria', 'value': 1638},
                    2: {'patient': 1, 'phylum': 'Actinobacteria', 'value': 569},
                    3: {'patient': 1, 'phylum': 'Bacteroidetes', 'value': 115},
                    4: {'patient': 2, 'phylum': 'Firmicutes', 'value': 433},
                    5: {'patient': 2, 'phylum': 'Proteobacteria', 'value': 1130},
                    6: {'patient': 2, 'phylum': 'Actinobacteria', 'value': 754},
                    7: {'patient': 2, 'phylum': 'Bacteroidetes', 'value': 555}})

data

"""
Explanation: Alternatively, we can create a DataFrame with a dict of dicts:
End of explanation
"""


data = data.T
data

"""
Explanation: We probably want this transposed:
End of explanation
"""


vals = data.value
vals

vals[5] = 0
vals

data

vals = data.value.copy()
vals[5] = 1000
data

"""
Explanation: Its important to note that the Series returned when a DataFrame is indexted is merely a view on the DataFrame, and not a copy of the data itself. So you must be cautious when manipulating this data:
End of explanation
"""


data.value[3] = 14
data

data['year'] = 2013
data

"""
Explanation: We can create or modify columns by assignment:
End of explanation
"""


data.treatment = 1
data

data.treatment

"""
Explanation: But note, we cannot use the attribute indexing method to add a new column:
End of explanation
"""


treatment = pd.Series([0]*4 + [1]*2)
treatment

data['treatment'] = treatment
data

"""
Explanation: Specifying a Series as a new columns cause its values to be added according to the DataFrame's index:
End of explanation
"""


month = ['Jan', 'Feb', 'Mar', 'Apr']
data['month'] = month

data['month'] = ['Jan']*len(data)
data

"""
Explanation: Other Python data structures (ones without an index) need to be the same length as the DataFrame:
End of explanation
"""


del data['month']
data

"""
Explanation: We can use del to remove columns, in the same way dict entries can be removed:
End of explanation
"""


data.values

"""
Explanation: We can extract the underlying data as a simple ndarray by accessing the values attribute:
End of explanation
"""


df = pd.DataFrame({'foo': [1,2,3], 'bar':[0.4, -1.0, 4.5]})
df.values

"""
Explanation: Notice that because of the mix of string and integer (and NaN) values, the dtype of the array is object. The dtype will automatically be chosen to be as general as needed to accomodate all the columns.
End of explanation
"""


data.index

"""
Explanation: Pandas uses a custom data structure to represent the indices of Series and DataFrames.
End of explanation
"""


data.index[0] = 15

"""
Explanation: Index objects are immutable:
End of explanation
"""


bacteria2.index = bacteria.index

bacteria2

"""
Explanation: This is so that Index objects can be shared between data structures without fear that they will be changed.
End of explanation
"""


!cat data/microbiome.csv

"""
Explanation: Importing data
A key, but often under-appreciated, step in data analysis is importing the data that we wish to analyze. Though it is easy to load basic data structures into Python using built-in tools or those provided by packages like NumPy, it is non-trivial to import structured data well, and to easily convert this input into a robust data structure:
genes = np.loadtxt("genes.csv", delimiter=",", dtype=[('gene', '|S10'), ('value', '&lt;f4')])

Pandas provides a convenient set of functions for importing tabular data in a number of formats directly into a DataFrame object. These functions include a slew of options to perform type inference, indexing, parsing, iterating and cleaning automatically as data are imported.
Let's start with some more bacteria data, stored in csv format.
End of explanation
"""


mb = pd.read_csv("data/microbiome.csv")
mb

"""
Explanation: This table can be read into a DataFrame using read_csv:
End of explanation
"""


pd.read_csv("data/microbiome.csv", header=None).head()

"""
Explanation: Notice that read_csv automatically considered the first row in the file to be a header row.
We can override default behavior by customizing some the arguments, like header, names or index_col.
End of explanation
"""


mb = pd.read_table("data/microbiome.csv", sep=',')

"""
Explanation: read_csv is just a convenience function for read_table, since csv is such a common format:
End of explanation
"""


mb = pd.read_csv("data/microbiome.csv", index_col=['Taxon','Patient'])
mb.head()

"""
Explanation: The sep argument can be customized as needed to accomodate arbitrary separators. For example, we can use a regular expression to define a variable amount of whitespace, which is unfortunately very common in some data formats: 
sep='\s+'

For a more useful index, we can specify the first two columns, which together provide a unique index to the data.
End of explanation
"""


pd.read_csv("data/microbiome.csv", skiprows=[3,4,6]).head()

"""
Explanation: This is called a hierarchical index, which we will revisit later in the tutorial.
If we have sections of data that we do not wish to import (for example, known bad data), we can populate the skiprows argument:
End of explanation
"""


pd.read_csv("data/microbiome.csv", nrows=4)

"""
Explanation: Conversely, if we only want to import a small number of rows from, say, a very large data file we can use nrows:
End of explanation
"""


data_chunks = pd.read_csv("data/microbiome.csv", chunksize=15)

mean_tissue = {chunk.Taxon[0]:chunk.Tissue.mean() for chunk in data_chunks}
    
mean_tissue

"""
Explanation: Alternately, if we want to process our data in reasonable chunks, the chunksize argument will return an iterable object that can be employed in a data processing loop. For example, our microbiome data are organized by bacterial phylum, with 15 patients represented in each:
End of explanation
"""


!cat data/microbiome_missing.csv

pd.read_csv("data/microbiome_missing.csv").head(20)

"""
Explanation: Most real-world data is incomplete, with values missing due to incomplete observation, data entry or transcription error, or other reasons. Pandas will automatically recognize and parse common missing data indicators, including NA and NULL.
End of explanation
"""


pd.isnull(pd.read_csv("data/microbiome_missing.csv")).head(20)

"""
Explanation: Above, Pandas recognized NA and an empty field as missing data.
End of explanation
"""


pd.read_csv("data/microbiome_missing.csv", na_values=['?', -99999]).head(20)

"""
Explanation: Unfortunately, there will sometimes be inconsistency with the conventions for missing data. In this example, there is a question mark "?" and a large negative number where there should have been a positive integer. We can specify additional symbols with the na_values argument:
End of explanation
"""


mb_file = pd.ExcelFile('data/microbiome/MID1.xls')
mb_file

"""
Explanation: These can be specified on a column-wise basis using an appropriate dict as the argument for na_values.
Microsoft Excel
Since so much financial and scientific data ends up in Excel spreadsheets (regrettably), Pandas' ability to directly import Excel spreadsheets is valuable. This support is contingent on having one or two dependencies (depending on what version of Excel file is being imported) installed: xlrd and openpyxl (these may be installed with either pip or easy_install).
Importing Excel data to Pandas is a two-step process. First, we create an ExcelFile object using the path of the file:
End of explanation
"""


mb1 = mb_file.parse("Sheet 1", header=None)
mb1.columns = ["Taxon", "Count"]
mb1.head()

"""
Explanation: Then, since modern spreadsheets consist of one or more "sheets", we parse the sheet with the data of interest:
End of explanation
"""


mb2 = pd.read_excel('data/microbiome/MID2.xls', sheetname='Sheet 1', header=None)
mb2.head()

"""
Explanation: There is now a read_excel conveneince function in Pandas that combines these steps into a single call:
End of explanation
"""


baseball = pd.read_csv("data/baseball.csv", index_col='id')
baseball.head()

"""
Explanation: There are several other data formats that can be imported into Python and converted into DataFrames, with the help of buitl-in or third-party libraries. These include JSON, XML, HDF5, relational and non-relational databases, and various web APIs. These are beyond the scope of this tutorial, but are covered in Python for Data Analysis.
Pandas Fundamentals
This section introduces the new user to the key functionality of Pandas that is required to use the software effectively.
For some variety, we will leave our digestive tract bacteria behind and employ some baseball data.
End of explanation
"""


player_id = baseball.player + baseball.year.astype(str)
baseball_newind = baseball.copy()
baseball_newind.index = player_id
baseball_newind.head()

"""
Explanation: Notice that we specified the id column as the index, since it appears to be a unique identifier. We could try to create a unique index ourselves by combining player and year:
End of explanation
"""


baseball_newind.index.is_unique

"""
Explanation: This looks okay, but let's check:
End of explanation
"""


pd.Series(baseball_newind.index).value_counts()

"""
Explanation: So, indices need not be unique. Our choice is not unique because some players change teams within years.
End of explanation
"""


baseball_newind.ix['wickmbo012007']

"""
Explanation: The most important consequence of a non-unique index is that indexing by label will return multiple values for some labels:
End of explanation
"""


player_unique = baseball.player + baseball.team + baseball.year.astype(str)
baseball_newind = baseball.copy()
baseball_newind.index = player_unique
baseball_newind.head()

baseball_newind.index.is_unique

"""
Explanation: We will learn more about indexing below.
We can create a truly unique index by combining player, team and year:
End of explanation
"""


baseball.reindex(baseball.index[::-1]).head()

"""
Explanation: We can create meaningful indices more easily using a hierarchical index; for now, we will stick with the numeric id field as our index.
Manipulating indices
Reindexing allows users to manipulate the data labels in a DataFrame. It forces a DataFrame to conform to the new index, and optionally, fill in missing data if requested.
A simple use of reindex is to alter the order of the rows:
End of explanation
"""


id_range = range(baseball.index.values.min(), baseball.index.values.max())
baseball.reindex(id_range).head()

"""
Explanation: Notice that the id index is not sequential. Say we wanted to populate the table with every id value. We could specify and index that is a sequence from the first to the last id numbers in the database, and Pandas would fill in the missing data with NaN values:
End of explanation
"""


baseball.reindex(id_range, method='ffill', columns=['player','year']).head()

baseball.reindex(id_range, fill_value='mr.nobody', columns=['player']).head()

"""
Explanation: Missing values can be filled as desired, either with selected values, or by rule:
End of explanation
"""


baseball.shape

baseball.drop([89525, 89526])

baseball.drop(['ibb','hbp'], axis=1)

"""
Explanation: Keep in mind that reindex does not work if we pass a non-unique index series.
We can remove rows or columns via the drop method:
End of explanation
"""


# Sample Series object
hits = baseball_newind.h
hits

# Numpy-style indexing
hits[:3]

# Indexing by label
hits[['womacto01CHN2006','schilcu01BOS2006']]

"""
Explanation: Indexing and Selection
Indexing works analogously to indexing in NumPy arrays, except we can use the labels in the Index object to extract values in addition to arrays of integers.
End of explanation
"""


hits['womacto01CHN2006':'gonzalu01ARI2006']

hits['womacto01CHN2006':'gonzalu01ARI2006'] = 5
hits

"""
Explanation: We can also slice with data labels, since they have an intrinsic order within the Index:
End of explanation
"""


baseball_newind[['h','ab']]

baseball_newind[baseball_newind.ab>500]

"""
Explanation: In a DataFrame we can slice along either or both axes:
End of explanation
"""


baseball_newind.ix['gonzalu01ARI2006', ['h','X2b', 'X3b', 'hr']]

baseball_newind.ix[['gonzalu01ARI2006','finlest01SFN2006'], 5:8]

baseball_newind.ix[:'myersmi01NYA2006', 'hr']

"""
Explanation: The indexing field ix allows us to select subsets of rows and columns in an intuitive way:
End of explanation
"""


baseball_newind.xs('myersmi01NYA2006')

"""
Explanation: Similarly, the cross-section method xs (not a field) extracts a single column or row by label and returns it as a Series:
End of explanation
"""


hr2006 = baseball[baseball.year==2006].xs('hr', axis=1)
hr2006.index = baseball.player[baseball.year==2006]

hr2007 = baseball[baseball.year==2007].xs('hr', axis=1)
hr2007.index = baseball.player[baseball.year==2007]

hr2006 = pd.Series(baseball.hr[baseball.year==2006].values, index=baseball.player[baseball.year==2006])
hr2007 = pd.Series(baseball.hr[baseball.year==2007].values, index=baseball.player[baseball.year==2007])

hr_total = hr2006 + hr2007
hr_total

"""
Explanation: Operations
DataFrame and Series objects allow for several operations to take place either on a single object, or between two or more objects.
For example, we can perform arithmetic on the elements of two objects, such as combining baseball statistics across years:
End of explanation
"""


hr_total[hr_total.notnull()]

"""
Explanation: Pandas' data alignment places NaN values for labels that do not overlap in the two Series. In fact, there are only 6 players that occur in both years.
End of explanation
"""


hr2007.add(hr2006, fill_value=0)

"""
Explanation: While we do want the operation to honor the data labels in this way, we probably do not want the missing values to be filled with NaN. We can use the add method to calculate player home run totals by using the fill_value argument to insert a zero for home runs where labels do not overlap:
End of explanation
"""


baseball.hr - baseball.hr.max()

"""
Explanation: Operations can also be broadcast between rows or columns.
For example, if we subtract the maximum number of home runs hit from the hr column, we get how many fewer than the maximum were hit by each player:
End of explanation
"""


baseball.ix[89521]["player"]

stats = baseball[['h','X2b', 'X3b', 'hr']]
diff = stats - stats.xs(89521)
diff[:10]

"""
Explanation: Or, looking at things row-wise, we can see how a particular player compares with the rest of the group with respect to important statistics
End of explanation
"""


stats.apply(np.median)

stat_range = lambda x: x.max() - x.min()
stats.apply(stat_range)

"""
Explanation: We can also apply functions to each column or row of a DataFrame
End of explanation
"""


slg = lambda x: (x['h']-x['X2b']-x['X3b']-x['hr'] + 2*x['X2b'] + 3*x['X3b'] + 4*x['hr'])/(x['ab']+1e-6)
baseball.apply(slg, axis=1).apply(lambda x: '%.3f' % x)

"""
Explanation: Lets use apply to calculate a meaningful baseball statistics, slugging percentage:
$$SLG = \frac{1B + (2 \times 2B) + (3 \times 3B) + (4 \times HR)}{AB}$$
And just for fun, we will format the resulting estimate.
End of explanation
"""


baseball_newind.sort_index().head()

baseball_newind.sort_index(ascending=False).head()

baseball_newind.sort_index(axis=1).head()

"""
Explanation: Sorting and Ranking
Pandas objects include methods for re-ordering data.
End of explanation
"""


baseball.hr.order(ascending=False)

"""
Explanation: We can also use order to sort a Series by value, rather than by label.
End of explanation
"""


baseball[['player','sb','cs']].sort_index(ascending=[False,True], by=['sb', 'cs']).head(10)

"""
Explanation: For a DataFrame, we can sort according to the values of one or more columns using the by argument of sort_index:
End of explanation
"""


baseball.hr.rank()

"""
Explanation: Ranking does not re-arrange data, but instead returns an index that ranks each value relative to others in the Series.
End of explanation
"""


pd.Series([100,100]).rank()

"""
Explanation: Ties are assigned the mean value of the tied ranks, which may result in decimal values.
End of explanation
"""


baseball.hr.rank(method='first')

"""
Explanation: Alternatively, you can break ties via one of several methods, such as by the order in which they occur in the dataset:
End of explanation
"""


baseball.rank(ascending=False).head()

baseball[['r','h','hr']].rank(ascending=False).head()

"""
Explanation: Calling the DataFrame's rank method results in the ranks of all columns:
End of explanation
"""


# Write your answer here

"""
Explanation: Exercise
Calculate on base percentage for each player, and return the ordered series of estimates.
$$OBP = \frac{H + BB + HBP}{AB + BB + HBP + SF}$$
End of explanation
"""


baseball_h = baseball.set_index(['year', 'team', 'player'])
baseball_h.head(10)

"""
Explanation: Hierarchical indexing
In the baseball example, I was forced to combine 3 fields to obtain a unique index that was not simply an integer value. A more elegant way to have done this would be to create a hierarchical index from the three fields.
End of explanation
"""


baseball_h.index[:10]

baseball_h.index.is_unique

baseball_h.ix[(2007, 'ATL', 'francju01')]

"""
Explanation: This index is a MultiIndex object that consists of a sequence of tuples, the elements of which is some combination of the three columns used to create the index. Where there are multiple repeated values, Pandas does not print the repeats, making it easy to identify groups of values.
End of explanation
"""


mb = pd.read_csv("data/microbiome.csv", index_col=['Taxon','Patient'])

mb.head(10)

mb.index

"""
Explanation: Recall earlier we imported some microbiome data using two index columns. This created a 2-level hierarchical index:
End of explanation
"""


mb.ix['Proteobacteria']

"""
Explanation: With a hierachical index, we can select subsets of the data based on a partial index:
End of explanation
"""


frame = pd.DataFrame(np.arange(12).reshape(( 4, 3)), 
                  index =[['a', 'a', 'b', 'b'], [1, 2, 1, 2]], 
                  columns =[['Ohio', 'Ohio', 'Colorado'], ['Green', 'Red', 'Green']])

frame

"""
Explanation: Hierarchical indices can be created on either or both axes. Here is a trivial example:
End of explanation
"""


frame.index.names = ['key1', 'key2']
frame.columns.names = ['state', 'color']
frame

"""
Explanation: If you want to get fancy, both the row and column indices themselves can be given names:
End of explanation
"""


frame.ix['a']['Ohio']

frame.ix['b', 2]['Colorado']

"""
Explanation: With this, we can do all sorts of custom indexing:
End of explanation
"""


mb.swaplevel('Patient', 'Taxon').head()

"""
Explanation: Additionally, the order of the set of indices in a hierarchical MultiIndex can be changed by swapping them pairwise:
End of explanation
"""


mb.sortlevel('Patient', ascending=False).head()

"""
Explanation: Data can also be sorted by any index level, using sortlevel:
End of explanation
"""


foo = pd.Series([np.nan, -3, None, 'foobar'])
foo

foo.isnull()

"""
Explanation: Missing data
The occurence of missing data is so prevalent that it pays to use tools like Pandas, which seamlessly integrates missing data handling so that it can be dealt with easily, and in the manner required by the analysis at hand.
Missing data are represented in Series and DataFrame objects by the NaN floating point value. However, None is also treated as missing, since it is commonly used as such in other contexts (e.g. NumPy).
End of explanation
"""


bacteria2

bacteria2.dropna()

bacteria2[bacteria2.notnull()]

"""
Explanation: Missing values may be dropped or indexed out:
End of explanation
"""


data

data.dropna()

"""
Explanation: By default, dropna drops entire rows in which one or more values are missing.
End of explanation
"""


data.dropna(how='all')

"""
Explanation: This can be overridden by passing the how='all' argument, which only drops a row when every field is a missing value.
End of explanation
"""


data.ix[7, 'year'] = np.nan
data

data.dropna(thresh=4)

"""
Explanation: This can be customized further by specifying how many values need to be present before a row is dropped via the thresh argument.
End of explanation
"""


data.dropna(axis=1)

"""
Explanation: This is typically used in time series applications, where there are repeated measurements that are incomplete for some subjects.
If we want to drop missing values column-wise instead of row-wise, we use axis=1.
End of explanation
"""


bacteria2.fillna(0)

data.fillna({'year': 2013, 'treatment':2})

"""
Explanation: Rather than omitting missing data from an analysis, in some cases it may be suitable to fill the missing value in, either with a default value (such as zero) or a value that is either imputed or carried forward/backward from similar data points. We can do this programmatically in Pandas with the fillna argument.
End of explanation
"""


data

"""
Explanation: Notice that fillna by default returns a new object with the desired filling behavior, rather than changing the Series or  DataFrame in place (in general, we like to do this, by the way!).
End of explanation
"""


_ = data.year.fillna(2013, inplace=True)
data

"""
Explanation: We can alter values in-place using inplace=True.
End of explanation
"""


bacteria2.fillna(method='bfill')

bacteria2.fillna(bacteria2.mean())

"""
Explanation: Missing values can also be interpolated, using any one of a variety of methods:
End of explanation
"""


baseball.sum()

"""
Explanation: Data summarization
We often wish to summarize data in Series or DataFrame objects, so that they can more easily be understood or compared with similar data. The NumPy package contains several functions that are useful here, but several summarization or reduction methods are built into Pandas data structures.
End of explanation
"""


baseball.mean()

"""
Explanation: Clearly, sum is more meaningful for some columns than others. For methods like mean for which application to string variables is not just meaningless, but impossible, these columns are automatically exculded:
End of explanation
"""


bacteria2

bacteria2.mean()

"""
Explanation: The important difference between NumPy's functions and Pandas' methods is that the latter have built-in support for handling missing data.
End of explanation
"""


bacteria2.mean(skipna=False)

"""
Explanation: Sometimes we may not want to ignore missing values, and allow the nan to propagate.
End of explanation
"""


extra_bases = baseball[['X2b','X3b','hr']].sum(axis=1)
extra_bases.order(ascending=False)

"""
Explanation: Passing axis=1 will summarize over rows instead of columns, which only makes sense in certain situations.
End of explanation
"""


baseball.describe()

"""
Explanation: A useful summarization that gives a quick snapshot of multiple statistics for a Series or DataFrame is describe:
End of explanation
"""


baseball.player.describe()

"""
Explanation: describe can detect non-numeric data and sometimes yield useful information about it.
End of explanation
"""


baseball.hr.cov(baseball.X2b)

"""
Explanation: We can also calculate summary statistics across multiple columns, for example, correlation and covariance.
$$cov(x,y) = \sum_i (x_i - \bar{x})(y_i - \bar{y})$$
End of explanation
"""


baseball.hr.corr(baseball.X2b)

baseball.ab.corr(baseball.h)

baseball.corr()

"""
Explanation: $$corr(x,y) = \frac{cov(x,y)}{(n-1)s_x s_y} = \frac{\sum_i (x_i - \bar{x})(y_i - \bar{y})}{\sqrt{\sum_i (x_i - \bar{x})^2 \sum_i (y_i - \bar{y})^2}}$$
End of explanation
"""


mb.head()

mb.sum(level='Taxon')

"""
Explanation: If we have a DataFrame with a hierarchical index (or indices), summary statistics can be applied with respect to any of the index levels:
End of explanation
"""


mb.to_csv("mb.csv")

"""
Explanation: Writing Data to Files
As well as being able to read several data input formats, Pandas can also export data to a variety of storage formats. We will bring your attention to just a couple of these.
End of explanation
"""


baseball.to_pickle("baseball_pickle")

"""
Explanation: The to_csv method writes a DataFrame to a comma-separated values (csv) file. You can specify custom delimiters (via sep argument), how missing values are written (via na_rep argument), whether the index is writen (via index argument), whether the header is included (via header argument), among other options.
An efficient way of storing data to disk is in binary format. Pandas supports this using Python’s built-in pickle serialization.
End of explanation
"""


pd.read_pickle("baseball_pickle")

"""
Explanation: The complement to to_pickle is the read_pickle function, which restores the pickle to a DataFrame or Series:
End of explanation
"""


