ท้าทาย

รับรายการตัวเลขคำนวณค่าเบี่ยงเบนมาตรฐานประชากรของรายการ

ใช้สมการต่อไปนี้เพื่อคำนวณค่าเบี่ยงเบนมาตรฐานของประชากร:

อินพุต

อินพุตจะเป็นรายการจำนวนเต็มในรูปแบบใด ๆ (รายการสตริง ฯลฯ ) ตัวอย่างบางส่วน:

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

ตัวเลขจะเป็นจำนวนเต็มเสมอ

อินพุตจะเป็น STDIN หรืออาร์กิวเมนต์ของฟังก์ชัน

เอาท์พุต

ผลลัพธ์จะต้องเป็นตัวเลขทศนิยม

กฎระเบียบ

คุณอาจใช้ฟังก์ชันในตัวเพื่อค้นหาส่วนเบี่ยงเบนมาตรฐาน

คำตอบของคุณอาจเป็นโปรแกรมเต็มรูปแบบหรือฟังก์ชั่น

ตัวอย่าง

10035, 436844, 42463, 44774 => 175656.78441352615

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

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

การชนะ

โปรแกรมหรือฟังก์ชันที่สั้นที่สุดชนะ

ลีดเดอร์บอร์ด

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>

คลิป 3

.sk

.sเป็นส่วนเบี่ยงเบนมาตรฐานแยกวิเคราะห์การป้อนข้อมูลในรูปแบบk{1,2,3}

Mathematica, 24 22 ไบต์

Nice, Mathematica มี built-in StandardDevi...oh ... ที่คำนวณค่าเบี่ยงเบนมาตรฐานตัวอย่างไม่ใช่ค่าเบี่ยงเบนมาตรฐานประชากร

แต่ถ้าเราใช้Variance... โอ้ ... ข้อตกลงเดียวกัน

แต่ยังมีอีกในตัวที่เกี่ยวข้อง:

CentralMoment[#,2]^.5&

เย้. :)

วิธีนี้ใช้ได้กับ 22 ไบต์:

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

และนี่สำหรับ 27:

N@RootMeanSquare[#-Mean@#]&

อ็อกเทฟ, 14 ไบต์

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

ลองบนideone

kdb + , 3 ไบต์

dev

หนึ่งใน derviates APL ต้องมีนี้เป็นในตัว

ทดสอบการทำงาน

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 ไบต์

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

สิ่งนี้จะกำหนดฟังก์ชั่นรถไฟที่ไม่มีชื่อและ monadic ซึ่งเทียบเท่ากับฟังก์ชั่นต่อไปนี้

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

ลองพวกเขาออนไลน์บนTryAPL

มันทำงานอย่างไร

รหัสประกอบด้วยรถไฟหลายขบวน

M←+/÷≢

สิ่งนี้กำหนด monadic 3-train (fork) Mที่ดำเนินการ+/(ผลรวมขององค์ประกอบทั้งหมด) และ≢(ความยาว) สำหรับอาร์กิวเมนต์ที่ถูกต้องจากนั้นนำไปใช้÷(การหาร) กับผลลัพธ์โดยส่งคืนค่าเฉลี่ยเลขคณิตของอินพุต

M×M

นี่คือส้อมที่ใช้อีกMอาร์กิวเมนต์ขวาซ้ำนี้เป็นครั้งที่สองและนำไปใช้×(ผลิตภัณฑ์) เพื่อผลลัพธ์ที่ได้กลับมาμ 2

×⍨-(M×M)

นี่เป็นอีกทางเลือกหนึ่งที่คำนวณค่าเฉลี่ยของเลขคณิตตามที่อธิบายไว้ก่อนหน้านี้ใช้×⍨(ผลิตภัณฑ์ที่มีตัวเอง) กับอาร์กิวเมนต์ที่ถูกต้องและนำไปใช้กับผล-ต่าง ( ในที่สุด)

สำหรับการป้อนข้อมูล(x ₁ , ... , x _N ) , ฟังก์ชั่นนี้ผลตอบแทน(x ₁ - μ ² , ... , x _N - μ ² )

*∘.5∘M

ฟังก์ชั่นนี้ประกอบไปคือนำไปใช้ในการโต้แย้งสิทธิของมันแล้วM *∘.5หลังใช้การโต้แย้งที่ถูกต้องเพื่อนำแผนที่ไปใช้aกับa*0.5(รากที่สองของa)

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

สุดท้ายเรามี monadic 2-train (atop) ซึ่งใช้ฟังก์ชั่นที่เหมาะสมก่อนจากนั้นทางซ้ายไปยังผลลัพธ์โดยคำนวณค่าเบี่ยงเบนมาตรฐานดังนี้

R, 41 40 39 36 30 28 ไบต์

รหัส

ขอบคุณบีกเกอร์ , Alex A.และMickyTเป็นจำนวนมาก

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

รหัสเก่า

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) 

นี่ควรจะทำให้ค่าเบี่ยงเบนมาตรฐานของประชากร

Pyth, 20 19 17 13 ไบต์

@.O^R2-R.OQQ2

ขอบคุณ @FryAmTheEggman สำหรับการตีกอล์ฟขนาด 4 ไบต์!

ลองออนไลน์

มันทำงานอย่างไร

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

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

ขอบคุณ @aditsu สำหรับการเล่นกอล์ฟ 1 ไบต์!

ลองใช้ออนไลน์ในล่าม CJam

มันทำงานอย่างไร

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 ไบต์

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

วิธีการที่แตกต่างกันเล็ก ๆ น้อย ๆ มากกว่าเดนนิสแก้ปัญหา Dyalog APL สิ่งนี้ควรใช้กับการติดตั้ง APL ใด ๆ

