deletions | additions       diff --git a/untitled.tex b/untitled.tex  index e0560c9..f1c6ff2 100644  --- a/untitled.tex  +++ b/untitled.tex 
...
\section{Methods}  \subsection{Data simulation}    We used fastsimcoal2 \cite{Excoffier_2013} to simulate ancestral trees and $3000$ haplotypes of $4000$ equispaced SNVs in a $2$ million base pair(Mbp) genomic region with recombination rate $1\times10^{-8}$per base pair per generation in a haploid population of constant effective size, $N_e=6200$ \cite{tenesa2007recent}. \cite{Tenesa_2007}.  We then randomly paired these 3000 haplotypes into 1500 individuals.   Once we randomly paired haplotypes into 1500 individuals, the disease status was assigned to the individuals based on randomly sampled susceptibility SNVs from the middle region of $950kbp-1050kbp$. For the disease susceptibility variants, the number of copies of the derived allele increases disease risk according to a logistic regression model: