Skip to content

Commit bfd920a

Browse files
committed
Added missing normalization step
1 parent eef7954 commit bfd920a

File tree

1 file changed

+8
-7
lines changed

1 file changed

+8
-7
lines changed

structure_threader/wrappers/maverick_wrapper.py

Lines changed: 8 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -182,12 +182,13 @@ def maverick_normalization(x_mean, x_sd, klist, draws=1e6, limit=95):
182182
z = np.zeros([len(x_mean), draws])
183183

184184
for i in range(z.shape[0]):
185-
z[i] = np.sort(np.array(
186-
[np.exp(rnorm(x_mean[i], x_sd[i])) for _ in range(draws)]))
185+
y = np.array([np.exp(rnorm(x_mean[i], x_sd[i])) for _ in range(draws)])
186+
z[i] = np.sort(y/sum(y))
187187

188-
norm_res = dict((i, {"norm_mean": np.mean(z[i]),
189-
"lower_limit": np.percentile(z[i], 2.5),
190-
"upper_limit": np.percentile(z[i], 97.5)}) for i, k in
191-
enumerate(klist))
188+
norm_res = dict(
189+
(i, {"norm_mean": np.mean(z[i]),
190+
"lower_limit": np.percentile(z[i], 2.5),
191+
"upper_limit": np.percentile(z[i], 97.5)})
192+
for i, k in enumerate(klist))
192193

193-
return norm_res
194+
return norm_res

0 commit comments

Comments
 (0)