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 : [email protected]
> 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
[email protected]
https://lists.sourceforge.net/lists/listinfo/maxima-lang-fr
------------------------------------------------------------------------------
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
[email protected]
https://lists.sourceforge.net/lists/listinfo/maxima-lang-fr
