メニュー
Hardy-Weinberg equilibrium expected genotype counts and heterozygosity calculation. Use when computing expected genotype frequencies, expected heterozygosity, or observed heterozygosity from allele counts.
Codex または Claude でインストール この Prompt をコピーして Codex、Claude、または他のアシスタントに貼り付けると、Skill ページを確認してインストールできます。
SOC 職業分類に基づく
Interact with research artifacts running in separate Docker containers via artifact-runner. Execute commands through HTTP API, read files, and verify artifact functionality.
Evaluate research artifacts running in separate Docker containers via artifact-runner. Access artifacts through HTTP API, execute commands, read files, and analyze PDFs.
Interact with artifact containers via HTTP API for paper evaluation tasks. Execute commands, read files, and list directories in remote artifact environments.
Evaluate research artifacts against NDSS badge criteria (Available, Functional, Reproduced) by checking DOI, documentation, exercisability, and reproducibility requirements.
Extract text, tables, and structured information from PDF documents using pdfplumber, PyPDF2, or pdftotext command-line tools.
Calculate protein dihedral angle (phi/psi) from atomic coordinates using the torsion formula. Use when computing backbone dihedral angles from atomic positions or validating protein geometry.
| name | population-genetics |
| description | Hardy-Weinberg equilibrium expected genotype counts and heterozygosity calculation. Use when computing expected genotype frequencies, expected heterozygosity, or observed heterozygosity from allele counts. |
Given observed genotype counts (AA, Aa, aa), calculate expected genotype counts under Hardy-Weinberg equilibrium and compare heterozygosity.
For alleles A and a with frequencies p and q = 1 - p:
If genotype counts are given, compute allele frequencies first:
Then expected counts = N × genotype frequency.
# Read input: AA count, Aa count, aa count
with open("/root/input.txt") as f:
lines = f.read().strip().split("\n")
obs_AA = int(lines[0].strip())
obs_Aa = int(lines[1].strip())
obs_aa = int(lines[2].strip())
N = obs_AA + obs_Aa + obs_aa
# Compute allele frequencies
p = (2 * obs_AA + obs_Aa) / (2 * N)
q = 1.0 - p
# Expected genotype counts under HWE
exp_AA = N * p * p
exp_Aa = N * 2 * p * q
exp_aa = N * q * q
# Expected heterozygosity = 2 * p * q
exp_het = 2 * p * q
# Observed heterozygosity = Aa / N
obs_het = obs_Aa / N
# Output
with open("/root/output.txt", "w") as f:
f.write(f"Expected AA: {exp_AA:.4f}\n")
f.write(f"Expected Aa: {exp_Aa:.4f}\n")
f.write(f"Expected aa: {exp_aa:.4f}\n")
f.write(f"Expected heterozygosity: {exp_het:.4f}\n")
f.write(f"Observed heterozygosity: {obs_het:.4f}\n")
Expected AA: X
Expected Aa: X
Expected aa: X
Expected heterozygosity: X.XXXX
Observed heterozygosity: X.XXXX
Round to 4 decimal places.