Wyzwanie

Biorąc pod uwagę listę liczb, oblicz odchylenie standardowe populacji na liście.

Użyj następującego równania, aby obliczyć odchylenie standardowe populacji:

Wejście

Na wejściu pojawi się lista liczb całkowitych w dowolnym formacie (lista, łańcuch itp.). Kilka przykładów:

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

Liczby zawsze będą liczbami całkowitymi.

Dane wejściowe będą do argumentów STDIN lub funkcji.

Wynik

Dane wyjściowe muszą być liczbą zmiennoprzecinkową.

Zasady

Możesz użyć wbudowanych funkcji, aby znaleźć odchylenie standardowe.

Twoja odpowiedź może być pełnym programem lub funkcją.

Przykłady

10035, 436844, 42463, 44774 => 175656.78441352615

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

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

Zwycięski

Najkrótszy program lub funkcja wygrywa.

Tabela liderów

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>

Klip , 3

.sk

.sjest odchyleniem standardowym, kanalizuje dane wejściowe w formularzu {1,2,3}.

Mathematica, 24 22 bajtów

Fajnie, Mathematica ma wbudowane StandardDevi...... och, które oblicza odchylenie standardowe próbki, a nie odchylenie standardowe populacji.

Ale co, jeśli użyjemy Variance... och ... tej samej umowy.

Ale jest jeszcze jedno powiązane wbudowane:

CentralMoment[#,2]^.5&

Tak :)

Działa to również dla 22 bajtów:

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

A to dla 27:

N@RootMeanSquare[#-Mean@#]&

Oktawa, 14 bajtów

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

Wypróbuj na ideone .

kdb + , 3 bajty

dev

Jeden z derwiatów APL musiał mieć to jako wbudowane.

Testowe uruchomienie

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 bajtów

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

Definiuje nienazwany, monadyczny ciąg funkcji, który jest równoważny z następującą funkcją.

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

Wypróbuj je online na TryAPL .

Jak to działa

Kod składa się z kilku pociągów.

M←+/÷≢

Definiuje monadyczny 3-pociąg (rozwidlenie), Mktóry wykonuje +/(sumę wszystkich elementów) i ≢(długość) dla właściwego argumentu, a następnie stosuje ÷(dzielenie) do wyników, zwracając średnią arytmetyczną danych wejściowych.

M×M

Jest to kolejny rozwidlenie, które stosuje Msię do właściwego argumentu, powtarza to po raz drugi i stosuje ×(produkt) do wyników, zwracając μ ² .

×⍨-(M×M)

Jest to kolejny rozwidlenie, które oblicza kwadrat średniej arytmetycznej, jak wyjaśniono wcześniej, stosuje ×⍨(produkt sam w sobie) do właściwego argumentu i ostatecznie stosuje -(różnicę) do wyników.

Dla danych wejściowych (x ₁ ,…, x _N ) funkcja zwraca (x ₁ - μ ² ,…, x _N - μ ² ) .

*∘.5∘M

Ta skomponowana funkcja ma Mzatem zastosowanie do jej właściwego argumentu *∘.5. Ten drugi używa poprawnego argumentu curry do zastosowania danych wejściowych ana mapie a*0.5(pierwiastek kwadratowy z a).

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

Wreszcie mamy ten monadyczny 2-ciąg (na szczycie), który najpierw stosuje prawą funkcję, a następnie lewą do wyniku, obliczając odchylenie standardowe w następujący sposób.

R, 41 40 39 36 30 28 bajtów

kod

Dzięki zlewce , Alexowi A. i MickyT za wiele bajtów.

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

stare kody

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) 

Powinno to dać odchylenie standardowe populacji.

Pyth, 20 19 17 13 bajtów

@.O^R2-R.OQQ2

Dzięki @FryAmTheEggman za grę w golfa na 4 bajtach!

Wypróbuj online.

Jak to działa

        .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 bajtów

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

Dzięki @aditsu za grę w golfa z 1 bajtu!

Wypróbuj online w interpretatorze CJam .

Jak to działa

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 bajty

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

Trochę inne podejście niż rozwiązanie Dennis Dyalog APL . Powinno to działać z dowolną implementacją APL.

