* Stats 506
* Josh Errickson
* 9/19/23
* This is the *output* from the Stata where I ran these commands
* interactively, before moving over to the do-file.

. mata: 3
  3

. mata:
------------------------------------------------- mata (type end to
> exit) ------------------------------------------------------------
: end
--------------------------------------------------------------------

. mata:
------------------------------------------------- mata (type end to
> exit) ------------------------------------------------------------
: 5 -4
  1

: sqrt(2)
  1.414213562

: x =  4

: x
  4

: x + 2
  6

: y = 3

: x*y
  12

: M = (1 , 2/ 4, 7)

: M
        1    2    3
    +----------------+
  1 |   1   .5    7  |
    +----------------+

: M = (1 , 2 \ 4, 7)

: M
       1   2
    +---------+
  1 |  1   2  |
  2 |  4   7  |
    +---------+

: N = (1, 4, 6 \ 2, 0, 4)

: M
       1   2
    +---------+
  1 |  1   2  |
  2 |  4   7  |
    +---------+

: N
       1   2   3
    +-------------+
  1 |  1   4   6  |
  2 |  2   0   4  |
    +-------------+

: M*N
        1    2    3
    +----------------+
  1 |   5    4   14  |
  2 |  18   16   52  |
    +----------------+

: M - N
                 <istmt>:  3200  conformability error
--------------------------------------------------------------------
r(3200);

. mata:
------------------------------------------------- mata (type end to
> exit) ------------------------------------------------------------
: M
       1   2
    +---------+
  1 |  1   2  |
  2 |  4   7  |
    +---------+

: N
       1   2   3
    +-------------+
  1 |  1   4   6  |
  2 |  2   0   4  |
    +-------------+

: N'
       1   2
    +---------+
  1 |  1   2  |
  2 |  4   0  |
  3 |  6   4  |
    +---------+

: Y = M*N

: Y
        1    2    3
    +----------------+
  1 |   5    4   14  |
  2 |  18   16   52  |
    +----------------+

: M
       1   2
    +---------+
  1 |  1   2  |
  2 |  4   7  |
    +---------+

: N
       1   2   3
    +-------------+
  1 |  1   4   6  |
  2 |  2   0   4  |
    +-------------+

: M , N
       1   2   3   4   5
    +---------------------+
  1 |  1   2   1   4   6  |
  2 |  4   7   2   0   4  |
    +---------------------+

: M , N \ M , N
       1   2   3   4   5
    +---------------------+
  1 |  1   2   1   4   6  |
  2 |  4   7   2   0   4  |
  3 |  1   2   1   4   6  |
  4 |  4   7   2   0   4  |
    +---------------------+

: I(6)
[symmetric]
       1   2   3   4   5   6
    +-------------------------+
  1 |  1                      |
  2 |  0   1                  |
  3 |  0   0   1              |
  4 |  0   0   0   1          |
  5 |  0   0   0   0   1      |
  6 |  0   0   0   0   0   1  |
    +-------------------------+

: J(4, 6, -7)
        1    2    3    4    5    6
    +-------------------------------+
  1 |  -7   -7   -7   -7   -7   -7  |
  2 |  -7   -7   -7   -7   -7   -7  |
  3 |  -7   -7   -7   -7   -7   -7  |
  4 |  -7   -7   -7   -7   -7   -7  |
    +-------------------------------+

: J(10, 1, 1)
        1
     +-----+
   1 |  1  |
   2 |  1  |
   3 |  1  |
   4 |  1  |
   5 |  1  |
   6 |  1  |
   7 |  1  |
   8 |  1  |
   9 |  1  |
  10 |  1  |
     +-----+

: X
                 <istmt>:  3499  X not found
--------------------------------------------------------------------
r(3499);

. mata:
------------------------------------------------- mata (type end to
> exit) ------------------------------------------------------------
: M
       1   2
    +---------+
  1 |  1   2  |
  2 |  4   7  |
    +---------+

: diag(M)
[symmetric]
       1   2
    +---------+
  1 |  1      |
  2 |  0   7  |
    +---------+

: diagonal(M)
       1
    +-----+
  1 |  1  |
  2 |  7  |
    +-----+

: diagonal(M)'
       1   2
    +---------+
  1 |  1   7  |
    +---------+

: end
--------------------------------------------------------------------

.