Part 1.
      DSTAT ; rhs = * $
      CREATE ; loginc=log(income)$
      CALC ; xbar=xbr(loginc);sdv=sdv(loginc)$
      CREATE ; normal=rnn(xbar,sdv)$
      KERNEL;rhs=loginc,normal$

      NAMELIST ; demogrfc = age, female, married $
      NAMELIST ; years = year1984,year1985,year1986,year1987,year1988,year1991 $

Part 2.
      REGRESS ; lhs = income ; rhs = one, demogrfc, years $
      CALC ; rsq0 = rsqrd $
      REGRESS ; lhs = income ; rhs = one, demogrfc, years, educ $
      CALC ; rsq1 = rsqrd $
      CALC ; list ; fstat = ((rsq1 - rsq0)/1) / ((1-rsq1)/(n-kreg)) $

Part 3.
      SAMPLE ; 1 - 100 $  (use the first 100 observations)
      REGRESS ; lhs = income ; rhs = one,demogrfc $
      REGRESS ; lhs = income ; rhs = one,demogrfc ; alg=LAD ; nbt=25 $
      REGRESS ; lhs = income ; rhs = one,demogrfc ; alg=LAD ; nbt=25 $
      CALC ; ran (12345) $ (use a 5 digit odd number)
      REGRESS ; lhs = income ; rhs = one,demogrfc ; alg=LAD ; nbt=25 $
      CALC ; ran ( 12345) $ (use a 5 digit odd number)=
      REGRESS ; lhs = income ; rhs = one,demogrfc ; alg=LAD ; nbt=25 $

Part 4.
      SAMPLE  ;all $
      REGRESS ; lhs = income ; rhs = one,demogrfc,years$
      REGRESS ; lhs = income ; rhs = one,demogrfc,years ; cluster=id $

Part 5. 
      SAMPLE  ; all $
      REGRESS ; lhs = income ; rhs = one,demogrfc,years ;cluster=id
              ; test: years = 0 $
      REGRESS ; lhs = income ; rhs = one,demogrfc,years ; cluster=id$
      MATRIX  ; by=b(5:10) ; vy=varb(5:10,5:10) $
      MATRIX  ; list ; wld = by'<vy>by $

Part 6.
      REGRESS ; Lhs = loginc ; rhs = one,age, educ, female, age*educ,
            age^2, age*female, educ*female $
      PARTIAL ; effects: age | educ = 12,14,16,18,20 $
      PARTIAL ; effects: educ & age = 25(3)64 ; plot(ci) $
      PARTIAL ; effects: female | educ = 12,16,20 & age = 25(5)64 ; plot $
      PARTIAL ; simulate: female | educ = 12,16,20 & age = 25(5)64 ; plot $

Part 7
      CREATE ; v2 = loginc^2 ; v3 = v2*loginc ; v4 = v3*loginc $
      CALC      ; sv = sqr (xbr (v2)) ;  m3 = xbr(v3) ; m4 = xbr (v4) $
      CALC      ; List ; Chisq = N *( ( m3/sv^3)^2 / 6  +  ((m4/sv^4) - 3)^2 /24 )  $