Mata
Y = (2\ 1\ 1\ 1\ 5\ 4\ 7\ 6\ 7\ 8\ 3\ 3\ 6\ 6\ 10\ 9\ 6\ 6\ 9\ 10)
X = (2,4\2,4\1,4\1,3\3,6\4,6\5,3\5,4\7,3\6,3\4,5\3,5\6,9\6,8\8,6\9,7\10,5\9,5\4,7\4,7)
Stata
mat Y = (2\ 1\ 1\ 1\ 5\ 4\ 7\ 6\ 7\ 8\ 3\ 3\ 6\ 6\ 10\ 9\ 6\ 6\ 9\ 10)
mat X = (2,4\2,4\1,4\1,3\3,6\4,6\5,3\5,4\7,3\6,3\4,5\3,5\6,9\6,8\8,6\9,7\10,5\9,5\4,7\4,7)
SAS IML
Y = {2, 1, 1, 1, 5, 4, 7, 6, 7, 8, 3, 3, 6, 6, 10, 9, 6, 6, 9, 10};
X = {
 2 4,
 2 4,
 1 4,
 1 3,
 3 6,
 4 6,
 5 3,
 5 4,
 7 3,
 6 3,
 4 5,
 3 5,
 6 9,
 6 8,
 8 6,
 9 7,
10 5,
 9 5,
 4 7,
 4 7};