Convert Stata matrices to numpy arrays

mat2nparray.ado

capture program drop mat2nparray
program define mat2nparray 
    version 11.0
    syntax namelist(min=1), SAVing(str) [ Format(str) APPend REPlace ]
    if "`format'"=="" local format "%16.0g"
    local saving: subinstr local saving "." ".", count(local ext)
    if !`ext' local saving "`saving'.py"
    tempname myfile
    file open `myfile' using "`saving'", write text `append' `replace'
    file write `myfile' "import numpy as np" _n _n

    foreach mat of local namelist {
        mkarray `mat' `myfile' `format'
    }
    file write `myfile' "class Bunch(dict):" _n
    file write `myfile' "    def __init__(self, **kw):" _n
    file write `myfile' "        dict.__init__(self, kw)" _n
    file write `myfile' "        self.__dict__  = self" _n _n _n
    file write `myfile' "results = Bunch("
    foreach mat of local namelist {
        file write `myfile' "`mat'=`mat', "
    }
    file write `myfile' ")" _n _n
file close `myfile'
end

capture program drop mkarray
program define mkarray

    args mat myfile fmt
    local nrows = rowsof(`mat')
    local ncols = colsof(`mat')
    local i 1
    local j 1
    file write `myfile' "`mat' = np.array(["
    local justifyn = length("`mat' = np.array([")
    forvalues i=1/`nrows' {
        forvalues j = 1/`ncols' {
            if `i' > 1 | `j' > 1 { // then we need to indent
                forvalues k=1/`justifyn' {
                    file write `myfile' " "
                }
            }
            if `i' < `nrows' | `j' < `ncols' {
                if `mat'[`i',`j'] == . {
                    file write `myfile' "np.nan" ", " _n
                }
                else {
                    file write `myfile' `fmt' (`mat'[`i',`j']) ", " _n
                }
            }
            else {
                if `mat'[`i',`j'] == . {
                    file write `myfile' "np.nan"
                }
                else {
                    file write `myfile' `fmt' (`mat'[`i',`j'])
                }
            }
        }
    } 

    if `nrows' == 1 | `ncols' == 1 {
        file write `myfile' "])" _n _n
    }
    else {
        file write `myfile' "]).reshape(`nrows',`ncols')" _n _n
    }
end