Tantangan

Diberikan daftar angka, hitung simpangan baku populasi dari daftar tersebut.

Gunakan persamaan berikut untuk menghitung simpangan baku populasi:

Memasukkan

Input akan daftar bilangan bulat dalam format apa pun (daftar, string, dll.). Beberapa contoh:

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

Angka-angka akan selalu berupa bilangan bulat.

Masukan akan ke STDIN atau argumen fungsi.

Keluaran

Output harus berupa angka floating point.

Aturan

Anda dapat menggunakan fungsi bawaan untuk menemukan standar deviasi.

Jawaban Anda dapat berupa program lengkap atau fungsi.

Contohnya

10035, 436844, 42463, 44774 => 175656.78441352615

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

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

Kemenangan

Program atau fungsi terpendek menang.

Papan peringkat

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

.sadalah standar deviasi, kmem-parsing input dalam bentuk {1,2,3}.

Mathematica, 24 22 byte

Bagus, Mathematica memiliki built-in StandardDevi...oh ... yang menghitung standar deviasi sampel, bukan standar deviasi populasi.

Tetapi bagaimana jika kita menggunakan Variance... oh ... kesepakatan yang sama.

Tapi ada built-in lain yang terkait:

CentralMoment[#,2]^.5&

Yay. :)

Ini juga berfungsi untuk 22 byte:

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

Dan ini untuk 27:

N@RootMeanSquare[#-Mean@#]&

Oktaf, 14 byte

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

Cobalah di ideone .

kdb + , 3 byte

dev

Salah satu yang diturunkan APL harus memiliki ini sebagai bawaan .

Uji coba

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←+/÷≢

Ini mendefinisikan kereta fungsi monadik tanpa nama, yang setara dengan fungsi berikut.

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

Cobalah secara online di TryAPL .

Bagaimana itu bekerja

Kode ini terdiri dari beberapa kereta.

M←+/÷≢

Ini mendefinisikan 3-kereta monadik (garpu) Myang mengeksekusi +/(jumlah semua elemen) dan ≢(panjang) untuk argumen yang benar, kemudian menerapkan ÷(pembagian) untuk hasilnya, mengembalikan rata-rata aritmatika dari input.

M×M

Ini adalah garpu lain yang berlaku Muntuk argumen yang tepat, mengulangi ini untuk kedua kalinya, dan berlaku ×(produk) dengan hasil, kembali μ ² .

×⍨-(M×M)

Ini adalah garpu lain yang menghitung kuadrat dari rata-rata aritmatika seperti yang dijelaskan sebelumnya, berlaku ×⍨(produk dengan dirinya sendiri) untuk argumen yang benar, dan akhirnya berlaku -(perbedaan) untuk hasilnya.

Untuk input (x ₁ , ..., x _N ) , fungsi ini mengembalikan (x ₁ - μ ² , ..., x _N - μ ² ) .

*∘.5∘M

Fungsi tersusun ini berlaku Muntuk argumen yang benar, lalu *∘.5. Yang terakhir menggunakan currying argumen kanan untuk menerapkan input peta ake a*0.5(akar kuadrat dari a).

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

Akhirnya, kita memiliki kereta 2 monadik (di atas), yang menerapkan fungsi kanan terlebih dahulu, lalu kiri ke hasilnya, menghitung standar deviasi sebagai berikut.

R, 41 40 39 36 30 28 byte

kode

Berkat gelas kimia , Alex A. dan MickyT untuk banyak byte.

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

kode lama

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) 

Ini harus menghasilkan deviasi standar populasi.

Pyth, 20 19 17 13 byte

@.O^R2-R.OQQ2

Terima kasih kepada @FryAmTheEggman untuk bermain golf 4 byte!

Cobalah online.

Bagaimana itu bekerja

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

Terima kasih kepada @aditsu karena bermain golf 1 byte!

Cobalah online di juru bahasa CJam .

