Add function to compute R(n, epsilon) in METE
Add function to compute P(n, epsilon) in SSNT
Add function to plot heat map for individual-level joint P(N, m)
Fix minor bug