capture clear
input byte y  byte x
4 -2
1 -1
0  0
1  1
4  2
end

list, noobs


/* Chapter 6, problem 2a: Linear fit */
quietly {
  regress y x
  local b0=_b[_cons]
  local b1=_b[x]
  #delimit ;
  graph twoway (scatter y x, 
                  msymbol(circle)       /* Change shape of markers     */
                  msize(large)          /* Make markers larger         */
                  mlwidth(medium)       /* Change marker outline width */
                  mlcolor(maroon)       /* Change marker outline color */
                  mfcolor(ebblue))      /* Change marker fill color    */
               (lfit    y x, lw(thick)
    text(3.6 -1.6 "Y= `b0' + `b1'(X) + e"))
  ;
  # delimit cr
}

/* Chapter 6, problem 2d: Polynomial fit */
  quietly {
  regress y x
  local b0=_b[_cons]
  local b1=_b[x]
  #delimit ;
  graph twoway (scatter y x, 
                  msymbol(circle)       /* Change shape of markers     */
                  msize(large)          /* Make markers larger         */
                  mlwidth(medium)       /* Change marker outline width */
                  mlcolor(maroon)       /* Change marker outline color */
                  mfcolor(ebblue))      /* Change marker fill color    */
               (qfit    y x, lw(thick)
    text(3.6 -1.6 "Y=X^2"))
  ;
  # delimit cr
}