version with:
1. phtest with fixed threshold on different user-session data
2. comparison with pelt and cpm result from R
3. show the change-point locations with moving average
4. cusum func from timeseries package