From f5e085b5c1a7c8c02c8eae7b584256cc2ea55aed Mon Sep 17 00:00:00 2001
From: Le Roux Erwan <erwan.le-roux@irstea.fr>
Date: Fri, 20 Sep 2019 16:26:58 +0200
Subject: [PATCH] [TREND ANALYSIS] add verbose argument when creating hypercube
 visualizer. add shape analysis from the MLE estimator, and the shapes are
 indeed between -0.4 and 0.4

---
 .../utils_hypercube.py                        |   7 +-
 .../shape_prior_check/__init__.py             |   0
 .../analyse_shape_from_some_experiment.py     |  10 +
 .../shape_prior_check/hist_values_shape.png   | Bin 0 -> 10172 bytes
 .../shape_from_some_experiment.txt            | 396 ++++++++++++++++++
 .../some_experiment_EVAN.py                   |   4 +-
 experiment/paper1_steps/utils.py              |   6 +-
 .../abstract_gev_trend_test.py                |   7 +
 8 files changed, 424 insertions(+), 6 deletions(-)
 create mode 100644 experiment/paper1_steps/poster_EVAN2019/shape_prior_check/__init__.py
 create mode 100644 experiment/paper1_steps/poster_EVAN2019/shape_prior_check/analyse_shape_from_some_experiment.py
 create mode 100644 experiment/paper1_steps/poster_EVAN2019/shape_prior_check/hist_values_shape.png
 create mode 100644 experiment/paper1_steps/poster_EVAN2019/shape_prior_check/shape_from_some_experiment.txt
 rename experiment/paper1_steps/poster_EVAN2019/{ => shape_prior_check}/some_experiment_EVAN.py (95%)

