sas中怎样做probit回归

data salary;
input yr inc : @@;
format yr year4.;
label yr="年份" inc="工资";
cards;
1985 1383.0 1986 1659.7 1987 1881.1 1988 2316.2 1989 2618.5 1990 2951.1 1991 3323.9 1992 3939.2 1993 4916.2 1994 6656.4 1995 8255.8 1996 9249.4 1997 9602.4 1998 9540.2 1999 10155.9 2000 10954.7 2001 12205.4 2002 13638.1 2003 15329.6 2004 17615.0 2005 20627.1 2006 24262.3 2007 29471.5 2008 35289.5 2009 40288.2
run;
proc corr data=salary;
by yr;
run;
proc reg data=salary;
model inc=yr/p;
run;
quit;