Bagaimana itu bekerja

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*⍨⍵-+/⍵÷≢⍵)÷≢⍵}

Pendekatan yang sedikit berbeda dari solusi APL Dennis 'Dyalog . Ini harus bekerja dengan implementasi APL apa pun.

Ini menciptakan fungsi monadik tanpa nama yang menghitung vektor ( x - μ ) ² sebagai 2*⍨⍵-+/⍵÷≢⍵, membaginya dengan N ( ÷≢⍵), mengambil jumlah vektor ini menggunakan +/, dan kemudian mengambil akar kuadrat ( .5*⍨).

Cobalah online

Julia, 26 19 byte

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

Ini menciptakan fungsi tanpa nama yang menerima array dan mengembalikan float.

Tidak disatukan, saya kira:

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

Saya meminjam algoritma untuk mendapatkan simpangan baku populasi dari simpangan baku sampel dari sini .

Solusi terpendek yang dapat saya temukan tanpa augment(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)

Mudah, kecuali mungkin fungsi panjang kustom saya sum(n>>[1])untuk menipu sistem tipe ketat Haskell.

Python 3.4+, 30 byte

from statistics import*;pstdev

Mengimpor fungsi bawaan pstdev, misalnya

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

Jelly , tidak bersaing

11 byte Jawaban ini tidak bersaing, karena menggunakan bahasa yang mengatasi tantangan.

S÷L
Dz_²ÇN½

Ini adalah terjemahan langsung dari jawaban APL saya ke Jelly. Cobalah online!

Bagaimana itu bekerja

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=:+/%#

Ini adalah terjemahan langsung dari jawaban APL saya ke J.

Cobalah online!

Simplex v.0.5 , 43 byte

Cuz. Saya benar-benar perlu bermain golf satu byte lagi.

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

Kode:

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

Penjelasan:

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

Contoh:

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

Cobalah online di sini

Perl5, 39 38

 16 untuk skrip
+22 untuk Msakelar
+ 1 untuk Esakelar
= 39

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

Diuji dalam Strawberry 5.20.2.

Oh, tapi kemudian saya menyadari bahwa Anda mengatakan jawaban kami bisa berfungsi, bukan program. Dalam hal itu,

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

baru saja 38. Diuji dalam Strawberry 5.20.2 sebagai

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

Mengambil input sebagai daftar

Terima kasih @ 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)

penjelasan

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

hasil

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

Hanya implementasi langsung dari persamaan di 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 (sebagai fungsi)

Saya hanya memperhatikan itu diizinkan menjadi fungsi. Aduh, itu mengurangi jawaban saya secara dramatis. Ini mendefinisikan fungsi Syang mengambil array dan mengembalikan deviasi standar populasi. Bacalah yang lain untuk penjelasan, tetapi lewati bagian parsing. Saya tidak ingin melakukannya lagi.

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

Sebagai sebuah program, 212 byte

Sayangnya, saya harus mengambil daftar input sebagai string dan menguraikannya sendiri. Ini menambahkan lebih dari 100 byte jawaban, jadi jika beberapa format input selain daftar yang dipisahkan koma diizinkan, saya akan senang mendengarnya. Perhatikan juga bahwa karena VALbuggy, memiliki ruang sebelum koma atau mengekstraksi string akan merusak program. Setelah koma atau di awal string baik-baik saja.

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)

Tidak digabungkan dan dijelaskan:

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)

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

Fungsi m () akan mengembalikan rata-rata daftar dalam input. Kedua fungsi pada pengembalian kesalahan% i konstanta imajiner konstan (-1). Kode untuk pengujian dan hasil. [Tapi hasilnya jika itu ok, itu adalah bagian nyata dari satu bilangan kompleks]

(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

Cobalah online!

Diambil l, daftar bilangan bulat, dan N, jumlah bilangan bulat yang ada.

Pyt , 13 byte

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

Menerapkan rumus untuk deviasi standar