add interval module for interval estimation + 1 example pllr_interval.py, document all of that