Défi

À partir d'une liste de nombres, calculez l'écart type de population de la liste.

Utilisez l'équation suivante pour calculer l'écart type de la population:

Contribution

L'entrée sera une liste d'entiers dans n'importe quel format (liste, chaîne, etc.). Quelques exemples:

56,54,89,87
67,54,86,67

Les nombres seront toujours des entiers.

L'entrée se fera en STDIN ou en arguments de fonction.

Production

La sortie doit être un nombre à virgule flottante.

Règles

Vous pouvez utiliser des fonctions intégrées pour trouver l'écart type.

Votre réponse peut être soit un programme complet soit une fonction.

Exemples

10035, 436844, 42463, 44774 => 175656.78441352615

45,67,32,98,11,3 => 32.530327730015607

1,1,1,1,1,1 => 0.0

Gagnant

Le programme ou la fonction la plus courte l'emporte.

Classement

var QUESTION_ID=60901,OVERRIDE_USER=30525;function answersUrl(e){return"http://api.stackexchange.com/2.2/questions/"+QUESTION_ID+"/answers?page="+e+"&pagesize=100&order=desc&sort=creation&site=codegolf&filter="+ANSWER_FILTER}function commentUrl(e,s){return"http://api.stackexchange.com/2.2/answers/"+s.join(";")+"/comments?page="+e+"&pagesize=100&order=desc&sort=creation&site=codegolf&filter="+COMMENT_FILTER}function getAnswers(){jQuery.ajax({url:answersUrl(answer_page++),method:"get",dataType:"jsonp",crossDomain:!0,success:function(e){answers.push.apply(answers,e.items),answers_hash=[],answer_ids=[],e.items.forEach(function(e){e.comments=[];var s=+e.share_link.match(/\d+/);answer_ids.push(s),answers_hash[s]=e}),e.has_more||(more_answers=!1),comment_page=1,getComments()}})}function getComments(){jQuery.ajax({url:commentUrl(comment_page++,answer_ids),method:"get",dataType:"jsonp",crossDomain:!0,success:function(e){e.items.forEach(function(e){e.owner.user_id===OVERRIDE_USER&&answers_hash[e.post_id].comments.push(e)}),e.has_more?getComments():more_answers?getAnswers():process()}})}function getAuthorName(e){return e.owner.display_name}function process(){var e=[];answers.forEach(function(s){var r=s.body;s.comments.forEach(function(e){OVERRIDE_REG.test(e.body)&&(r="<h1>"+e.body.replace(OVERRIDE_REG,"")+"</h1>")});var a=r.match(SCORE_REG);a&&e.push({user:getAuthorName(s),size:+a[2],language:a[1],link:s.share_link})}),e.sort(function(e,s){var r=e.size,a=s.size;return r-a});var s={},r=1,a=null,n=1;e.forEach(function(e){e.size!=a&&(n=r),a=e.size,++r;var t=jQuery("#answer-template").html();t=t.replace("{{PLACE}}",n+".").replace("{{NAME}}",e.user).replace("{{LANGUAGE}}",e.language).replace("{{SIZE}}",e.size).replace("{{LINK}}",e.link),t=jQuery(t),jQuery("#answers").append(t);var o=e.language;/<a/.test(o)&&(o=jQuery(o).text()),s[o]=s[o]||{lang:e.language,user:e.user,size:e.size,link:e.link}});var t=[];for(var o in s)s.hasOwnProperty(o)&&t.push(s[o]);t.sort(function(e,s){return e.lang>s.lang?1:e.lang<s.lang?-1:0});for(var c=0;c<t.length;++c){var i=jQuery("#language-template").html(),o=t[c];i=i.replace("{{LANGUAGE}}",o.lang).replace("{{NAME}}",o.user).replace("{{SIZE}}",o.size).replace("{{LINK}}",o.link),i=jQuery(i),jQuery("#languages").append(i)}}var ANSWER_FILTER="!t)IWYnsLAZle2tQ3KqrVveCRJfxcRLe",COMMENT_FILTER="!)Q2B_A2kjfAiU78X(md6BoYk",answers=[],answers_hash,answer_ids,answer_page=1,more_answers=!0,comment_page;getAnswers();var SCORE_REG=/<h\d>\s*([^\n,]*[^\s,]),.*?(\d+)(?=[^\n\d<>]*(?:<(?:s>[^\n<>]*<\/s>|[^\n<>]+>)[^\n\d<>]*)*<\/h\d>)/,OVERRIDE_REG=/^Override\s*header:\s*/i;

