add wensitan.R for computing weekly trend