diff --git a/notebooks/HiC-linksize-empirical-distribution.ipynb b/notebooks/HiC-linksize-empirical-distribution.ipynb
new file mode 100644
index 0000000..4842bf1
--- /dev/null
+++ b/notebooks/HiC-linksize-empirical-distribution.ipynb
@@ -0,0 +1,322 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# HiC link size empirical distributions for datasets in ALLHiC paper\n",
+ "\n",
+ "We plot the link size (x-axis) vs link density (y-axis) for each of the four datasets used in the ALLHiC paper. We then try to infer the coeffiencients in the power law model.\n",
+ "\n",
+ "$$Y=A \\times X^B$$\n",
+ "\n",
+ "Based on `allhic extract`, the coefficients for the four datasets are:\n",
+ "- Ler0 (AAGCTT): $Y=0.00162 \\times X^{-0.8135}$\n",
+ "- Rice (GATC): $Y=0.000347 \\times X^{-0.8042}$\n",
+ "- AP85 Chr1A (AAGCTT): $Y=1.62e^{-5} \\times X^{-0.8164}$\n",
+ "- Molokai Chr01A (AAGCTT): $Y=0.000325 \\times X^{-0.8844}$"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Dependencies"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "%matplotlib inline\n",
+ "import matplotlib.pyplot as plt\n",
+ "import seaborn as sns\n",
+ "sns.set_style(\"darkgrid\")\n",
+ "\n",
+ "import pandas as pd"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Load data"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 7,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " | \n",
+ " #Bin | \n",
+ " BinStart | \n",
+ " BinSize | \n",
+ " NumLinks | \n",
+ " TotalSize | \n",
+ " LinkDensity | \n",
+ "
\n",
+ " \n",
+ " \n",
+ " \n",
+ " | 0 | \n",
+ " 0 | \n",
+ " 2048 | \n",
+ " 90 | \n",
+ " 45072 | \n",
+ " 119653414 | \n",
+ " 0.000004 | \n",
+ "
\n",
+ " \n",
+ " | 1 | \n",
+ " 1 | \n",
+ " 2138 | \n",
+ " 95 | \n",
+ " 49670 | \n",
+ " 119652784 | \n",
+ " 0.000004 | \n",
+ "
\n",
+ " \n",
+ " | 2 | \n",
+ " 2 | \n",
+ " 2233 | \n",
+ " 99 | \n",
+ " 44248 | \n",
+ " 119652119 | \n",
+ " 0.000004 | \n",
+ "
\n",
+ " \n",
+ " | 3 | \n",
+ " 3 | \n",
+ " 2332 | \n",
+ " 103 | \n",
+ " 44274 | \n",
+ " 119651426 | \n",
+ " 0.000004 | \n",
+ "
\n",
+ " \n",
+ " | 4 | \n",
+ " 4 | \n",
+ " 2435 | \n",
+ " 108 | \n",
+ " 45417 | \n",
+ " 119650705 | \n",
+ " 0.000004 | \n",
+ "
\n",
+ " \n",
+ "
\n",
+ "