Citat:
Kul respons!Mycket intressant artikel. Att ta hänsyn till osäkerheterna i de olika faktorerna är egentligen ganska självklart, men ändå något som sällan görs.
Frågan här är väl om de kraftigt feluppskattat någon av dessa osäkerheter (givet vårt nuvarande vetande enligt litteraturen). Vad vi vet om uppkomsten av flercelligt liv på jorden, är att det tog väldigt lång tid (så att solen hann bli medelålders). Angående uppkomsten av intelligens vet vi ingenting egentligen.
Får jag gissa, så är vi ensamma i vår galax... Men encelligt liv torde vara hyffsat vanligt.
Det har även skrivits artiklar som ansätter hög sannolikhet för att vi är bland de första. Endast 8% av universums förväntade stjärnpopulation finns eller har funnits i nuvarande epok. Liv i universum är i så fall något som kommer att bli vanligare i framtiden.
PS När du körde dina simuleringar, hur valde du standardavvikelserna? Utgick du från artikeln, eller något eget?
Frågan här är väl om de kraftigt feluppskattat någon av dessa osäkerheter (givet vårt nuvarande vetande enligt litteraturen). Vad vi vet om uppkomsten av flercelligt liv på jorden, är att det tog väldigt lång tid (så att solen hann bli medelålders). Angående uppkomsten av intelligens vet vi ingenting egentligen.
Får jag gissa, så är vi ensamma i vår galax... Men encelligt liv torde vara hyffsat vanligt.
Det har även skrivits artiklar som ansätter hög sannolikhet för att vi är bland de första. Endast 8% av universums förväntade stjärnpopulation finns eller har funnits i nuvarande epok. Liv i universum är i så fall något som kommer att bli vanligare i framtiden.
PS När du körde dina simuleringar, hur valde du standardavvikelserna? Utgick du från artikeln, eller något eget?
Den "toy model" som jag har kollat upp ordentligt med egna beräkningar, är alltså den de redogör för på sid 3-4, för att illustrera den statistiska effekt som är en huvudpoäng med arXiv-artikeln. (Denna inledande modell är alltså inte avsedd till att beskriva verkligheten. Detta kommer senare i artikeln.)
I denna förenklade modell beror sannolikheten för Extra Terrestrial Intelligence p=P(ETI) vid en godtycklig stjärna på 9 slumpfaktorer:
p = f1•f2•...•f9
där alla fi är (a) oberoende och (b) plant fördelade mellan 0 och 0.2. Därmed blir väntevärdet
E(p) = E(f1)•E(f2)•...•E(f9) = 1•10^(-9)
så i en galax med N=10^11 stjärnor kan vi förvänta oss E(N•p)=100 ETI. I artikeln påstås även att med denna modell är det 21.45% chans för 0 ETI i galaxen. Om alla f1 däremot är precis 0.1 så blir det förstås samma E(p) medan sannolikheten för 0 ETI blir drastiskt lägre: 3.7•10^(-44).
Mitt mål var
A: att öht förstå vad som menas och hur de har kommit fram till siffrorna
B: kollräkna för att i bästa fall (!) upptäcka något fel
A. Vad menas egentligen? Mitt svar:
Det finns något korrekt värde för p, det är bara det att det är okänt för oss. Givet det rätta värdet för p, och antalet stjärnor N i en galax eller det observerbara universum etc, så ges antalet m ETI av en Poissonfördelning
pm = exp(-N•p)(N•p)^m/m!
Så med t ex p=10^(-9) och N=10^11 blir
p0 = e^(-100) = 3.72•10^(-44).
Stämmer bra med artikeln.
Det sannolikhetsfördelningarna för f1 - f9 modellerar är vår ignorans om deras rätta värden. Så t ex är ett f1 mellan 0.01 och 0.02 lika sannolik som ett värde mellan 0.06 och 0.07 enligt vad vi vet just nu. Så enligt vad vi vet just nu, vad är det sannolika värdet för p0? Svar:
E(p0) = E(exp(-N•f1•f2•...•f9)),
medelvärdesbildat över de plana fördelningarna för f1, f2, ..., f9.
B. Kollräkning av artikelns resultat:
Detta kan beräknas med Monte Carlo-simuleringar. Vid EN simulering slumpar man alltså ut 9 st tal mellan 0 och 0.2 och multiplicerar ihop. I t ex Excel kan detta göras med funktionen 0.2*SLUMP() för varje faktor fi. Med Excel kan vi sedan göra iaf några hundra sådana beräkningar, och sedan räkna ut vad p0 blir i medel. Hyfsat lätt för alla att göra, iaf om man har någon erfarenhet av formler och statistik i Excel.
Men om kan programmera i något programmeringsspråk kan man förstås göra mycket mer. På det sättet gjorde jag en miljon simuleringar, som vid den senaste körningen gav 21.36 % chans för 0 ETI i galaxen, och med 102.9 ETI som väntevärde. Hyfsat nära artikelns värden på 21.45% resp 100.
Slutligen kan man också analysera detta analytiskt. Dvs beräkna sannolikhetsfördelningen för p, givet fördelningarna för de 9 faktorerna. Matematiskt är det faltningar som behöver göras. Resultatet kan sammanfattas med
p = 0.2^9•e^z ; z<0
där sannolikhetsfördelningen för z ges av
h(z) = e^z•z^8/8! ; z<0.
T ex blir då
E(p) = ∫ 0.2^9•e^(2z)•z^8/8! dz
som efter några P.I. ger E(p)=10^(-9), precis som det ska.
E(p0) kan då beräknas med
E(p0) = ∫ e^(N•p)•h(z)dz
Löste detta numeriskt med det Mathematica-liknande programmet Maxima (med dess function quad_qagi) och fick
E(p0) = 21.426913 %
med en påstådd relativ noggrannhet på 10^(-8). Avviker faktiskt något från artikelns 21.45% men det får nog anses respektabelt (där mitt värde alltså är det mer korrekta).
__________________
Senast redigerad av nerdnerd 2018-07-24 kl. 16:16.
Senast redigerad av nerdnerd 2018-07-24 kl. 16:16.