discontinuity code snippets in r

create a regression discontinuity model in r