DATA quakes;

INFILE 'quake.dat';
INPUT year n_over_7;

proc means;

proc genmod	; model	n_over_7 = year	/
   LINK	= identity dist	= Poisson;		/* additive	model=>	difference in average COUNT	*/

proc genmod	; model	n_over_7 = year	/
   LINK	= identity dist	= Normal;		/* additive	model=>	difference in average COUNT	*/
										/* but asking that SD's	be estimated separately	*/

proc genmod	; model	n_over_7 = year	/
   LINK	= log	   dist	= Poisson;		/* multiplicative model= ratio of average COUNTS */

proc genmod	; model	n_over_7 = year	/
   LINK	= log	   dist	= Normal;		/* multiplicative model=> ratio	of average COUNTS */
										/* but asking that SD's	be estimated separately	*/


RUN;