Show years for which the circulation data are missing
. tabulate circ year if circ == ., missing
Create icirc containing a linear interpolation of circ on year for missing values of
circ and perform this calculation separately for each magazine
. by magazine: ipolate circ year, gen(icirc)