Skip to content
Snippets Groups Projects
Commit 394ff5d1 authored by Martin Renoult's avatar Martin Renoult
Browse files

Update Kalman_filter

parent 91b06d49
No related branches found
No related tags found
No related merge requests found
...@@ -94,10 +94,8 @@ stdprior = np.sqrt(np.diag(samp_var)) ...@@ -94,10 +94,8 @@ stdprior = np.sqrt(np.diag(samp_var))
gausspri = np.random.normal(loc=pri_mn[0], scale=stdprior[0], size=1000) gausspri = np.random.normal(loc=pri_mn[0], scale=stdprior[0], size=1000)
gausspost = np.random.normal(loc=pos_mn[0], scale=stdfinal[0], size=1000) gausspost = np.random.normal(loc=pos_mn[0], scale=stdfinal[0], size=1000)
post_stats_90_pri = stats.quantiles(gausspri, qlist=(5, 95)) post_stats_90_pri = np.percentile(gausspri, q=(5, 95))
post_stats_90_pri = [ v for v in post_stats_90_pri.values()] post_stats_90_post = np.percentile(gausspost, q=(5, 95))
post_stats_90_post = stats.quantiles(gausspost, qlist=(5, 95))
post_stats_90_post = [ v for v in post_stats_90_post.values()]
plt.xlabel('Climate Sensitivity') plt.xlabel('Climate Sensitivity')
plt.ylabel('LGM tropical (20S - 30N) temperature change') plt.ylabel('LGM tropical (20S - 30N) temperature change')
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment