diff --git a/README.md b/README.md
index 9eba85c..9b7bc09 100644
--- a/README.md
+++ b/README.md
@@ -6,9 +6,18 @@ This package contains tools for causal analysis using observational (rather than
 
 Assuming you have pip installed, just run
 ```
-pip install causality 
+pip install causality
 ```
 
+## [Causal Analysis](https://github.com/akelleh/causality/tree/master/causality/analysis)
+
+The simplest interface to this package is probably through the `CausalDataFrame` object in [`causality.analysis.CausalDataFrame`](https://github.com/akelleh/causality/blob/master/causality/analysis/dataframe.py#L8). This is just an extension of the `pandas.DataFrame` object, and so it inherits the same methods.
+
+The `CausalDataFrame` current supports two kinds of causal analysis. First, it has a `CausalDataFrame.zmean` method. This method lets you control for a set of variables, `z`, when you're trying to estimate the effect of a discrete variable `x` on a continuous variable, `y`. It supports both returning the `y` estimates at each `x` value, as well as providing bootstrap error bars. For more details, check out the readme [here](https://github.com/akelleh/causality/tree/master/causality/analysis).
+
+The second kind of analysis supported is plotting to show the effect of discrete or continuous `x` on continous `y` while controlling for `z`. You can do this with the `CausalDataFrame.zplot` method. For details, check out the readme [here](https://github.com/akelleh/causality/tree/master/causality/analysis).
+
+
 ## Measuring Causal Effects
 
 the [`causality.estimation`](https://github.com/akelleh/causality/tree/master/causality/estimation) module contains tools for estimating causal effects from observational and experimental data. Most tools are parametric, like `PropensityScoreMatching`, and can be found in `causality.estimation.parametric`. Other models are non-parametric, and rely on directly estimating densities and using the g-estimation approach.
@@ -51,10 +60,10 @@ graph = ic_algorithm.search(X, variable_types)
 Now, we have the inferred graph stored in `graph`.  In this graph, each variable is a node (named from the DataFrame columns), and each edge represents statistical dependence between the nodes that can't be eliminated by conditioning on the variables specified for the search.  If an edge can be oriented with the data available, the arrowhead is indicated in `'arrows'`.  If the edge also satisfies the local criterion for genuine causation, then that directed edge will have `marked=True`.  If we print the edges from the result of our search, we can see which edges are oriented, and which satisfy the local criterion for genuine causation:
 ```python
 >>> graph.edges(data=True)
-[('x2', 'x1', {'arrows': [], 'marked': False}), 
- ('x2', 'x4', {'arrows': ['x4'], 'marked': False}), 
- ('x3', 'x1', {'arrows': [], 'marked': False}), 
- ('x3', 'x4', {'arrows': ['x4'], 'marked': False}), 
+[('x2', 'x1', {'arrows': [], 'marked': False}),
+ ('x2', 'x4', {'arrows': ['x4'], 'marked': False}),
+ ('x3', 'x1', {'arrows': [], 'marked': False}),
+ ('x3', 'x4', {'arrows': ['x4'], 'marked': False}),
  ('x4', 'x5', {'arrows': ['x5'], 'marked': True})]
 ```
 
diff --git a/causality/analysis/README.md b/causality/analysis/README.md
new file mode 100644
index 0000000..36eecc8
--- /dev/null
+++ b/causality/analysis/README.md
@@ -0,0 +1,153 @@
+# Analysis
+
+This module contains tools for using the Robin's G-Formula and arbitrary machine learning estimators to estimate and plot causal effects. By "causal effect" we mean the distribution or conditional expectation of Y given X, controlling for an admissable set of covariates, Z, to make the effect identifiable. For a primer on choosing these Z variables, check out the article [here](https://medium.com/@akelleh/a-technical-primer-on-causality-181db2575e41).
+
+More intuitively, you want to estimate the effect of X on Y, but you know you need to control for some set of confounders, Z, to get the true effect. Otherwise, you expect there to be confounding bias.
+
+# The `CausalDataFrame`
+
+The `CausalDataFrame` is an extension of the `pandas.DataFrame`, so you can intialize it as you normally would intialize a `pandas.DataFrame`, e.g.
+
+```python
+import numpy as np
+import pandas as pd
+import matplotlib.pyplot as pp
+from causality.analysis.dataframe import CausalDataFrame
+
+N = 1000
+
+z = np.random.normal(1., size=N)
+x = np.random.binomial(1, p=1./(1. + np.exp(-z/.1)))
+y = x + z + np.random.normal(size=N)
+
+# It's easy to create a data frame
+df = CausalDataFrame({'x': x, 'y': y, 'z': z})
+```
+Here, we've created a dataset where `x` has a direct effect on `y`, but a third variable `z` has a direct effect on both `x` and `y`. The result will be that when `z` is higher, both `x` and `y` will tend to be higher due to the influence of `z`. This correlation is not due to the causal dependence between `x` and `y`, and is instead due to confounding by `z`. We can make a causal plot that controls for `z` so we can see the true dependence between `x` and `y` easily with the `CausalDataFrame`:
+
+
+```python
+# and the interface to zplot is basically the same as the pandas.DataFrame.plot method!
+df.zplot(x='x', y='y', z_types={'z': 'c'}, z=['z'], kind='bar', bootstrap_samples=500); pp.ylabel("$E[Y|do(X=x)]$"); pp.show()
+
+```
+![The causal estimate](./img/discrete_zplot.png)
+
+This `zplot` method passes args and kwargs to the `plot` method of the `pandas.DataFrame`, so you can pass it all of its usual formatting options. We'll give a more complete summary of all of its particular methods below.
+
+ You can also still use all of the usual `DataFrame` methods, for example to get a naive plot for comparison.
+
+ ```python
+df.groupby('x').mean().reset_index().plot(x='x', y='y', kind='bar'); pp.ylabel("$E[Y|X=x]$"); pp.show()
+ ```
+ ![The naive estimate](./img/discrete_zplot_naive.png)
+
+The correct answer in this example is that if you intervene to set the value of `x` to `x=0`, you'll find (on average) `y=1`. If you set `x=1`, you'll find (on average) `y=2`. You can see the causal `zplot` method finds the correct answer, within the 95% confidence level. You can see naive observational estimate has much lower `y` at `x=0`!
+
+## The `CausalDataFrame.zplot` method
+
+If you have continous `x`, you can instead use the `kind='line'` argument. Let's generate a similar data set to see how this works. This time, let's have two confounding variables.
+
+```python
+N = 1000
+lower = -1
+upper = 1
+z1 = np.random.uniform(lower, upper, size=N)
+z2 = np.random.uniform(lower, upper, size=N)
+x  = np.random.uniform(lower, upper, size=N) + (z1 + z2)/2.
+z  = z1 + z2
+y  = np.random.normal(size=N) - x + 2.* z
+X  = CausalDataFrame({'x': x, 'y': y, 'z1': z1, 'z2': z2})
+```
+We can see from this data generating process that the true relationship, holding the confounders `z` constant, is a negative relationship between `x` and `y` that is linear with a slope of `-1`. If we just plot the (confounded) data, we can't see any relationship at all:
+
+```python
+X.plot(x='x', y='y', style='bo', alpha=0.2, kind='scatter')
+```
+ ![The naive estimate](img/continuous_zplot_naive.png)
+
+We can control for the `z` variables, and recover the negative relationship!
+
+```python
+X.zplot(x='x', y='y', z=['z1', 'z2'], z_types={'z1': 'c', 'z2': 'c'}, kind='line')
+```
+ ![The random forest estimate](img/continuous_zplot_random_forest.png)
+
+ Unfortunately, the relationship is very noisy. The model used by default to do the controlling is a random forest model. It won't be the best model for every problem, and doesn't work here as well as kernel regression. Those are the two typed models that are currently supported for automatic controlling. You can switch to kernel density regression by specifying `model_type='kernel'`.
+
+```python
+X.zplot(x='x', y='y', z=['z1', 'z2'], z_types={'z1': 'c', 'z2': 'c'}, kind='line', model_type='kernel')
+```
+ ![The kernel estimate](img/continuous_zplot_kernel.png)
+
+ You're free to use other models you define yourself, as well. The models can be fitted or not. If the model is not fitted, you should pass the model object through the `model` kwarg.
+
+ ```python
+ from sklearn.linear_model import LinearRegression
+
+treatment = 'x'
+outcome = 'y'
+confounders = ['z1', 'z2']
+X.zplot(x='x', y='y', z=confounders, z_types={'z1': 'c', 'z2': 'c'}, kind='line', model=LinearRegression)
+```
+![linear regression model results](./img/continuous_zplot_linear.png)
+
+
+If you'd like to pass a fit model, for example as you might if you're fitting a complicated model like a neural network, you can pass it through the `fitted_model` kwarg. Here's a simple multi-layer perceptron, just to give an example.
+
+```python
+from sklearn.neural_network import MLPRegressor
+model = MLPRegressor(hidden_layer_sizes=(128,128,128), max_iter=100, learning_rate_init=0.01)
+
+treatment = 'x'
+outcome = 'y'
+confounders = ['z1', 'z2']
+model.fit(X[[treatment] + confounders], X[outcome])
+# requirement: model.predict(X[[treatment] + confounders]) yields a numpy array of scalar predictions for y, dimension (n_samples,)
+X.zplot(x='x', y='y', z=confounders, z_types={'z1': 'c', 'z2': 'c'}, kind='line', fitted_model=model)
+```
+
+![MLP model results](./img/continuous_zplot_mlp.png)
+
+With discrete variables, it's a little easier to bootstrap error bars. We have some kwargs available to give a little extra control over the bootstrapping process. Note that we use the normal approximation for the bootstrap confidence intervals. Using the percentile approach tended to give overly narrow intervals.
+
+Returning to the discrete data example from before,
+```python
+N = 1000
+
+z = np.random.normal(1., size=N)
+x = np.random.binomial(1, p=1./(1. + np.exp(-z/.1)))
+y = x + z + np.random.normal(size=N)
+
+# It's easy to create a data frame
+df = CausalDataFrame({'x': x, 'y': y, 'z': z})
+
+# and the interface to zplot is basically the same as the pandas.DataFrame.plot method!
+df.zplot(x='x', y='y', z_types={'z': 'c'}, z=['z'], kind='bar', bootstrap_samples=500); pp.ylabel("$E[Y|do(X=x)]$"); pp.show()
+```
+![discrete plot with bootstrap parameter](./img/discrete_zplot_bootstrap.png)
+
+The default number of samples for the bootstrap is 500 samples. Much less than that tended to give overly narrow intervals. I'd encourage you to test them yourself with a simulation if you're planning to vary this parameter!
+
+You can also adjust the confidence level for your error bars. The default is the 95% confidence level. As we decrease the confidence level, you'll see the error bars shrink.
+```
+df.zplot(x='x', y='y', z=['z'], z_types={'z': 'c'}, kind='bar', bootstrap_samples=500, confidence_level=0.80)
+```
+![discrete plot with bootstrap parameter and 80%CL](./img/discrete_zplot_bootstrap_80CL.png)
+
+You might also like to just get the values and error bars out from these plots. You can do that for the discrete plot with the `zmean` method.
+
+## The `CausalDataFrame.zmean` method
+
+The interface for the `zmean` method is exactly like the `zplot` method. You can pass models and bootstrap parameters in the same way. The return value is a new dataframe.
+
+```python
+df.zmean(x='x', y='y', z=['z'], z_types={'z': 'c'}, bootstrap_samples=500, confidence_level=0.95)
+```
+![zmean results](img/zmean_results.png)
+
+Leaving off the `bootstrap_samples` kwarg while specifying the `confidence_level` will cause `zmean` to default to `bootstrap_samples=500`. Leaving off both keyword arguments will result in no confidence intervals being given:
+```python
+df.zmean(x='x', y='y', z=['z'], z_types={'z': 'c'})
+```
+![zmean results](img/zmean_results_no_CI.png)
diff --git a/causality/analysis/img/continuous_zplot_kernel.png b/causality/analysis/img/continuous_zplot_kernel.png
new file mode 100644
index 0000000..c30b13e
Binary files /dev/null and b/causality/analysis/img/continuous_zplot_kernel.png differ
diff --git a/causality/analysis/img/continuous_zplot_linear.png b/causality/analysis/img/continuous_zplot_linear.png
new file mode 100644
index 0000000..f7971bc
Binary files /dev/null and b/causality/analysis/img/continuous_zplot_linear.png differ
diff --git a/causality/analysis/img/continuous_zplot_mlp.png b/causality/analysis/img/continuous_zplot_mlp.png
new file mode 100644
index 0000000..e4ae808
Binary files /dev/null and b/causality/analysis/img/continuous_zplot_mlp.png differ
diff --git a/causality/analysis/img/continuous_zplot_naive.png b/causality/analysis/img/continuous_zplot_naive.png
new file mode 100644
index 0000000..0b133ec
Binary files /dev/null and b/causality/analysis/img/continuous_zplot_naive.png differ
diff --git a/causality/analysis/img/continuous_zplot_random_forest.png b/causality/analysis/img/continuous_zplot_random_forest.png
new file mode 100644
index 0000000..aaab03a
Binary files /dev/null and b/causality/analysis/img/continuous_zplot_random_forest.png differ
diff --git a/causality/analysis/img/discrete_zplot.png b/causality/analysis/img/discrete_zplot.png
new file mode 100644
index 0000000..22ea596
Binary files /dev/null and b/causality/analysis/img/discrete_zplot.png differ
diff --git a/causality/analysis/img/discrete_zplot_bootstrap.png b/causality/analysis/img/discrete_zplot_bootstrap.png
new file mode 100644
index 0000000..22ea596
Binary files /dev/null and b/causality/analysis/img/discrete_zplot_bootstrap.png differ
diff --git a/causality/analysis/img/discrete_zplot_bootstrap_80CL.png b/causality/analysis/img/discrete_zplot_bootstrap_80CL.png
new file mode 100644
index 0000000..46f2b93
Binary files /dev/null and b/causality/analysis/img/discrete_zplot_bootstrap_80CL.png differ
diff --git a/causality/analysis/img/discrete_zplot_naive.png b/causality/analysis/img/discrete_zplot_naive.png
new file mode 100644
index 0000000..2e49980
Binary files /dev/null and b/causality/analysis/img/discrete_zplot_naive.png differ
diff --git a/causality/analysis/img/zmean_results.png b/causality/analysis/img/zmean_results.png
new file mode 100644
index 0000000..2d2b4aa
Binary files /dev/null and b/causality/analysis/img/zmean_results.png differ
diff --git a/causality/analysis/img/zmean_results_no_CI.png b/causality/analysis/img/zmean_results_no_CI.png
new file mode 100644
index 0000000..df30613
Binary files /dev/null and b/causality/analysis/img/zmean_results_no_CI.png differ