In [1]:

    

%run src/ssy_discretized_test.py


    

In [2]:

    

import pandas as pd
import matplotlib.pyplot as plt


    

In [3]:

    

ssy = SSY()


    

In [4]:

    

D=3


    

In [5]:

    

test_val_spec_rad(ssy, K=D, I=D, J=D)


    

    
Out[5]:





0.9995003158424812

In [6]:

    

f = mc_factory(ssy, K=D, I=D, J=D)


    

In [7]:

    

n_vals = np.array([750, 1000])
m_vals = np.array([6000, 8000, 10000, 12000, 14000])


    

In [8]:

    

k = 1000
draws = np.empty(k)
means = np.empty((len(n_vals), len(m_vals)))
stds = np.empty((len(n_vals), len(m_vals)))

for n_i, n in enumerate(n_vals):
    for m_i, m in enumerate(m_vals):
        print(f'Calculating n={n}, m={m}')
        for i in range(k):
            draws[i] = f(n=n, m=m)
        means[n_i, m_i] = draws.mean()
        stds[n_i, m_i] = draws.std()


    

    




Calculating n=750, m=6000
Calculating n=750, m=8000
Calculating n=750, m=10000
Calculating n=750, m=12000
Calculating n=750, m=14000
Calculating n=1000, m=6000
Calculating n=1000, m=8000
Calculating n=1000, m=10000
Calculating n=1000, m=12000
Calculating n=1000, m=14000

In [9]:

    

means_strings = means.round(6).astype(str)
n_strings = n_vals.astype(str)

start_table = r"""
\begin{table}
\centering
\begin{tabular}{llll}
"""

m_table = ' & m = '.join(m_vals.astype(str))
m_table = '    & m = ' + m_table + r' \\' + '\n' + r'\hline \hline' '\n'

end_table = r"""
\end{tabular}
\end{table}
"""

row_string = ''
for row in range(len(n_strings)):
    
    temp_means = ' & '.join(means_strings[row, :])
    
    x = ['{:f}'.format(item) for item in stds[row, :]]    
    temp_stds = '(' + ') & ('.join(x) + ')'
    row_string += f'n = {n_strings[row]} & ' + temp_means + r' \\' + '\n'
    row_string += '         & ' + temp_stds + r' \\' + '\n'
    row_string += r'\hline' '\n'

print(start_table + m_table + row_string + end_table)


    

    




\begin{table}
\centering
\begin{tabular}{llll}
    & m = 6000 & m = 8000 & m = 10000 & m = 12000 & m = 14000 \\
\hline \hline
n = 750 & 0.999551 & 0.999551 & 0.999552 & 0.999552 & 0.99955 \\
         & (0.000018) & (0.000017) & (0.000014) & (0.000014) & (0.000012) \\
\hline
n = 1000 & 0.999535 & 0.999536 & 0.999537 & 0.999538 & 0.999539 \\
         & (0.000018) & (0.000019) & (0.000017) & (0.000017) & (0.000018) \\
\hline

\end{tabular}
\end{table}

In [ ]:

    




    

In [ ]: