Sfida

Dato un elenco di numeri, calcolare la deviazione standard della popolazione dell'elenco.

Utilizzare la seguente equazione per calcolare la deviazione standard della popolazione:

Ingresso

L'input sarà un elenco di numeri interi in qualsiasi formato (elenco, stringa, ecc.). Qualche esempio:

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

I numeri saranno sempre numeri interi.

L'input sarà su STDIN o argomenti di funzione.

Produzione

L'output deve essere un numero in virgola mobile.

Regole

È possibile utilizzare le funzioni integrate per trovare la deviazione standard.

La tua risposta può essere un programma completo o una funzione.

Esempi

10035, 436844, 42463, 44774 => 175656.78441352615

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

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

vincente

Vince il programma o la funzione più breve.

Classifica

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è la deviazione standard, kanalizza l'input nel modulo {1,2,3}.

Mathematica, 24 22 byte

Bene, Mathematica ha un StandardDevi...oh ... incorporato che calcola la deviazione standard del campione, non la deviazione standard della popolazione.

E se Varianceusassimo ... oh ... stesso affare.

Ma c'è ancora un altro built-in correlato:

CentralMoment[#,2]^.5&

Sìì. :)

Questo funziona anche per 22 byte:

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

E questo per 27:

N@RootMeanSquare[#-Mean@#]&

Ottava, 14 byte

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

Provalo su ideone .

kdb + , 3 byte

dev

Uno dei derviati APL doveva avere questo come incorporato.

Prova

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 byte

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

Questo definisce un treno di funzioni monadico senza nome, che equivale alla seguente funzione.

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

Provali online su TryAPL .

Come funziona

Il codice è composto da diversi treni.

M←+/÷≢

Questo definisce un monadico 3-train (fork) Mche esegue +/(somma di tutti gli elementi) e ≢(lunghezza) per l'argomento giusto, quindi applica ÷(divisione) ai risultati, restituendo la media aritmetica dell'input.

M×M

Questo è un altro fork che si applica Mall'argomento giusto, lo ripete una seconda volta e applica ×(prodotto) ai risultati, restituendo μ ² .

×⍨-(M×M)

Questo è un altro fork che calcola il quadrato della media aritmetica come spiegato prima, applica ×⍨(prodotto con se stesso) all'argomento giusto e infine applica -(differenza) ai risultati.

Per l'ingresso (x ₁ ,…, x _N ) , questa funzione restituisce (x ₁ - μ ² ,…, x _N - μ ² ) .

*∘.5∘M

Questa funzione composta si applica quindi Mal suo argomento giusto *∘.5. Quest'ultimo utilizza il giusto argomento curry per applicare l'input della mappa aa a*0.5(radice quadrata di a).

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

Infine, abbiamo questo monadico 2-treno (in cima), che applica prima la funzione destra, quindi la sinistra al suo risultato, calcolando la deviazione standard come segue.

R, 41 40 39 36 30 28 byte

codice

Grazie a Becher , Alex A. e MickyT per molti byte.

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

vecchi codici

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) 

Ciò dovrebbe produrre la deviazione standard della popolazione.

Pyth, 20 19 17 13 byte

@.O^R2-R.OQQ2

Grazie a @FryAmTheEggman per il golf off 4 byte!

Provalo online.

Come funziona

        .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 byte

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

Grazie a @aditsu per giocare a golf a 1 byte!

Provalo online nell'interprete CJam .

Come funziona

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 byte

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

Un approccio leggermente diverso rispetto alla soluzione Dyalog APL di Dennis . Questo dovrebbe funzionare con qualsiasi implementazione APL.

Questo crea una funzione monadica senza nome che calcola il vettore ( x - µ ) ² come 2*⍨⍵-+/⍵÷≢⍵, lo divide per N ( ÷≢⍵), prende la somma di questo vettore usando +/e quindi prende la radice quadrata ( .5*⍨).

Provalo online

Julia, 26 19 byte

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

Questo crea una funzione senza nome che accetta un array e restituisce un float.

Ungolfed, immagino:

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 byte

stdDev(augment(Ans,{mean(Ans

Ho preso in prestito l'algoritmo per ottenere la deviazione standard della popolazione dalla deviazione standard del campione qui .

La soluzione più breve che potrei trovare senza augment(è di 9 byte:

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

Haskell, 61 byte

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

Semplice, tranne forse la mia funzione di lunghezza personalizzata sum(n>>[1])per ingannare il rigido sistema di tipi di Haskell.

Python 3.4+, 30 byte

from statistics import*;pstdev

Importa la funzione integrata pstdev, ad es

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

JavaScript (ES6), 73 byte

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

Gelatina , non competitiva

11 byte Questa risposta non è competitiva, poiché utilizza un linguaggio che postdatizza la sfida.

S÷L
Dz_²ÇN½

Questa è una traduzione diretta della mia risposta APL a Jelly. Provalo online!

Come funziona

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 byte

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

Questa è una traduzione diretta della mia risposta APL a J.

Provalo online!

Simplex v.0.5 , 43 byte

Solo perche '. Ho davvero bisogno di giocare a golf con questo byte in più.

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 byte

Codice:

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).

Spiegazione:

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

Esempio:

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

Provalo online qui

Perl5, 39 38

 16 per lo script
+22 per lo Mswitch
+ 1 per lo Eswitch
= 39

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

Testato su Fragola 5.20.2.

Oh, ma poi ho capito che hai detto che le nostre risposte possono essere funzioni anziché programmi. In quel caso,

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

ha solo 38. Testato in Strawberry 5.20.2 come

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

Python, 57 byte

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

Accetta input come elenco

Grazie @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)

spiegazione

<#
$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
#>

risultato

:\>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 byte

Solo un'implementazione diretta dell'equazione 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 byte (come funzione)

Ho appena notato che può essere una funzione. Whoops, questo riduce drasticamente la mia risposta. Questo definisce una funzione Sche accetta un array e restituisce la deviazione standard della popolazione. Vai a leggere l'altro per una spiegazione, ma salta la parte di analisi. Non voglio farlo di nuovo.

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

Come programma, 212 byte

Sfortunatamente, devo prendere l'elenco di input come stringa e analizzarlo da solo. Ciò aggiunge oltre 100 byte alla risposta, quindi se è consentito un formato di input diverso da un elenco separato da virgole, sarei felice di ascoltarlo. Si noti inoltre che poiché VALè errato, avere uno spazio prima della virgola o trascinare la stringa interrompe il programma. Dopo la virgola o all'inizio della stringa va bene.

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)

Ungolfed e spiegato:

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)

Assioma, 137 byte

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 funzione m () restituisce la media dell'elenco in input. Entrambe le funzioni sull'errore restituiscono% i la costante immaginaria sqrt (-1). Codice per test e risultati. [ma il risultato se è ok, è la parte reale di un numero complesso]

(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 byte

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

Provalo online!

Contiene lun elenco di numeri interi e Nil numero di numeri interi presenti.

Pyt , 13 byte

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

Implements the formula for standard deviation