nulldata 100
genr x = uniform()
loop 500 --progressive 
  genr sdu = 2 + 5*x
  genr u = normal() * sdu
  genr y = 10 + 20*x + u
  ols y 0 x
  hccm y 0 x
  hsk y 0 x
  genr w = 1.0/sdu
  wls w y 0 x
endloop