Subversion Repositories wimsdev

!if $wims_read_parm!=slib_header
  !goto proc
!endif
slib_title=Quantile
slib_parms=3\
 ,data\
 [0.1,0.25,0.5,0.75,0.9],ordres des quantiles demandés (réels entre 0 et 1)\
 2,option 1 ou 2

slib_out= Quantiles
slib_comment=Option 1 : donne la j-ème\
  valeur de data ordonné par ordre croissant avec j=ceil(np)\
  si n est la taille de data et p l'ordre du quantile demandé.\
  Option 2 : la formule diffère seulement\
  lorsque np est un entier; dans ce cas, la valeur donnée est\
  la moyenne des j-ème et (j+1)-ème valeurs de data ordonné\
  par ordre croissant.
slib_author=Sophie, Lemaire

slib_example=[0,1,1,7,8,9,4],[0.25,0.5,0.75,0.9,0.1],1\
[0,1,1,7,8,9,4],[0.25,0.5,0.75,0.9,0.1],2\
[-1.5,1,1.5,7,8,9,4,5,6,4],[0,0.2,0.5,0.7,1],1\
[-1.5,1,1.5,7,8,9,4,5,6,4],[0,0.2,0.5,0.7,1],2\
[-2],[0.5]
!exit
:proc
!reset slib_out
!distribute item $wims_read_parm into slib_data,slib_order,slib_option
!set slib_data=!declosing $slib_data
!set slib_order=!declosing $slib_order
!default slib_order=0.1,0.25,0.5,0.75,0.9
slib_data=!sort numeric items $slib_data

slib_datacnt=!itemcnt $slib_data
slib_ordercnt=!itemcnt $slib_order
slib_data=!nospace $slib_data
slib_quant=!values ceil(x*$slib_datacnt)+1 for x in $slib_order
slib_data=$(slib_data[1]),$slib_data,$(slib_data[-1])

!default slib_option=2
!if $slib_option=1
  slib_quantile=!item $slib_quant of $slib_data
!endif

!if $slib_option=2
  slib_quant2 = $slib_quant
  !for slib_i = 1 to $slib_ordercnt
    slib_aux = $[$(slib_order[$slib_i])*$slib_datacnt]
    slib_aux = $[ceil($slib_aux)-$slib_aux]
    !if $slib_aux=0
      slib_quant2=!replace item number $slib_i by $[$(slib_quant[$slib_i])+1] in $slib_quant2
    !endif
  !next
  slib_quant=!item $slib_quant of $slib_data
  slib_quant2=!item $slib_quant2 of $slib_data
  slib_quantile=!exec pari ([$slib_quant]+[$slib_quant2])*0.5
!endif

slib_out=$slib_quantile