Simulator for genetisk drift i populationer er et værktøj til at simulere og visualisere neutral genetisk drift ved hjælp af den diploide Wright–Fisher-model. Ved at opsætte en simuleret avlspopulation kan man følge og sammenligne reproducerbare allelfrekvensforløb på tværs af flere uafhængige populationer for at se, hvordan genetisk variation går tabt eller fikseres over generationer. Værktøjet hjælper med at forstå, hvordan tilfældig stikprøveudtagning påvirker allelfrekvenser over tid, og viser simuleret heterozygoti, fiksering og tabsrater uden at hævde at forudsige reelle evolutionære resultater.
Wright–Fisher-modellen og matematisk grundlag
Den matematiske kerne i simulatoren bygger på en neutral diploid Wright–Fisher-model, hvor hver generation udgør en tilfældig stikprøve af 2Ne allelkopier trukket fra den foregående generation.
Ved hver generation udtager Wright–Fisher-overgangen det næste antal alleler og dividerer det med de tilgængelige genkopier: Xₜ₊₁ ~ Binomial(2Ne, pₜ); pₜ₊₁ = Xₜ₊₁ ÷ (2Ne).
Før stikprøveudtagningen beskriver disse momenter rækkevidden af mulige næste generationer: E[Xₜ₊₁] = 2Ne·pₜ; Var(Xₜ₊₁) = 2Ne·pₜ(1 − pₜ); Var(pₜ₊₁) = pₜ(1 − pₜ) ÷ (2Ne).
Hver binomialtrækning bruger nøjagtig Bernoulli-stikprøveudtagning; der anvendes ingen normalfordelingsapproksimation. Beregninger bevarer JavaScript-dobbeltpræcision. Frekvenser ligger i intervaller af 1 ÷ (2Ne); visningen afrunder til højst seks decimaler.
Populationsopsætning og parametre
I brugerfladen under "Populationsopsætning" kan du konfigurere simuleringen ved hjælp af følgende parametre:
- Scenarier: Foruddefinerede indstillinger, herunder "Lille population", "Stor population" og "Sjælden allel".
- Effektiv populationsstørrelse, Ne: Angives i avlindivider. Skal være et helt tal fra 1 til 100.000.
- Startallelfrekvens, p₀: Angives som en andel, 0–1. Skal være mellem 0 og 1, inklusive begge grænser.
- Generationer: Angives i generationer. Skal være et helt tal fra 1 til 2.000.
- Uafhængige populationer: Angives i gentagelser. Skal være et helt tal fra 1 til 200.
- Tilfældighedsfrø: Et tekstfelt på maksimalt 100 tegn, der bruges til at reproducere de simulerede forløb. Knappen "Nyt frø" genererer et tilfældigt input.
For at beskytte browserens ydeevne gælder der en øvre grænse for beregningsarbejdet. Dette setup planlægger ‹draws› allel-kopitrækninger, hvilket overskrider browsergrænsen på ‹limit›. Den maksimale grænse er 5.000.000 planlagte allel-kopitrækninger, beregnet som 2Ne × generationer × gentagelser. Hvis denne grænse overskrides, blokerer simulatoren kørslen.
Fortolkning af resultater og diagrammer
Når simuleringen kører, viser værktøjet statusbeskeder som "Udtager stikprøver til de næste generationer…" efterfulgt af "Simulerede ‹replicates› populationer gennem ‹generations› generationer med frøet ”‹seed›”.". Hvis du ændrer parametre, opfordres du til at "Skift populationsscenarie, eller kør dette eksempel med frø igen.".
Resultaterne visualiseres i et diagram med titlen "Allelfrekvensforløb" og ARIA-labelen "Allelfrekvens efter generation for gentagne simulerede populationer". Diagrammets akser og referencelinjer er defineret således:
- X-aksen: "Generationer".
- Y-aksen: "Allelfrekvens p".
- Referencelinje: "Forventet p = p₀".
- Sti-indikator: "Viser
‹shown›ud af‹total›forløb; alle gentagelser er inkluderet i opsummeringen.".
Under diagrammet præsenteres en række opsummerende statistikker for de simulerede populationer:
- Fikseret ved p = 1: Andelen af populationer, hvor allelen har opnået fuld fiksering.
- Tabt ved p = 0: Andelen af populationer, hvor allelen er helt tabt.
- Ikke absorberet ved T: Populationer, der stadig fluktuerer ved simuleringens afslutning.
- Gennemsnitligt endeligt p: Den gennemsnitlige allelfrekvens på tværs af alle replikater ved det sidste tidspunkt.
- Forventet H ved T: Den teoretiske forventede heterozygoti.
- Gennemsnitligt simuleret H ved T: Den faktiske gennemsnitlige heterozygoti observeret i simuleringen.
- Gennemsnitlig absorptionsgeneration: Det gennemsnitlige generationstal for fiksering eller tab (viser "ikke nået", hvis ingen populationer rammer grænserne).
Formel, indsættelse og præcision
For at gøre beregningerne gennemskuelige viser simulatoren de konkrete matematiske substitutioner for den første generation baseret på dine input:
- Første generation bruger 2 ×
‹population›=‹copies›allel-kopitrækninger. - E[X₁] =
‹copies›×‹frequency›=‹mean›A-kopier. - For den første generation er Var(X₁) =
‹copies›×‹frequency›× (1 −‹frequency›) =‹variance›; SD(p₁) =‹sd›. - Første plottede gentagelse udtog X₁ =
‹sampled›A-kopier, så p₁ =‹sampled›÷‹copies›=‹frequency›.
Heterozygoti og tab af genetisk variation
En af de mest centrale indsigter ved genetisk drift er det gradvise tab af heterozygoti (andelen af uens allelpar i populationen) over tid.
Den forventede andel af uens allelpar falder over t generationer som følger: H₀ = 2p₀(1 − p₀); Hₜ = H₀(1 − 1 ÷ (2Ne))ᵗ.
Dette viser, at hastigheden af det genetiske tab er omvendt proportional med den effektive populationsstørrelse (Ne). I mindre populationer falder heterozygotien hurtigt, og alleler fikseres eller tabes hurtigere end i store populationer.
Metodens antagelser og begrænsninger
Denne side implementerer den neutrale diploide Wright–Fisher-model beskrevet af Nature Education og Stanford Human Genetics-lærebogen. Kilder blev kontrolleret den 6. august 2026.
Modellen bygger på en række strenge forenklinger af virkeligheden:
- Konstant effektiv populationsstørrelse.
- Diskrete, ikke-overlappende generationer.
- Tilfældig, uafhængig stikprøveudtagning.
- Selektiv neutralitet (ingen naturlig selektion).
- Ét ukoblet locus med to alleler.
- Ingen mutation eller migration.
Brug ikke dette som en prognose, hvis selektion, populationsændring, ubalance i køn eller reproduktiv succes, overlappende generationer, indavl, kobling, genflow eller populationsstruktur har betydning. En ingeniørmæssig sikkerhedsfaktor gælder ikke for denne stokastiske biologiske model.
Behandling af personoplysninger
Dine parametre og simulerede populationer bliver i denne browser. Al databehandling og generering af tilfældige forløb sker lokalt på din enhed, og ingen data uploades eller sendes eksternt.
Ofte stillede spørgsmål (FAQ)
Skal jeg indtaste den faktiske eller den effektive populationsstørrelse?
Indtast den effektive populationsstørrelse, Ne: størrelsen på en idealiseret avlspopulation, der ville drifte med samme hastighed som den population, du ønsker at studere. Den faktiske populationsstørrelse kan være meget større, når kønsfordeling, reproduktiv succes, populationsudsving eller populationsstruktur er ujævn.
Er den neutrale sandsynlighed for endelig fiksering lig med startfrekvensen?
Ja, i denne neutrale model uden mutation eller migration har en allel, der starter med frekvensen p₀, en endelig fikseringssandsynlighed på p₀ og en tabssandsynlighed på 1 − p₀. De viste rater dækker dog kun det indtastede antal generationer og gentagelser, så de vil variere med frøet.
Hvorfor kan den gennemsnitlige allelfrekvens forblive tæt på p₀, mens heterozygotien falder?
Neutral drift har ingen foretrukken retning på tværs af mange uafhængige populationer, så den gennemsnitlige frekvens forbliver forventet p₀. Individuelle populationer spreder sig dog mod 0 eller 1, og begge grænser har nul heterozygoti. Variation går derfor tabt, selvom gennemsnittet på tværs af populationerne forbliver næsten uændret.
Kan dette forudsige, hvad en reel population vil gøre?
Nej. Den viser resultater under en bevidst idealiseret, neutral Wright–Fisher-model. Reelle forudsigelser kræver forsvarlige estimater af den effektive størrelse samt eventuel selektion, mutation, migration, ændret populationsstørrelse, overlappende generationer, kobling og populationsstruktur, som har betydning for det specifikke tilfælde.