body{text-align:left!important}#answer-list,#language-list{padding:10px;width:290px;float:left}table thead{font-weight:700}table td{padding:5px}

<script src="https://ajax.googleapis.com/ajax/libs/jquery/2.1.1/jquery.min.js"></script> <link rel="stylesheet" type="text/css" href="//cdn.sstatic.net/codegolf/all.css?v=83c949450c8b"> <div id="answer-list"> <h2>Leaderboard</h2> <table class="answer-list"> <thead> <tr><td></td><td>Author</td><td>Language</td><td>Size</td></tr></thead> <tbody id="answers"> </tbody> </table> </div><div id="language-list"> <h2>Winners by Language</h2> <table class="language-list"> <thead> <tr><td>Language</td><td>User</td><td>Score</td></tr></thead> <tbody id="languages"> </tbody> </table> </div><table style="display: none"> <tbody id="answer-template"> <tr><td>{{PLACE}}</td><td>{{NAME}}</td><td>{{LANGUAGE}}</td><td>{{SIZE}}</td><td><a href="{{LINK}}">Link</a></td></tr></tbody> </table> <table style="display: none"> <tbody id="language-template"> <tr><td>{{LANGUAGE}}</td><td>{{NAME}}</td><td>{{SIZE}}</td><td><a href="{{LINK}}">Link</a></td></tr></tbody> </table>

Clip , 3

.sk

.s est l'écart type, k analyse l'entrée dans le formulaire {1,2,3}.

Mathematica, 24 22 octets

Nice, Mathematica a une fonction intégrée StandardDevi... oh ... qui calcule l'écart type de l'échantillon, pas l'écart type de la population.

Mais si nous utilisons Variance ... oh ... même accord.

Mais il existe encore une autre fonction intégrée:

CentralMoment[#,2]^.5&

Yay. :)

Cela fonctionne également pour 22 octets:

Mean[(#-Mean@#)^2]^.5&

Et cela pour 27:

N@RootMeanSquare[#-Mean@#]&

Octave, 14 octets

g=@(a)std(a,1)

Essayez-le sur ideone .

kdb + , 3 octets

dev

L'un des derviats APL devait avoir ceci comme intégré.

Essai

q)dev 56, 54, 89, 87
16.53028
q)f:dev
q)f 10035, 436844, 42463, 44774
175656.8
q)f 45,67,32,98,11,3
32.53033

Dyalog APL, 24 23 21 20 19 17 octets

*∘.5∘M×⍨-M×M←+/÷≢

Ceci définit un train de fonctions monadique sans nom, qui est équivalent à la fonction suivante.

{.5*⍨M(×⍨⍵)-M⍵×(M←{(+/⍵)÷≢⍵})⍵}

Essayez-les en ligne sur TryAPL .

Comment ça fonctionne

Le code se compose de plusieurs trains.

M←+/÷≢

Ceci définit un 3-train monadique (fourche) Mqui exécute +/(somme de tous les éléments) et ≢(longueur) pour l'argument de droite, puis applique ÷(division) aux résultats, renvoyant la moyenne arithmétique de l'entrée.

M×M

C'est un autre fork qui s'applique Mau bon argument, le répète une deuxième fois et applique ×(produit) aux résultats, renvoyant μ ² .

×⍨-(M×M)

C'est encore une autre fourchette qui calcule le carré de la moyenne arithmétique comme expliqué précédemment, applique ×⍨(produit avec lui-même) au bon argument, et enfin applique -(différence) aux résultats.

Pour l'entrée (x ₁ ,…, x _N ) , cette fonction renvoie (x ₁ - μ ² ,…, x _N - μ ² ) .

*∘.5∘M

Cette fonction composée s'applique Mà son argument de droite, puis*∘.5 . Ce dernier utilise le currying d'argument droit pour appliquer l'entrée de carte aà a*0.5(racine carrée de a).

(*∘.5∘M)(×⍨-(M×M))

Enfin, nous avons ce train monadique à 2 trains (au sommet), qui applique la fonction droite d'abord, puis la gauche à son résultat, calculant l'écart type comme suit.

R, 41 40 39 36 30 28 octets

