Rev 20 | Rev 4161 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed
Rev | Author | Line No. | Line |
---|---|---|---|
20 | reyssat | 1 | |
2 | !if $wims_read_parm!=slib_header |
||
3 | !goto proc |
||
4 | !endif |
||
5 | slib_title=Quantile of Normal law |
||
6 | slib_parms=5\ |
||
7 | ,[data]\ |
||
8 | 0,mean\ |
||
9 | 1,ecart-type\ |
||
10 | -5,x1\ |
||
11 | 5,x2 |
||
12 | |||
13 | slib_author=Bernadette PERRIN-RIOU |
||
14 | slib_out=quantile at data |
||
15 | slib_comment=data can be a list of reals ; default $x1 and x2 are (mean+/- 5)/ecart-type |
||
16 | slib_example=[0.2,0.4,0.5],0,1 |
||
3265 | bpr | 17 | slib_require=octave |
20 | reyssat | 18 | !exit |
19 | |||
20 | :proc |
||
21 | |||
22 | !distribute item $wims_read_parm into slib_data,slib_mean,slib_var,slib_x1,slib_x2 |
||
23 | slib_out= |
||
24 | !default slib_mean=0 |
||
25 | !default slib_var=1 |
||
26 | !default slib_x1=$[($slib_mean-5)*($slib_var)] |
||
27 | !default slib_x2=$[($slib_mean+5)*($slib_var)] |
||
28 | slib_data=!declosing $slib_data |
||
29 | |||
30 | slib_out=!exec octave normal_inv([$slib_data], $slib_mean,$slib_var^2) |
||
31 | |||
32 | slib_out=!words2items $slib_out |
||
33 | |||
34 | slib_out=!trim $slib_out |