diff --git a/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/utils_hypercube.py b/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/utils_hypercube.py
index c756b0ab..267d4e35 100644
--- a/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/utils_hypercube.py
+++ b/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/utils_hypercube.py
@@ -32,11 +32,13 @@ def load_quantity_visualizer(quantity_hypercube_class, altitudes, last_starting_
 def load_altitude_visualizer(altitude_hypercube_class, altitudes, last_starting_year, nb_data_reduced_for_speed,
                              only_first_one, save_to_file, study_classes, trend_test_class
                              , exact_starting_year=None, first_starting_year=1958,
-                             orientations=None):
+                             orientations=None,
+                             verbose=True):
     visualizers = [StudyVisualizer(study, temporal_non_stationarity=True, verbose=False, multiprocessing=True)
                    for study in study_iterator_global(study_classes=study_classes, only_first_one=only_first_one,
                                                       altitudes=altitudes,
-                                                      orientations=orientations)]
+                                                      orientations=orientations,
+                                                      verbose=verbose)]
     altitude_to_visualizer = OrderedDict(zip(altitudes, visualizers))
     visualizer = altitude_hypercube_class(altitude_to_visualizer,
                                           save_to_file=save_to_file,
@@ -45,6 +47,7 @@ def load_altitude_visualizer(altitude_hypercube_class, altitudes, last_starting_
                                           last_starting_year=last_starting_year,
                                           first_starting_year=first_starting_year,
                                           exact_starting_year=exact_starting_year,
+                                            verbose=verbose,
                                           )
     assert isinstance(visualizer, AltitudeHypercubeVisualizer)
     return visualizer
diff --git a/experiment/paper1_steps/poster_EVAN2019/shape_prior_check/__init__.py b/experiment/paper1_steps/poster_EVAN2019/shape_prior_check/__init__.py
new file mode 100644
index 00000000..e69de29b
diff --git a/experiment/paper1_steps/poster_EVAN2019/shape_prior_check/analyse_shape_from_some_experiment.py b/experiment/paper1_steps/poster_EVAN2019/shape_prior_check/analyse_shape_from_some_experiment.py
new file mode 100644
index 00000000..97c42f9a
--- /dev/null
+++ b/experiment/paper1_steps/poster_EVAN2019/shape_prior_check/analyse_shape_from_some_experiment.py
@@ -0,0 +1,10 @@
+import matplotlib.pyplot as plt
+
+# Load shapes
+lines = []
+with open('shape_from_some_experiment.txt') as f:
+    for l in f:
+        lines.append(float(l.split('\n')[0]))
+# Build his
+plt.hist(lines, bins=50,  histtype='step')
+plt.show()
\ No newline at end of file
diff --git a/experiment/paper1_steps/poster_EVAN2019/shape_prior_check/hist_values_shape.png b/experiment/paper1_steps/poster_EVAN2019/shape_prior_check/hist_values_shape.png
new file mode 100644
index 0000000000000000000000000000000000000000..2aa54e3a4c6b777ae861ae7a2195f278cd138626
GIT binary patch
literal 10172
zcmeHtc{r5&-~WwhU!*Kq+o%+0C_CenNDE&{l29fj`@RmNbINX%J))>oM99vJP8j>1
zohf7=vScur=Y3b_Ip;d(JkRgBzSs4;e$O92*X5p>`<}V)`}6s{-|yG*z61ZIrLlF>
zo=pgXY}GuYb`e2X2?)aSaRUeJInmrZ1TX8HPHJA-0Ds;aZU(^TjrL~@oDgLD4fLPo
z1vb?ZcFH@SzT&KFXYTBF!_f?}x#4VYZRc!lX)5As=ICT;XDcEt@skARsK_m6XL|)n
z$v?gzVdrQexz#3g4M9W@O|_pdxhK&2eC_qVK3ZhzFW%qbcYtvB%l<3km(w;(wFFtF
z{p-T%afJe@oMQ^d6hds>nvBZAWE7{4<UJ}7G|W;w<96+GFP}k^R#sYZ&4#lA`j3Om
z+lp#;z1q*}r@eb;K)6Wnw6NgZiA{d0Q>E2-tdO)@rPowpOvmT&_*D9EiEd=|O(_;c
zRcd+23r6L54q-(OsQ>@XpXT7X32FUkeVBm!PDR(rfOo&F*U-_?AyHTmEmhHoKs;li
z;)PC1vUbvO>SiuN>UBz(ppKqi)HDVmM7*)-$>WtaKTo1mZDw)b#IapSfYWz*{iVy7
zqem<e0#|cyuYjJO9*s)TJ%DJxqX;viqoOvdBBN#(*JfTH$1ZXpXJV0c$eECTZ*!jP
zy&d_IMf2KTo*12{Mdy=xBd4w(Iecp$;kN^bpU@3PKby{6&a#~SsC_CUYGIn+$~Jx<
zD~pyNvd;5Y<zqGuADcd1nWt^W^?vzTqtBvtY87K=c52enl8b<(sv@k)H3gr3F<lv*
z^(M0Xt0`d}D-zDjvf+QUi9NvjrQd5hG&wmrR@Tm#k@w5`?UlWSPPooDH!!!F9>>Z%
znNN1-iW7BNks<YW>p#{?k_DKvlFXf4qUT57;MZ1`W`-iPD+h^Kp|4-RX5lDd1q<C)
zbH@)Ns#Ck1ktr$}n?YV%b{~D>Z(0?;m(~;i@S*8g%d?reB&L!?@nm-Fam$NJD_`Ev
zw49TyxVM#mRt|%xF6?$j8XhUiW)zLTdhzn**c*K7^9v56O;5J<G%C}3c@<r)(@n@K
z9XPnFOEA6R$HXtE2k4vNPZ;|weQxY+@9gBQsHlJtX}-Y-;M|{_l{FJf?CS2`^ZNB`
zT1R|*e6ihigdlLZ$_>f4>EZUJyZQ1ERhDmU;r;1xk5GO_hD1$mZI2{`Mc+avf)rj1
zXDP63ecoB><+iev#Ppb-=#X6gI>7^Do(VQt*@wR0uO1SRmzURB>}KcBpuwaa5f>*e
z73AjTHn+94`E=#k80IByLR9m`tB3?i^2)6HV5qEub@5Ck|6s3EW2oX(!TwjTUM2f7
z*Rot`-RRV-saf}z$oJ#veJ?UHRE7tmBO}!`H8n3=tVgQy?rbT8JfJMj3@Tdl%h_Mw
z!Y`AKXC%lwnZCQjj&UAsIwUS0Y-m?{;#FGO;6Sis1EmTorfxWk`@ShMb8WOWRl6=h
zUizw`A+Gza<;=%lwhwz9MF{EVM3Iiyh8>@l2`p@;GXoV_?lS``Th<}rV(cC}Hy{W_
z6k$Pj?*HGm(S5^GaI^YB>`ymNMn*-QHZU-_q_0nD)yDSuFxp<~ro4FZ;{9g^mFXiS
z#(`9Y`YZ{F7*>Yxanbv|>=N(`UH7YR4&2ue3gnYEFSp7p<sM1$^-)vnb#A@j05SKL
zI`XzvtG1&6ox(RP?q5IXyA`ZOB&{5XK5DR&U&i{@_Oezj#bdDMa~L5+wSU*7Mp=^c
z;+IcXvP@6Eefw6;#l?lOpdpMuu7zEl=f(PXx3#xd`V0HSQVz2qjw#ovEFl4CL?jm#
zDQpu`;{AnVYmi>r<qD6HM?9RI9Y$_O2@aKfnA>eH+*iML42CIA1=PpOQZjv4CvYmp
zR#sMTA^1$}Y;7Bs0VhXo*#!wB2~Jwaj~`dl)zw`Y*Cz*c<ydIo@%RL{eh-JadfEJG
zw|DG3<yH6h56;#~5?>qUR~lG6a%rQarKJVm-OMw5+Z)2+&?u1*92|`6EpWhCwmf|S
z#{#*;D%rLx*|*`vv=aq5ISuoy<%b79w709nl^A9{c{T=d{D?^Wre(o0ZAHhSGyX(-
z?8_%l_U_xK9wy*X$B6k6$L4HERn(oWW!>G~u_}1a5J5$m1zHbJZHy$*Wum=e{6iXU
zqWyK2_tNZQdp>FR`+Hl1J4}3)hZ~cGV~jm5&&U3>mzu(c+|+x=@o}oJc+iW=FrMr!
zl!oNR{r2m|EYGjwi^H;HKDeXc6ISiz$)|&ppFKMWS!Oz17w0hY$Zgom8pbHziCimk
z9>dtBy?Pb2yu2)+#fCIF?V8sVK@jigf3NKSn`rSr{n>dK@`HotOktQhmoIahnVIPr
z83}4*SDhKO_1lF)3Y|tO-`sd_+VuEmjMH#kXRek0_rQng#^5T;OX>am@iI0+`o<my
zztv$3PJLy{jT{^hurxP9iO^QzBVyOw+d4V`OBDyYtupK9+7P6xLz2pNv-X3&N#Cbs
zZJ*hQE7><StvlcBMN{cqj6~??+8Byij1e~{=XUIhz4zETlc25qGUXLKCJDs*5a`Je
z=qpd(`){WN?H2#G7-?y^AJrX7y(M)7t%U$2l~7aX<PM`?FsVr3<g_c9QAaZxY90;-
zm38~WE4k^Cxnu~dpy1#XK(O~~XbG;Yto%zpbWg&4Uu>F25}g|4o12<;@XOkUkeA}g
zfshassE{r#M16gIp3k2@(_$btY!kQ<5!p=&$Uleuj;N@p=Js}OFZvYi4Ku^2I5|6e
zHkPPwXvlxz!Ub9k&<!)jVFYn5zO(In_(Nt^KGlr)y|OR#$k_D@Op_S4FdpH#wSxx_
z(rRnGy;+b^n@tNA2;?$ud43GzI@#HoZ)bAUgrwB|+K?wrKQnB0b~aYpLhI<2HzN18
z3GRVacbIscukPZKD@|Y_+|h>`2)DJ^Pbcvxx=n?kUx#{_Zd|HBTB^5GS1VtoQ-$%%
z^MP^F7V2=F^YJprguK7JQ_|jwgyY%y39y8;`2`j>j`D$uU>sC``G*^j#};SCEuh{X
z`sjwRXnj~HU2m#=&M=j+h3|?$iuhTl<D9DW9Rx%p3qS4mV#xq6Am?1rx_XvD4?U58
zfzwYlb05MIhW-Q&$Z_Ktzwyd?H|fy3vb}uk47+xCUF}{c(v-GAm7o>EuX3o8DA7pg
zMfz&c><i>1{C(rV{Zlv>#qg2`P$79n&c34`w>-Uouo`qkosRO~a<ri<$&AH9M2~ex
zRS}oLCHZ8p{-_*+blJ>CRS^rv!WY@>2`pjt|K-{L>56~ZHx@I^LP)d$;+B>A>~$;t
z-j<_b&AXq$ToH5+B7V4T5?!@zC0%8kfI{z>&$CrV<-?xm0jdP!lL#y2ZR)dQa6llK
zguemf9eu&D?kdF&zqSqluE)CLam)@vnM3Fw6XDDy=7%5$o?mBs1H6PZ8<ko9@!c;t
zfEI`P6CLHl_3?_6(k|m|hmIT>w_jSLlao>dzaa|S8xl!nxZx2$z*B%l5?*GTpu#Tw
zYYr$JId@J!kGBnar0AAs;<MbEG;C>#5PA=FCb4)A$5|*`Gs&|y*dnZ=BxBRP3oW&b
zg&*kVr2HM-6bVV~ZEd?aIZyd-qrF*e)({Q_tW~Zg3JU|{nfc@|Jd9I92nmNf-4N^;
zU<}5tzhtqygTy^ZA=*)W`dm6H_$(4QQ`AF)(3(C~Fc{WZ;^DNi+(A|uAB2=a<CF1q
zcz8IyaU0Uq_h3UATGQg~NPLo{MN@y8EY;hI;cH_5jveDTR8_DtTZ@8)L%V?me%N7r
zeBAWYZ+8uncu*W#+yBwfw(wx0=lDzawWSeUT`cC&l%mJ%P(VOuwpp!?UFPf8q0gT`
zU#QKq=_$e+BLr)4xNqsW^Sz!EuhCL-CEwO4!DBd%QmTFG)G0N!5CQqw7x_PMEWcqR
z5^e?A<22IntpM8j%*FF6x!WI?kx^Gu^YHM9jE=5^9L6opex7-4hZhH;%biJH^Q=q4
zVkJuFThJ^`Oa)e#2<7?UxoEZWByO*@r$7(nYzsc!J0u|y_Fv%FjjOP3hdYtFn^KSZ
zjrH_)10|Jz$bmc_IkXs|3Pj+{|4JSC-$t|3myhpa72Rxs(H`a4x@VWPMMH>DvFqm_
zF^Z&TW)jg@fLR>qvPcSr5Ukj(e~t6lv12p_I;Dx;v(@e^Gn5oa0H~sz9mZZZLa9Le
zNZ)8*A}r*#j+kSH+<-&G3-K})&ti=CK=ICU9&5Qm^7gpq+GWPcIrHWHc1jFs{uzj&
zQ`fv_DwcSG)l6(!09Ipl`*u)3z>`z{P*a~@h#zQra!N;6x0TdbIFiT@Dk${pM8}Ak
zIzqU1=${gR7=Ypwx~xG#LE9K%yqJA^`1y&;a-dI|ferbT78c5TodRm3xpVrP6zcb$
zs|PAsaadW)$jB%%CMMqmpsL97B)a^INXh4D!-1Eni9tYHRRD`yNu}A@Qg9~I>Hbo|
z1uDQl2wtJDS+}jMQeM7%xta|ymEN!o3Ev8F8+TsRw(xWGmE7BUO7vbQ2Pgrq1^AYy
zn$2Wi=9OYs8%7?=HH)W7{2xAi(7@L{c<=z!#f1IukTn8T24!v3m0?+o`<IUkNfY9d
zlEjvlmP7C&?=)O~hkdgt2so6h<oA180rDxb@gj$Ap_mhWG%z6GrEaQg=NdZYTV;E9
zaB|YAzP|Hgtzza434Aas7@O|5?qA-mcTxgMV7uEIF_(%PoRYp?Mz~Vq?r_w&L~cD7
z3o^fjogoc~Az<{5#J0g+$<TEg`?3e|RAkw3z@n;b;1G=}3zc@%ncT95$H5MP3+&;F
ztUJy$j0Uw_aq5ZeY;9>D^Et)AO@ST1A8Z*7y6CX}K<z(4hQrPcs?j!7uQEQ0DiLrl
z`_@0pdA|feNQ3`F?XTL%vLT#p@ZR>InUvZ&We!y}xsB@(5qZ`fhzRB{n+8@6&qv3K
z=}B!i%6&8D`B#oN83Bx1p`eT+PD#|(M^pOUWINcAkn&Q?NIx;ABH3?B1o1P*IiI#n
zFSI%PZuM#|mxnQTFh@}Ym<>^YVeEXWcb7{DiPir!*uTKwA0MFrd4IagGccZ(8K#q>
zej1dw(o8v*pe^~+J^tix(yK;1v##|mXlam7oh>apxx}<WDuQ{0osz^eOQc+<db@m<
zh|EQzfEHL6!t9J<LaKI>5Qx-i=OhTX7wPGB&1wjud<uwfAX<)5c)GWRZ)c2=tHH`5
z1sitb{XIb^W_x@4EgKsfa({1cLNI3RVBWMc$}klA(Dmmv(s`H?h<qG<u(b@mY)0uq
zx7=KvG>RNEbYCVN1LWum*#|#Q>!H7(zx*MD+%J0m0{)F!edG9G*^~F*=M^R}Xu~>@
zOgPoF_lP2miAUAcuN%1r(sLSuQ2p{b#>A(Mf$BVn{_ad!Sy>Hm6xuDGjr6p(aIfDc
z7<QUU`l)iOuyY$eZ?a!_?y?M6L^3QJRE-bv6OtD0N_y`_)gp>f7`8W0pB^}S_N=24
zUupFph7a=SNM{bAcpUhAFr*9?1XkR%4uWLA0;1hZ_Vu;y%6_)G++#-;+`4ruY76ix
zINv&Yt>7m#v4F~E!JQVYF7$j8c#vdDT0f{@X36YX??Ha*E%Mq-<5Yd3vLY&VQPfhJ
zt=gaEHr>yA;`XI$4i&p+QiNBlG)fTFcOaw^6kKlKl#0XklR$la7gCFSstBQc!~c!V
zd@UVO?ST9fS$O-_rKmC7uL4t)JRd29+GJuUgojSK;f+)`trMA->Y)n-JWTpP@jK<O
z*q8a%mx|^)&L?N0W)M=uyJhaAgmXs+`_q%iNy(i5j-ydl)ZO+GxjqSUdHPQ+7ZEOP
zP#t{+tGH`WBZeMJ@k=ZtTgYLwSXusV8Ahl+_8)s~HRmFuul<M*rKq0a{RJ{qxNRLs
zRbazSV?N53=7hG5s-r)94t9G3zsA+u=&o>Z5RBE3liS;R+NJo7bdF3Ro}9n)!Y_Au
zj04g#Neb$rdv>2beHt=P+qpnFqPiTrHi)eR@|E)T?eW{UZx7JPOnl;VcBJYN@Qmo_
z=vX-iQ>dMbFY=l5sn}4Hl>wfaT9Y*?8QWf7w)NXW!1JmA!S$o)PBtQ1Nh-?EwUd0>
zdZgmcAHDka&{WA>Lo8@VHMJE6ndub*NS5e65V+XjMvAD_0Ur7ac&oG>dDRsqLx8c2
zJXDn^c{qZDw-(FN`6p8@M@3xKO|k!K1ZP9{y^}FF5gC*}m4kt&U@eR>Rb$v_=z>n$
z%p*}ZZjE%cMdgZi^vaxGchc{y6rz0@vL1xd<5f>5C#Oq#dY>Kk5i;cO=Zfb4Ikni2
zR3qu>=^+6Dz$hwUYI^dg!K;%f`SKgCKl6KvWvd9Q=6%h}74@6adt<@kHe4s&7*bPH
z^BoC+k#-h1-~ho!9HbD!B+3-&(^O$_K`hB@3q5ukz!0|xDjcs26Xc&8t`8PwE>av8
zx-2Mh7Kwb!FWZ?rP||CaQ7B5(1%x+&f(Er{Cv~w>wmZ{zCC4krqA|>@HoAhAPZa`A
zgL*ZT{!*{O2@|HL5>zV;pq_W=79^Yt&{H>6EA+gmQ(aGRP|!f4*Hj>M7R3B@S&$)<
z2Us?!rP+7z>vs?)vmYOyH)h>IkjQ~Iwts#TVg9q_=%au1jUmFq^#DYt=UuJXPfmd=
zo^ASP7{&hQ-XDmPN<LP8)h}Ea&Y#W)VLJJEy(5^0s;28i{=P|I-y(J$;F4lgGoa>@
zXQwaZf85{_-{>~h9bve?!6<;T@Z3Vi(V(8L!En*wBW!(|$Zrl{g9o$uZ=)pJneYN|
zDh&Hjyhc8&aaSncr|iAJt%oApuUyA%7a#Lti&i3q-X|9%Y-0KlrUnk|_<K=8;kK%9
z@O9>%8+yN5?<63qS}taIy@{*l{@?;Rj;pon=Vq9>!aRsb7UA%ET45;1-WlPHpeC)(
zrAuv}9H$CLz8k4$L3YeFVzKS%2CY7&CgdUJS%knE<4r@|W7H^A=0H>%4hR#-%++aJ
z+Y23ZUBT!LD=I1ii+fUs4iOy+c(yQ6UteFiCn6$pSwZ(ZrFk4KI8F80OBOAoJlIxN
z9U(>mj2K98Y6t{LSNYrB&4UllNd$sKl7@cAO5EzicRl108yRtbRzVl$Qp1u-Li-5c
zTMn>e<kpd@%nWa%UENl!-d9^#<Nyj6n2DfF&L=9#H^>c!xKCzLOdAs9le)-xc{1Eb
zuQYw*slY=WrNoV5;JiPo)dL?omkLddB(dr12^N)OW5~)|;JJ$fJY-(fD;%eR<)L%&
z;)B7La(XW_q)!wjx5@srxXXm+i<qV)+ngnXKK;#(g!G(K2FtwtfFJ?8_-=bL@cDi)
zMBT@p#SBP7A0ow$O2SyRr+X7Qknj_lzUx+?n^BQ&kTU>A2eA=9<uvm((^px>x^qVo
zqjzLhdC&KgLA#)ez5%@4n0bvd1UY<mSy0t8fAH9BZ`9qeW3gLQ&#n_Gy{OOaF=hP(
zw^@`ed>_l1Oy_;UN*V&bf9z-%3|2Da<tDU!V`+aSXj-hAWi#!y!S>PO37aO)Bo+;>
zDtm<9nL3DUzNX(^Yv|OoOY%(<d;5JQ7o#P``917?8VGj*mw}|+A03{#flWr&1*>)4
z_WAl7r#xnCEK_r)czTps7--D4Eh+xm>-0~_m9VbmmuDqD)S1NM+z%f+)6Z=XAR?1b
zE*g_F#T)y&Z+~ItNbI#wUUz+a)a)`|Y_@m*CVHQ?oI(4+by|_g<R|!|<mZuXVzm!F
z%L48H=?dSnIf4Be;@Oj&$?d_Rr7mGI5k_|x%Z_XaJUUD@I!ZgGzp>0XT~=cG_I-=T
zX^u-*JEE{@Ij0x3=W{cPcMV9gixjWBF4@k{gU_v!Zqh6Ceyw;PJE6w?3ESLJH1)3+
z<XwcxorKbBc1H6b_BE>N56$LgC^C;`y3Akt2Ce63Sqh*tvA8mV6=vhu5`gA_oiC%u
zn+UB6J<atWYof$)&<Sw>>%0QMaSuv3Q9J{qc<%Q&z%rK_2T>!0->wAHpg3K!_*v4I
zFZGe@)s-Idsv$`O4C;qBwOK8#tb#!&o`!^E<f-77)d0v){CYvR{5?fs4<{$|5_@TL
zG-`){##L;{D&t{FVd0A&VdgF0_6grkJ_g-|>9q6|yF!jmP?r2Slu%`Lvc{xVnKmsJ
zlliIdv)Y<Z_7eWf;B-&i`R9hh`Lc&=4lGZ)fr3rZAclbta~OI#J^55022IDrL>PR`
zl{)+si7|Dh!2X5=WgU`z6Rd8Gd0nhG@jg#AiOpRK^d343v<@QlPKCX{{EUj~pdimm
zCd&^#YTc7}48UzxdE2Mh`AcEv8V6(NmF{|;2Y{dK$(IC95jCwdE82x9t<JZi372Kn
zk+H4Ejn=&z=%YrdHgu#SBO{$VFYD;s&*jrbTQ5K3v$$C)Pc0m;c?waUjvOw>T{KQy
zSPrR9TzB!>HGv4x3qiS79eH+Wb0FQQIPbm2ZwPW%w95_IzYprl4io_;eV1w-nosjo
z{BZ4UOJ*v+=Q$p{(<xX)HByK$R^`Z%BLiaCCCy$pDjr1pXbQ1Xek37cq+06&m4wHQ
zx4mr6G-<(8Q{+UIBL-XM;Dw0082dnZGcO{dSG~3ckm=P8{h$vBvVJT2EWlr(>N?9J
zG{**(Co+lI+1Vgw0#HG_w#v{RC4q7VIwbUiblvo5ad!S13x?`?)R#r?gX`I|XQ-3A
z2=UajK$WepHUPcqy*8IbL|xXQ>WIO7@~TZL1G~03=m6MJ32nE5JDVi}m}}F_hUTYZ
zqcs>0XrB(C1%2<Sd*$JKPbEKpK3qTv&zKnRsHR*_zrqREN6o;S9w)PZMz{$RnPe3o
zdNH*kr%{DpTlwoHU^kml8j;k5=Q(z@4d)YTw0CSA)0KtOnC#5jKLyI(9S-BVvdw0I
z?36<n(NWIjhC!jD8R|7p6-{KIEic?ac_1ppp$|&eq5DD=HU)6QSTrQK({eio9$IAt
zqO=C}Dg?BR3yM6?qsDcPdHpwAqCL|@81xSa2K{nq^EfGad>S9eEOkQKaiP^69bt9z
zQ|IC>rQP7lM+FAdLA$w^L>o&@YS^`Vw^^PYwSNq6;zi4?8*NUJ28X*&P0jqnE4?lu
z*^Q}==T*kCC_=Tncz7(p>5L4;Z0<kJV-y(&P7>oQI6W3<D1)R^9AESOs_z?BIQ}Y4
z!8#tO=x`=;IWtGW;w;>;3N3}GQnN)!Np7xcKk=qj`A|X~W(lH7%4Bt1UmCv1d>E(f
z<!TNtOCfKR?XX~(vF+;eSCG7%CTp`XJ^3Hn+SHO%@TP;6p_CFTczC>M_YS4mv1CSU
zgi|(UsP*B)hsmI52hzlzyncNQby-iGI6><H+n)>GWcjQvr}SIESp^q}0kLoxLr=BR
z8jkhuMo;IzO)v<oCv&U1rwNl7aGbnM!3|b<GWyxuV(c>*NvJt2m<JY=B-$ki4Mp!`
zWod3O&1Cg9c%S=w_KS|(Bvy;9&JcY=gjYv3#8z}v1iyVkJ2U@DxL2smRjST(YiFl>
z+p6|_Y;$3Ryl9duMbNDr9=A~88c#@o9jG09)TrqAo;`aUG5ZXL1j~fR{jn=IV98uO
z<(qlxwHb}fl}2|6Zqw%E(|t-VH)@A5OHYAF7r+t>fPx7A)mL}Fv;A}_y%O3&rm&>n
zI4&{bDf(!K8=29EK0z=r->Th6mrK|;jV$~tzt7dM*Zw`29CN5n20n$IdtW4spVvXt
zkFs)$4igg$O-*%34u2ZzA)b>cJOBadb|i%--C8JuNQ3g0Q8HVDTD9p{bB`?3r+r`H
z#Kzz<L=v7+3B4T!Yp~8bFvOps)e^JteGMj|`UHDTm&!uybHQ^ERl~8GXfhSKSiw*f
z9UxKizg0H1l`-8g)Hs8-YE<8oF_opYy=q$Z6fps~s5q_Xa;9<L*v9o}97Ap67nI<<
z@A1)Zy_E0bdtt{Dp6?Xma1Y?gAV8M5U`djkv{_B$`E%#aNvI#Y#be~GGaPMbYf1Kj
zwMNrHQF-H#<l2Cwf^}L(hU?BVvXEQw<2c{tRNorth?_zBB?WUvv_|0^KbHA5dQYRW
zZ8*DDq87b)9PrF*VGW0O64*bYR?SLR(yD<ugm4^>i8sax<co!Q{lgDF3>bI})HNq5
zKp&-H+^=9h=;;$Rq4R4#J5(LFS)xEB2gFMqJipN+i7p2zQ9U!yv?Ay{;Dk98m^kqN
z#GRnpSU}~Lg5ZQ28nFbe$M~gB+kHDpOo$ccP3Zo`6hfVsVuGxUFPz9ImRe=d8_`=`
z^4!FpnnoD{s#lLDyU*7&6PLw@Xs?Wvc=8=A#nN#L7P>L(8>ED}#)vb>F*gGnIr6Z$
zcx^AV*aeM6Py{`%b~Dj;Z3(Ax51BMJQWkB*GQCCjg?V}~X!&z<2R4;qnR8eR5OMK`
zB_tZIIX9n<gHGmKsthbc-DFpe;v};ZM>HcXR!Evd2T~e+U|^WMHbHh&a`-e+gOMK$
zkq%ovf%1>fK3`bdKV(K|1EJB^5sj!GTE<x*$SlPG$mc({xc;!#{$IUYVj_tRN&!x(
S9A^;t(>$%EmT~g>o&NytZ){-z

literal 0
HcmV?d00001

diff --git a/experiment/paper1_steps/poster_EVAN2019/shape_prior_check/shape_from_some_experiment.txt b/experiment/paper1_steps/poster_EVAN2019/shape_prior_check/shape_from_some_experiment.txt
new file mode 100644
index 00000000..813b0898
--- /dev/null
+++ b/experiment/paper1_steps/poster_EVAN2019/shape_prior_check/shape_from_some_experiment.txt
@@ -0,0 +1,396 @@
+-0.17986546997305747
+-0.21512821907488117
+-0.22839258514334876
+-0.2498225571434785
+-0.10453885381717745
+-0.11536301350989701
+-0.21098040403900203
+-0.3050659290512925
+-0.15925873857758652
+-0.17268190177825643
+-0.2594056844372453
+-0.28039798550349143
+-0.06722834840436546
+-0.1479118601104086
+-0.21433354559531342
+-0.28211290009969725
+-0.0887680562068616
+-0.2208205997056717
+-0.10494743088832467
+-0.11452045521531878
+0.028945433687082932
+0.0067015122458568815
+-0.1472801324090879
+-0.22098011094150238
+-0.0610051579422065
+-0.09022142187635976
+-0.13520802315098068
+-0.12338337567622575
+-0.011899752730240097
+-0.05310318658946517
+-0.06023136975812889
+-0.06563016191691054
+-0.03974170346172862
+-0.03827310510219756
+-0.21367056961008526
+-0.20203181697856493
+-0.21188535197847605
+-0.1639189163167375
+-0.10616033170445793
+-0.10950723984817709
+0.13299031367466324
+0.13458707050406582
+0.07547256999354447
+0.07553864770823433
+0.23444800760973622
+0.23630698649068987
+-0.17986546997305747
+-0.12857827749906253
+-0.22839258514334876
+-0.23992665528817214
+-0.10453885381717745
+-0.09501113151618684
+-0.21098040403900203
+-0.267924138055065
+-0.15925873857758652
+-0.1632068762563936
+-0.2594056844372453
+-0.29430696329291695
+-0.06722834840436546
+-0.1723186485631052
+-0.21433354559531342
+-0.2820272223403669
+-0.0887680562068616
+-0.1998469561489447
+-0.10494743088832467
+-0.10687122711559201
+0.028945433687082932
+-0.0035694057343077634
+-0.1472801324090879
+-0.16226667617602392
+-0.0610051579422065
+-0.1881631947332648
+-0.13520802315098068
+-0.18064106727025547
+-0.011899752730240097
+-0.044448559505585454
+-0.06023136975812889
+-0.0603700543349135
+-0.03974170346172862
+-0.046393073508512156
+-0.21367056961008526
+-0.18470833159640074
+-0.21188535197847605
+-0.2512896889755669
+-0.10616033170445793
+-0.11765252724444987
+0.13299031367466324
+0.11442126594429894
+0.07547256999354447
+0.06645974620355453
+0.23444800760973622
+0.23607921468817256
+-0.17986546997305747
+-0.11767421855242631
+-0.22839258514334876
+-0.24300684469910894
+-0.10453885381717745
+-0.03744201572720143
+-0.21098040403900203
+-0.24560126821798106
+-0.15925873857758652
+-0.16152378998653077
+-0.2594056844372453
+-0.3007158748567143
+-0.06722834840436546
+-0.15841519158251238
+-0.21433354559531342
+-0.2686522942568237
+-0.0887680562068616
+-0.14987108005772576
+-0.10494743088832467
+-0.11016932479851534
+0.028945433687082932
+-0.018940487347995708
+-0.1472801324090879
+-0.15248545102652994
+-0.0610051579422065
+-0.21025410178814968
+-0.13520802315098068
+-0.17856344747026348
+-0.011899752730240097
+-0.04010152067388792
+-0.06023136975812889
+-0.060863919099542585
+-0.03974170346172862
+-0.04705695832001752
+-0.21367056961008526
+-0.18568613549947116
+-0.21188535197847605
+-0.22060391513735975
+-0.10616033170445793
+-0.12036731892022036
+0.13299031367466324
+0.11117383548536468
+0.07547256999354447
+0.07252132278308279
+0.23444800760973622
+0.23727953535405227
+-0.07167408155700153
+-0.07284068855576813
+-0.09537461461592106
+-0.09129636040505239
+-0.03565240312207722
+-0.03632498473754695
+-0.03820593178143043
+-0.07121744555334567
+-0.04012032394334429
+-0.0421795471212244
+-0.38124963216329727
+-0.4072920662869968
+0.13908474912174423
+0.31108640183679825
+-0.06366732312024873
+-0.08142802214894124
+-0.04043456397190185
+-0.06209554218262815
+-0.10306101124433735
+-0.105061605902908
+0.26495290044739506
+0.2651693537670664
+-0.0554182614697758
+-0.12952225331810446
+-0.006287103759081175
+-0.023577383744064054
+0.09567628293191234
+0.19154502325337455
+0.14014740161102704
+0.09305836929837064
+0.06304171036103452
+0.06811512321802246
+0.24709760299667172
+0.24546539155981217
+0.16576276213650265
+0.16682452330147396
+-0.04716342309438759
+-0.050693380192965404
+-0.018339224787634383
+-0.04390430160195761
+0.034030477584533134
+0.03221808998045728
+-0.006686894148799695
+-0.003653548445140834
+-0.0409527995366322
+-0.04330533674810822
+-0.07167408155700153
+-0.07385322485385248
+-0.09537461461592106
+-0.08691089793940054
+-0.03565240312207722
+-0.035961500074552614
+-0.03820593178143043
+-0.04904968521387157
+-0.04012032394334429
+-0.040058422832121424
+-0.38124963216329727
+-0.4255775255857207
+0.13908474912174423
+0.32948472064054846
+-0.06366732312024873
+-0.06581014041501734
+-0.04043456397190185
+-0.0473115474051234
+-0.10306101124433735
+-0.1034503812140259
+0.26495290044739506
+0.26583070811307496
+-0.0554182614697758
+-0.10382805329318165
+-0.006287103759081175
+-0.06161039910393466
+0.09567628293191234
+0.1707671228679576
+0.14014740161102704
+0.10063236364861382
+0.06304171036103452
+0.07241723165747403
+0.24709760299667172
+0.2464524828762165
+0.16576276213650265
+0.16708560345285467
+-0.04716342309438759
+-0.04223688297742566
+-0.018339224787634383
+-0.05621326427186505
+0.034030477584533134
+0.015112983286377522
+-0.006686894148799695
+0.0084135985045764
+-0.0409527995366322
+-0.041330879541588206
+-0.07167408155700153
+-0.07610750819868832
+-0.09537461461592106
+-0.08595044610907085
+-0.03565240312207722
+-0.036908752736647336
+-0.03820593178143043
+-0.06473281826001054
+-0.04012032394334429
+-0.03788220693764348
+-0.38124963216329727
+-0.42807544746430654
+0.13908474912174423
+0.3384350785635435
+-0.06366732312024873
+-0.0725640855386424
+-0.04043456397190185
+-0.05175116106078634
+-0.10306101124433735
+-0.10222160920574865
+0.26495290044739506
+0.2604726253803136
+-0.0554182614697758
+-0.1147454497634801
+-0.006287103759081175
+-0.08823807685369725
+0.09567628293191234
+0.18364291387325196
+0.14014740161102704
+0.05730792075979853
+0.06304171036103452
+0.07211510826670224
+0.24709760299667172
+0.24675719313983843
+0.16576276213650265
+0.1662832108850845
+-0.04716342309438759
+-0.036858800462481356
+-0.018339224787634383
+-0.06485242301330693
+0.034030477584533134
+-0.001833316760152273
+-0.006686894148799695
+0.008107799213101277
+-0.0409527995366322
+-0.039515631316688726
+0.023153980046429462
+0.028368437709736436
+0.02382835483828985
+0.05668510096794785
+0.05768326430840509
+0.07154526043214227
+0.030043102735118146
+0.04231645787397986
+-0.17175820929152086
+-0.14520180164338675
+0.0028562069807351372
+0.036598025744808194
+-0.029649574367349434
+-0.03599847625320862
+-0.044957908239197805
+-0.029960661741132214
+0.4561885125552172
+0.46062013837736504
+-0.012661463769511036
+-0.02273628741179249
+0.024635517950243817
+0.028985695188802826
+0.06861234003281175
+0.06896891203654665
+0.010549053211388493
+0.008453432608898706
+0.20904284330405942
+0.213052274459237
+0.07335130956548458
+-0.013503329370868473
+0.08037111672081598
+0.07975178731239568
+0.022768386837852123
+0.014691330687326493
+0.18829249971095505
+0.1834474831546366
+-0.09443055073415636
+-0.11433649025507156
+-0.017603453383083965
+-0.02450581033427974
+0.023153980046429462
+0.035206263268310826
+0.02382835483828985
+0.07921884186951642
+0.05768326430840509
+0.07108059383978341
+0.030043102735118146
+0.06285698052791111
+-0.17175820929152086
+-0.16311239604418018
+0.0028562069807351372
+0.09585120070226141
+-0.029649574367349434
+0.008434149867197758
+-0.044957908239197805
+-0.005172810269527511
+0.4561885125552172
+0.4611735283254629
+-0.012661463769511036
+-0.014583956392028854
+0.024635517950243817
+0.016769887124387435
+0.06861234003281175
+0.022837477304427034
+0.010549053211388493
+-0.02916284149819585
+0.20904284330405942
+0.17779321734729808
+0.07335130956548458
+-0.00035516640950627006
+0.08037111672081598
+0.02335057573261063
+0.022768386837852123
+0.0025935572691968556
+0.18829249971095505
+0.1843362900913051
+-0.09443055073415636
+-0.12516161155424058
+-0.017603453383083965
+-0.03341544379480277
+0.023153980046429462
+0.03056795906937061
+0.02382835483828985
+0.06871705590566193
+0.05768326430840509
+0.06826999150314064
+0.030043102735118146
+0.06251686121956282
+-0.17175820929152086
+-0.1583926069310994
+0.0028562069807351372
+0.09549477544547744
+-0.029649574367349434
+0.005458378484294713
+-0.044957908239197805
+-0.007512152372008501
+0.4561885125552172
+0.45922425158297886
+-0.012661463769511036
+-0.012099524333183728
+0.024635517950243817
+0.010166201478040926
+0.06861234003281175
+0.016880844841037745
+0.010549053211388493
+-0.023541225896109897
+0.20904284330405942
+0.16934269602225305
+0.07335130956548458
+0.002218562767923835
+0.08037111672081598
+0.013182217830514917
+0.022768386837852123
+-0.008410124971168734
+0.18829249971095505
+0.16467866697088315
+-0.09443055073415636
+-0.12436056994005414
+-0.017603453383083965
+-0.03351139880158266
\ No newline at end of file
diff --git a/experiment/paper1_steps/poster_EVAN2019/some_experiment_EVAN.py b/experiment/paper1_steps/poster_EVAN2019/shape_prior_check/some_experiment_EVAN.py
similarity index 95%
rename from experiment/paper1_steps/poster_EVAN2019/some_experiment_EVAN.py
rename to experiment/paper1_steps/poster_EVAN2019/shape_prior_check/some_experiment_EVAN.py
index 0dae7739..d83c7f69 100644
--- a/experiment/paper1_steps/poster_EVAN2019/some_experiment_EVAN.py
+++ b/experiment/paper1_steps/poster_EVAN2019/shape_prior_check/some_experiment_EVAN.py
@@ -21,11 +21,11 @@ def main_non_stationary_model_comparison():
     stop_loop = False
     for altitude in POSTER_ALTITUDES[:]:
         for trend_test_class in [GevLocationTrendTest, GevScaleTrendTest, GevLocationAndScaleTrendTest,
-                                 ComparisonAgainstMu, ComparisonAgainstSigma][::-1][:]:
+                                 ComparisonAgainstMu, ComparisonAgainstSigma][:3]:
             vizualiser = get_full_altitude_visualizer(Altitude_Hypercube_Year_Visualizer, altitude=altitude,
                                                       exact_starting_year=1958, reduce_strength_array=False,
                                                       trend_test_class=trend_test_class,
-                                                      )
+                                                      verbose=False)
             # vizualiser.save_to_file = False
             vizualiser.visualize_massif_trend_test_one_altitude(poster_plot=True, write_text_on_massif=False)
             if stop_loop:
diff --git a/experiment/paper1_steps/utils.py b/experiment/paper1_steps/utils.py
index eb7ff555..0634dbf5 100644
--- a/experiment/paper1_steps/utils.py
+++ b/experiment/paper1_steps/utils.py
@@ -13,7 +13,8 @@ def get_full_altitude_visualizer(altitude_hypercube_class, exact_starting_year=N
                                  trend_test_class = GevLocationTrendTest,
                                  offset_starting_year=10,
                                  study_class=CrocusSwe3Days,
-                                 orientation=None):
+                                 orientation=None,
+                                 verbose=True):
     altitudes, first_starting_year, last_starting_year, nb_data_reduced_for_speed, only_first_one, save_to_file, _ = get_full_parameters(
         altitude=altitude, offset_starting_year=offset_starting_year)
     if exact_starting_year is not None:
@@ -23,6 +24,7 @@ def get_full_altitude_visualizer(altitude_hypercube_class, exact_starting_year=N
                                           nb_data_reduced_for_speed, only_first_one, save_to_file, study_classes,
                                           trend_test_class, first_starting_year=first_starting_year,
                                           exact_starting_year=exact_starting_year,
-                                          orientations=[orientation])
+                                          orientations=[orientation],
+                                          verbose=verbose)
     visualizer.reduce_strength_array = reduce_strength_array
     return visualizer
diff --git a/experiment/trend_analysis/univariate_test/abstract_gev_trend_test.py b/experiment/trend_analysis/univariate_test/abstract_gev_trend_test.py
index 95a4345c..1db2c6ee 100644
--- a/experiment/trend_analysis/univariate_test/abstract_gev_trend_test.py
+++ b/experiment/trend_analysis/univariate_test/abstract_gev_trend_test.py
@@ -120,6 +120,13 @@ class AbstractGevTrendTest(AbstractUnivariateTest):
         return self.non_stationary_estimator.margin_function_fitted.get_gev_params(coordinate=np.array([1958]),
                                                                                    is_transformed=False)
 
+    @cached_property
+    def stationary_constant_gev_params(self) -> GevParams:
+        # Constant parameters correspond to any gev params
+        return self.stationary_estimator.margin_function_fitted.get_gev_params(coordinate=np.array([1958]),
+                                                                                   is_transformed=False)
+
+
     @property
     def test_trend_slope_strength(self):
         if self.crashed:
-- 
GitLab