Subversion Repositories wimsdev

!goto $wims_read_parm

:def

title=Weighted standard deviation

synonyme=weighted deviation, deviation weighted, weighted data deviation, standard deviation weighted

input=data2d

!exit

:proc

pari_header=dmat=Mat([$fml2]);w=sum(x=1,$cnt,dmat[x,2]);\

	mean=(dmat[,1]~*dmat[,2])/w;\

	dev=sqrt(sum(x=1,$cnt,dmat[x,2]*(dmat[x,1]-mean)^2)/(w));

result=!exec pari print(dev);\

	print(w);

!distribute lines $result into result, totweight

!exit

:output

<center>

Weighted deviation =

!htmlmath $result

</center> <p>

Data: $cnt weighted numbers { $fml3 } (total weight = $totweight).

<p>

(Deviation of the whole population.)

<p>

!exit