To tworzy nienazwaną funkcję monadyczną, która oblicza wektor ( x - µ ) ^2, ponieważ 2*⍨⍵-+/⍵÷≢⍵dzieli go przez N ( ÷≢⍵), pobiera sumę tego wektora za pomocą +/, a następnie przyjmuje pierwiastek kwadratowy ( .5*⍨).

Wypróbuj online

Julia, 26 19 bajtów

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

Tworzy to nienazwaną funkcję, która akceptuje tablicę i zwraca liczbę zmiennoprzecinkową.

Przypuszczam, że:

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 bajtów

stdDev(augment(Ans,{mean(Ans

Pożyczyłem algorytm, aby uzyskać odchylenie standardowe populacji z próbki odchylenia standardowego od tutaj .

Najkrótsze rozwiązanie bez którego mogę znaleźć augment(to 9 bajtów:

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

Haskell, 61 bajtów

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

Prosto, z wyjątkiem mojej funkcji niestandardowej długości, sum(n>>[1])która może oszukać system ścisłego pisania Haskella.

Python 3.4+, 30 bajtów

from statistics import*;pstdev

Importuje wbudowaną funkcję pstdev, np

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

JavaScript (ES6), 73 bajty

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

Galareta , niekonkurująca

11 bajtów Ta odpowiedź jest niekonkurencyjna, ponieważ używa języka, który stanowi datę późniejszą wyzwania.

S÷L
Dz_²ÇN½

To jest bezpośrednie tłumaczenie mojej odpowiedzi APL na galaretkę. Wypróbuj online!

Jak to działa

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 bajtów

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

To jest bezpośrednie tłumaczenie mojej odpowiedzi APL na J.

Wypróbuj online!

Simplex v.0.5 , 43 bajty

Po prostu. Naprawdę muszę zagrać w ten jeszcze jeden bajt.

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 bajtów

Kod:

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

Wyjaśnienie:

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

Przykład:

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

Wypróbuj online tutaj

Perl5, 39 38

 16 dla skryptu
+22 dla Mprzełącznika
+ 1 dla Eprzełącznika
= 39

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

Testowane w Strawberry 5.20.2.

Och, ale potem zdałem sobie sprawę, że powiedziałeś, że nasze odpowiedzi mogą być funkcjami zamiast programami. W tym wypadku,

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

ma tylko 38. Testowano w Strawberry 5.20.2 as

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

Python, 57 bajtów

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

Pobiera dane wejściowe jako listę

Dzięki @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)

wyjaśnienie

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

wynik

:\>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 bajtów

Prosta implementacja równania w Fortranie:

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 bajtów (jako funkcja)

Właśnie zauważyłem, że może to być funkcja. Ups, to dramatycznie ogranicza moją odpowiedź. Definiuje funkcję, Sktóra pobiera tablicę i zwraca odchylenie standardowe populacji. Przeczytaj wyjaśnienie drugiego, ale pomiń część dotyczącą analizy. Nie chcę tego więcej robić.

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

Jako program 212 bajtów

Niestety muszę wziąć listę wejściową jako ciąg i sam ją przeanalizować. Dodaje to do odpowiedzi ponad 100 bajtów, więc jeśli dozwolony jest jakiś format wejściowy inny niż lista oddzielona przecinkami, chętnie go usłyszę. Zauważ też, że ponieważ VALjest błędny, spacja przed przecinkiem lub przeciąganie łańcucha przerywa program. Po przecinku lub na początku łańcucha jest w porządku.

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)

Nie golfił i wyjaśnił:

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)

Aksjomat, 137 bajtów

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

Funkcja m () zwróciłaby średnią z listy na wejściu. Obie funkcje po błędzie zwracają% i urojoną stałą sqrt (-1). Kod testu i wyników. [ale wynik, jeśli jest w porządku, jest to prawdziwa część jednej liczby zespolonej]

(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 bajtów

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

Wypróbuj online!

Pobiera llistę liczb całkowitych i Nliczbę obecnych liczb całkowitych.

Pyt , 13 bajtów

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

Implementuje formułę odchylenia standardowego