par 30000
*
* This file takes msudata from mcom and generates vector of trends
*
read(msudata.dif) m1-m218 / dif
*
*                 There are 218 cols of 264 monthly observations
*
*                 Make time trend
sample 1 264
genr t=time(0)
*
*
*                 Make decadal trend vector
dim mtrend 218
*
*                 do trend calculation col by col
set nodoecho
do #=1,218
   ?genr m#=120*m#/100
   ?ols m# t /coef=beta
*
*                 write trend coef to mtrend vector
   sample # #
   genr mtrend=beta:1
   sample 1 264
*
   endo 
*
sample 1 218
write(mt1) t mtrend
stop