winter = c(1:60, 305:365)
wooster.winter = wooster
wooster.winter$data =  - wooster$data[c(winter, winter + 365,
winter + 365 * 2, winter + 365 *3),]
wooster.winter$data[,1]= - wooster.winter$data[,1]
wooster.winter$data <- cbind(wooster.winter$data,.year.number=c(rep(1:4,rep(121,4))))