code

Merci au bécher , Alex A. et MickyT pour beaucoup d'octets.

cat(sd(c(v=scan(),mean(v))))   

anciens codes

v=scan();n=length(v);sd(v)/(n/(n-1))**0.5
m=scan();cat(sqrt(sum(mean((m-mean(m))^2))))
m=scan();cat(mean((m-mean(m))^2)^.5) 

Cela devrait donner l'écart type de la population.

Pyth, 20 19 17 13 octets

@.O^R2-R.OQQ2

Merci à @FryAmTheEggman pour avoir joué au golf sur 4 octets!

Essayez-le en ligne.

Comment ça fonctionne

        .OQ    Compute the arithmetic mean of the input (Q).
      -R   Q   Subtract the arithmetic mean of all elements of Q.
   ^R2         Square each resulting difference.
 .O            Compute the arithmetic mean of the squared differences.
@           2  Apply square root.

CJam, 24 22 21 octets

q~_,_@_:+d@/f-:mh\mq/

Merci à @aditsu pour avoir joué au golf sur 1 octet!

Essayez-le en ligne dans l' interpréteur CJam .

Comment ça fonctionne

q~                    e# Read all input and evaluate it.
  _,                  e# Copy the array and push its length.
    _@                e# Copy the length and rotate the array on top.
      _:+d            e# Copy the array and compute its sum. Cast to Double.
          @/          e# Rotate the length on top and divide the sum by it.
            f-        e# Subtract the result (μ) from the array's elements.
              :mh     e# Reduce by hypotenuse.
                      e# a b mh -> sqrt(a^2 + b^2)
                      e# sqrt(a^2 + b^2) c mh -> sqrt(sqrt(a^2 + b^2)^2 + c^2)
                      e#                           = sqrt(a^2 + b^2 + c^2)
                      e# ⋮
                 \mq/ e# Divide the result by the square root of the length.

APL, 24 octets

{.5*⍨+/(2*⍨⍵-+/⍵÷≢⍵)÷≢⍵}

Une approche légèrement différente de la solution Dyalog APL de Dennis . Cela devrait fonctionner avec n'importe quelle implémentation APL.

Cela crée une fonction monadique sans nom qui calcule le vecteur ( x - µ ) ² as 2*⍨⍵-+/⍵÷≢⍵, le divise par N ( ÷≢⍵), prend la somme de ce vecteur en utilisant +/, puis prend la racine carrée (.5*⍨ ).

Essayez-le en ligne

Julia, 26 ans 19 octets

x->std([x;mean(x)])

Cela crée une fonction sans nom qui accepte un tableau et renvoie un flottant.

Non golfé, je suppose:

function f(x::Array{Int,1})
    # Return the sample standard deviation (denominator N-1) of
    # the input with the mean of the input appended to the end.
    # This corrects the denominator to N without affecting the
    # mean.
    std([x; mean(x)])
end

TI-BASIC, 7 octets

