weighted_mean(a,w):=a.w/lsum(x,x,w);
weighted_var_maple(a,w):=block([v1:lsum(x,x,w),v2:lsum(x^2,x,w),mu],mu:(a.w)/v1,v1/(v1^2-v2)*(w.(a-mu)^2));
weighted_var_matlab(a,w):=block([v1:lsum(x,x,w),mu],mu:(a.w)/v1,(w.(a-mu)^2)/v1);

weighted_var_maple([61,64,67,70,73],[5,18,42,27,8])$
%,numer;
12.02072173667888

weighted_var_matlab([61,64,67,70,73],[5,18,42,27,8])$
%,numer;
8.5275

D'après Wikipedia (http://en.wikipedia.org/wiki/Weighted_mean#Weighted_sample_variance),

la version de matlab est la variance biaisée, alors que la version de maple est la variance non

biaisée.

 

Jean-Claude Arbaut

> Message du 12/11/12 13:56
> De : "Robert SETIF"
> A : maxima-lang-fr@lists.sourceforge.net
> Copie à :
> Objet : [Maxima-lang-fr] variance d'une liste pondérée
>
> Bonjour.Je ne trouve pas de commande pour la variance d'une liste pondérée.Par exemple [61,64,67,70,73] pondérée par[5,18,42,27,8].
> 8 fois 61, 18 fois 64, ...
> XCAS trouve 8.5275 ainsi que MATLAB.
> Par contre Maple 16 trouve 12.0207.
> Il semmble qu'il n'y ait pas de commande pour la variance d'une liste pondérée dans Mathematica 6 , MuPad, MScilab.MERCI
>
>
>
>


------------------------------------------------------------------------------
Everyone hates slow websites. So do we.
Make your web apps faster with AppDynamics
Download AppDynamics Lite for free today:
http://p.sf.net/sfu/appdyn_d2d_nov


_______________________________________________
Maxima-lang-fr mailing list
Maxima-lang-fr@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/maxima-lang-fr