coxtv function

fit a Cox non-proportional hazards model