stdDev(augment(Ans,{mean(Ans

J'ai emprunté l'algorithme pour obtenir l'écart type de la population à partir de l'écart type de l'échantillon d' ici .

La solution la plus courte sans laquelle j'ai pu trouver augment(est de 9 octets:

stdDev(Ans√(1-1/dim(Ans

Haskell, 61 octets

d n=1/sum(n>>[1])
f a=sqrt$d a*sum(map((^2).(-)(d a*sum a))a)

Simple, sauf peut-être ma fonction de longueur personnalisée sum(n>>[1])pour tromper le système de type strict de Haskell.

Python 3.4+, 30 octets

from statistics import*;pstdev

Importe la fonction intégrée pstdev, par exemple

>>> pstdev([56,54,89,87])
16.53027525481654

JavaScript (ES6), 73 octets

a=>Math.sqrt(a.reduce((b,c)=>b+(d=c-eval(a.join`+`)/(l=a.length))*d,0)/l)

Gelée , non compétitive

11 octets Cette réponse n'est pas concurrente, car elle utilise un langage postdaté.

S÷L
Dz_²ÇN½

Ceci est une traduction directe de ma réponse APL à Jelly. Essayez-le en ligne!

Comment ça fonctionne

S÷L        Helper link. Argument: z (vector)

S          Compute the sum of z.
  L        Compute the length of z.
 ÷         Divide the former by the latter.
           This computes the mean of z.

Dz_²ÇN½    Main link. Argument: z (vector)

Ç          Apply the previous link, i.e., compute the mean of z.
 ²         Square the mean.
   ²       Square all number in z.
  _        Subtract each squared number from the squared mean.
    Ç      Take the mean of the resulting vector.
     N     Multiply it by -1.
      ½    Take the square root of the result.

J, 18 octets

[:%:@M*:-M*M=:+/%#

Ceci est une traduction directe de ma réponse APL à J.

Essayez-le en ligne!

Simplex v.0.5 , 43 octets

Juste parce que. J'ai vraiment besoin de jouer à cet octet de plus.

t[@u@RvR]lR1RD@wA@T@{j@@SR2ERpR}u@vR@TR1UEo   
t[      ]                                     ~~ Applies inner function to entire strip (left-to-right)
  @                                           ~~ Copies current value to register
   u                                          ~~ Goes up a strip level
    @                                         ~~ Dumps the register on the current byte
     R                                        ~~ Proceeds right (s1)
      v                                       ~~ Goes back down
       R                                      ~~ Proceeds right (s0)
                                              ~~ Go right until an empty byte is found
         lR1RD                                ~~ Push length, 1, and divide.
              @                               ~~ Store result in register (1/N)
               wA                             ~~ Applies A (add) to each byte, (right-to-left)
                 @T@                          ~~ Puts 1/N down, multiplies it, and copies it to the register
                    {          }              ~~ Repeats until a zero-byte is met
                     j@@                      ~~ inserts a new byte and places register on it
                        SR                    ~~ Subtract it from the current byte and moves right
                          2E                  ~~ Squares result
                            RpR               ~~ Moves to the recently-created cell, deletes it, and continues
                                u@v           ~~ takes 1/N again into register
                                   R@T        ~~ multiplies it by the new sum
                                      R1UE    ~~ takes the square root of previous
                                          o   ~~ output as number

Prolog (SWI), 119 octets

Code:

q(U,X,A):-A is(X-U)^2.
p(L):-sumlist(L,S),length(L,I),U is S/I,maplist(q(U),L,A),sumlist(A,B),C is sqrt(B/I),write(C).

Explication:

q(U,X,A):-A is(X-U)^2.   % calc squared difference of X and U
p(L):-sumlist(L,S),      % sum input list
      length(L,I),       % length of input list
      U is S/I,          % set U to the mean value of input list
      maplist(q(U),L,A), % set A to the list of squared differences of input and mean
      sumlist(A,B),      % sum squared differences list
      C is sqrt(B/I),    % divide sum of squares by length of list
      write(C).          % print answer

Exemple:

p([10035, 436844, 42463, 44774]).
175656.78441352615

Essayez-le en ligne ici

Perl5, 39 38

 16 pour le script
+22 pour le Mcommutateur
+ 1 pour le Ecommutateur
= 39

perl -MStatistics::Lite=:all -E"say stddevp@ARGV" .1 .2 300

Testé dans Strawberry 5.20.2.

Oh, mais j'ai réalisé que vous aviez dit que nos réponses pouvaient être des fonctions plutôt que des programmes. Dans ce cas,

{use Statistics::Lite":all";stddevp@_}

a juste 38. Testé dans Strawberry 5.20.2 comme

print sub{use Statistics::Lite":all";stddevp@_}->( .1, .2, 300)

Python, 57 octets

lambda l:(sum((x-sum(l)/len(l))**2for x in l)/len(l))**.5

Takes input as a list

Thanks @xnor

PowerShell, 122

:\>type stddev.ps1
$y=0;$z=$args -split",";$a=($z|?{$_});$c=$a.Count;$a|%{$y+=$_};$b=$y/$c;$a|%{$x+
=(($_-$b)*($_-$b))/$c};[math]::pow($x,0.5)

explanation

<#
$y=0                            init
$z=$args -split","              split delim ,
$a=($z|? {$_})                  remove empty items
$c=$a.Count                     count items
$a|%{$y+=$_}                    sum
$b=$y/$c                        average
$a|%{$x+=(($_-$b)*($_-$b))/$c}  sum of squares/count
[math]::pow($x,0.5)             result
#>

result

:\>powershell -nologo -f stddev.ps1 45,67,32,98,11,3
32.5303277300156

:\>powershell -nologo -f stddev.ps1 45,  67,32,98,11,3
32.5303277300156

:\>powershell -nologo -f stddev.ps1 45,  67,32, 98 ,11,3
32.5303277300156

:\>powershell -nologo -f stddev.ps1 10035, 436844, 42463, 44774
175656.784413526

:\>powershell -nologo -f stddev.ps1 1,1,1,1,1,1
0

Fortran, 138 bytes

Just a straightforward implementation of the equation in Fortran:

double precision function std(x)
integer,dimension(:),intent(in) :: x
std = norm2(dble(x-sum(x)/size(x)))/sqrt(dble(size(x)))
end function

SmileBASIC, 105 bytes (as a function)

Je viens de remarquer que c'est autorisé d'être une fonction. Oups, cela réduit considérablement ma réponse. Cela définit une fonction Squi prend un tableau et renvoie l'écart type de la population. Allez lire l'autre pour une explication, mais sautez la partie d'analyse. Je ne veux pas recommencer.

DEF S(L)N=LEN(L)FOR I=0TO N-1U=U+L[I]NEXT
U=1/N*U FOR I=0TO N-1T=T+POW(L[I]-U,2)NEXT RETURN SQR(1/N*T)END

En tant que programme, 212 octets

Malheureusement, je dois prendre la liste d'entrée comme une chaîne et l'analyser moi-même. Cela ajoute plus de 100 octets à la réponse, donc si un format d'entrée autre qu'une liste séparée par des virgules est autorisé, je serais heureux de l'entendre. Notez également que parce que VALc'est bogué, le fait d'avoir un espace avant la virgule ou de suivre la chaîne interrompt le programme. Après la virgule ou au début de la chaîne, c'est bien.

DIM L[0]LINPUT L$@L I=INSTR(O,L$,",")IF I>-1THEN PUSH L,VAL(MID$(L$,O,I-O))O=I+1GOTO@L ELSE PUSH L,VAL(MID$(L$,O,LEN(L$)-O))
N=LEN(L)FOR I=0TO N-1U=U+L[I]NEXT
U=1/N*U FOR I=0TO N-1T=T+POW(L[I]-U,2)NEXT?SQR(1/N*T)

Non golfé et expliqué:

DIM L[0]  'define our array
LINPUT L$ 'grab string from input

'parse list
'could've used something cleaner, like a REPEAT, but this was shorter
@L
I=INSTR(O,L$,",")                 'find next comma
IF I>-1 THEN                      'we have a comma
 PUSH L,VAL(MID$(L$,O,I-O))       'get substring of number, parse & store
 O=I+1                            'set next search location
 GOTO @L                          'go again
ELSE                              'we don't have a comma
 PUSH L,VAL(MID$(L$,O,LEN(L$)-O)) 'eat rest of string, parse & store
ENDIF                             'end

N=LEN(L) 'how many numbers we have

'find U
'sum all of the numbers, mult by 1/N
FOR I=0 TO N-1
 U=U+L[I]
NEXT
U=1/N*U

'calculate our popstdev
'sum(pow(x-u,2))
FOR I=0 TO N-1
 T=T+POW(L[I]-U,2)
NEXT
PRINT SQR(1/N*T) 'sqrt(1/n*sum)

Axiome, 137 octets

m(a:List Float):Complex Float==(#a=0=>%i;reduce(+,a)/#a)
s(a:List Float):Complex Float==(#a=0=>%i;n:=m(a);sqrt(m([(x-n)^2 for x in a])))

La fonction m () retournerait la moyenne de la liste en entrée. Les deux fonctions en cas d'erreur renvoient% i la constante imaginaire sqrt (-1). Code pour le test et les résultats. [mais le résultat si c'est ok, c'est la partie réelle d'un nombre complexe]

(6) -> s([45,67,32,98,11,3])
   (6)  32.5303277300 15604966

(7) -> s([10035,436844,42463,44774])
   (7)  175656.7844135261 4035

(8) -> s([1,1,1,1,1,1])
   (8)  0.0

Python 3 , 49 octets

lambda l,N:(sum((i-sum(l)/N)**2for i in l)/N)**.5

Essayez-le en ligne!

Takes in l, a list of integers, and N, the number of integers present.

Pyt, 13 bytes

←Đ↔Ł↔е-²Ʃ⇹/√

Implements the formula for standard deviation