สิ่งนี้จะสร้างฟังก์ชัน monadic ที่ไม่มีชื่อซึ่งคำนวณเวกเตอร์ ( x - µ ) ²ในขณะที่2*⍨⍵-+/⍵÷≢⍵หารสิ่งนี้ด้วยN ( ÷≢⍵) รับผลรวมของเวกเตอร์นี้โดยใช้+/แล้วใช้สแควร์รูท ( .5*⍨)

ลองออนไลน์

Julia, 26 19 ไบต์

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

สิ่งนี้จะสร้างฟังก์ชั่นที่ไม่มีชื่อที่ยอมรับอาร์เรย์และส่งกลับค่าลอย

Ungolfed ฉันเดา:

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 ไบต์

stdDev(augment(Ans,{mean(Ans

ผมยืมขั้นตอนวิธีการที่จะได้รับค่าเบี่ยงเบนมาตรฐานของประชากรจากส่วนเบี่ยงเบนมาตรฐานตัวอย่างจากที่นี่

ทางออกที่สั้นที่สุดที่ฉันสามารถค้นหาได้augment(คือ 9 ไบต์:

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

Haskell, 61 ไบต์

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

ตรงไปตรงมายกเว้นฟังก์ชั่นความยาวที่ฉันกำหนดเองsum(n>>[1])เพื่อหลอกลวงระบบประเภทที่เข้มงวดของ Haskell

Python 3.4+, 30 ไบต์

from statistics import*;pstdev

นำเข้าฟังก์ชั่น builtin pstdevเช่น

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

JavaScript (ES6), 73 ไบต์

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

เยลลี่ไม่ใช่การแข่งขัน

11 ไบต์ คำตอบนี้ไม่ใช่การแข่งขันเนื่องจากใช้ภาษาที่โพสต์ความท้าทาย

S÷L
Dz_²ÇN½

นี่คือการแปลโดยตรงของคำตอบ APL ของฉันไปที่ Jelly ลองออนไลน์!

มันทำงานอย่างไร

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 ไบต์

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

นี่เป็นการแปลโดยตรงของคำตอบ APL ของฉันต่อ J

ลองออนไลน์!

Simplex v.0.5 , 43 ไบต์

เพียงแค่เพราะว่า ฉันต้องการเล่นกอล์ฟนี้อีกหนึ่งไบต์

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 ไบต์

รหัส:

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

คำอธิบาย:

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

ตัวอย่าง:

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

ลองออนไลน์ได้ที่นี่

Perl5, 39 38

 16 สำหรับสคริปต์
+22 สำหรับMสวิตช์
+ 1 สำหรับEสวิตช์
= 39

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

ทดสอบในสตรอเบอร์รี่ 5.20.2

โอ้ แต่แล้วฉันก็ตระหนักว่าคุณตอบว่าคำตอบของเราสามารถทำงานได้แทนโปรแกรม ในกรณีนั้น,

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

มีเพียง 38. การทดสอบในสตรอเบอร์รี่ 5.20.2 เป็น

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

Python ขนาด 57 ไบต์

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

รับอินพุตเป็นรายการ

ขอบคุณ @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)

คำอธิบาย

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

ผลลัพธ์

:\>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 ไบต์

เพียงนำสมการไปใช้ในสมการใน 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 ไบต์ (เป็นฟังก์ชัน)

ฉันเพิ่งสังเกตว่ามันได้รับอนุญาตให้เป็นฟังก์ชั่น อ๊ะนั่นช่วยลดคำตอบลงอย่างมาก ฟังก์ชันนี้กำหนดฟังก์ชันSที่รับอาร์เรย์และส่งกลับค่าเบี่ยงเบนมาตรฐานของประชากร ไปอ่านอีกคำหนึ่งเพื่ออธิบาย แต่ข้ามส่วนการแยกวิเคราะห์ ฉันไม่ต้องการที่จะทำมันอีกครั้ง

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

ในฐานะโปรแกรม 212 ไบต์

น่าเสียดายที่ฉันต้องใช้รายการอินพุตเป็นสตริงและแยกวิเคราะห์ด้วยตัวเอง นี่จะเพิ่มคำตอบมากกว่า 100 ไบต์ดังนั้นหากรูปแบบการป้อนข้อมูลบางอย่างนอกเหนือจากรายการที่คั่นด้วยเครื่องหมายจุลภาคได้รับอนุญาตฉันยินดีที่จะได้ยิน นอกจากนี้โปรดทราบว่าเนื่องจากVALเป็นบั๊กซีการมีช่องว่างก่อนเครื่องหมายจุลภาคหรือต่อท้ายสตริงจะหยุดโปรแกรม หลังจากเครื่องหมายจุลภาคหรือที่จุดเริ่มต้นของสตริงเป็นเรื่องปกติ

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 และอธิบาย:

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)

ความจริง 137 ไบต์

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

ฟังก์ชั่น m () จะคืนค่าเฉลี่ยของรายการในอินพุต ทั้งฟังก์ชั่นที่มีข้อผิดพลาดจะคืนค่า% i ให้กับ sqrt จินตภาพคงที่ (-1) รหัสสำหรับการทดสอบและผลลัพธ์ [แต่ผลลัพธ์ถ้ามันโอเคมันเป็นส่วนจริงของจำนวนเชิงซ้อน]

(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 ไบต์

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

ลองออนไลน์!

ใช้lรายการของจำนวนเต็มและNจำนวนของจำนวนเต็ม

Pyt , 13 ไบต์

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

ใช้สูตรสำหรับการเบี่ยงเบนมาตรฐาน