G. scandens에서 부리 깊이는 전혀 유전되나요?
G. scandens에서 부리 깊이의 유전율이 낮아 보입니다. 관측된 유전율이 우연히 얻어진 값일 수도 있고, 실제로는 이 종에서 부리 깊이가 잘 유전되지 않을 수도 있어요. 이제 이 가설을 검정하겠습니다. 이를 위해 쌍 순열 검정을 수행할 거예요.
이 연습은 강의의 일부입니다
Python으로 하는 통계적 사고 (2부)
연습 안내
- 유전율 복제본 배열을 초기화하세요. 10,000쌍의 순열 복제본을 생성합니다.
- 복제본을 생성하는
for루프를 작성하세요.np.random.permutation()을 사용해bd_parent_scandens배열을 섞으세요.- 이전 연습에서 작성한
heritability()함수를 사용해, 섞인 배열과bd_offspring_scandens배열 사이의 유전율을 계산하세요. 결과를 복제본 배열에 저장하세요.
- p-value를, 이전 연습에서 계산한 관측값
heritability_scandens보다 큰 복제본의 개수로 계산하세요.
실습형 인터랙티브 연습
이 예제를 이 샘플 코드를 완성하여 풀어보세요.
# Initialize array of replicates: perm_replicates
perm_replicates = ____
# Draw replicates
for i in range(10000):
# Permute parent beak depths
bd_parent_permuted = ____
perm_replicates[i] = ____
# Compute p-value: p
p = np.sum(____ >= ____) / len(____)
# Print the p-value
print('p-val =', p)