ForensicStatistics7.2: การคำนวณค่าทางสถิติ กรณีญาติพี่น้อง (lecture ศ.นพ.ธานินทร์ ภู่พัฒน์ (ตอนที่ 2)
ตอนที่ 2
การอบรมเครือข่ายนิติพันธุศาสตร์แห่งประเทศไทย
“Forensic Genetic Examination Workshop”
วันที่ 15-16 กันยายน 2554 อาคารสัจธรรม ชั้น 3 สถาบันนิติเวชวิทยา โรงพยาบาลตำรวจ
เรื่อง “สถิติความสัมพันธ์เครือญาติโดยดีเอ็นเอ” ศ.นพ. ธานินทร์ ภู่พัฒน์
พิจารณาค่า Probability ในส่วนต่างๆ
1.1 กรณี Full sibling ของ Pr(aa/ab)
|
Pairs of IBD alleles |
Pr(Zi) |
Pr(Zi) sibling โอกาสเกิด IBDกรณีต่างๆ |
Pr(aa/ab, Zi) |
Product
|
|
0 กรณีไม่ share อัลลีล |
Pr(Z0) |
1/4 |
Pa2 (ไม่ได้รับการถ่ายทอดอัลลีลตรงกันมาเลย โอกาสที่จะมีอัลลีลเป็น aa เท่ากับ Pa2) |
(1/4)Pa2 |
|
1a กรณี share อัลลีล a |
(1/2)Pr(Z1a) |
1/4 |
Pa (ได้รับการถ่ายทอด 1 อัลลีล คือ a โอกาสที่จะได้อัลลีล a อีกเท่ากับ Pa) |
(1/4)Pa |
|
1b กรณี share อัลลีล b |
(1/2)Pr(Z1b) |
1/4 |
0 (ได้รับถ่ายทอดอัลลีล b มาแล้วโอกาสที่จะเป็น aa เท่ากับ 0) |
0 |
|
2 กรณี share 2 อัลลีล |
Pr(Z2) |
1/4 |
0 (ได้รับถ่ายทอดอัลลีล ab มาแล้วโอกาสที่จะเป็น aa เท่ากับ 0) |
0 |
|
|
|
|
|
(¼)Pa(1+Pa) |
หมายเหตุ Pr(Z1) เป็น 1/2 แบ่งเป็นโอกาสที่จะเป็น IBD 1a เท่ากับ 1/4 และ IBD 1bเท่ากับ 1/4
Product = Pr(Zi) x Pr(aa/ab, Zi)
1.2 กรณี Full sibling ของ Pr(aa/aa)
|
Pairs of IBD alleles |
Pr(Zi) |
Pr(Zi) sibling โอกาสเกิด IBDกรณีต่างๆ |
Pr(aa/aa, Zi) |
Product
|
|
0 กรณีไม่ share อัลลีล |
Pr(Z0) |
1/4 |
Pa2 (ไม่ได้รับการถ่ายทอดอัลลีลตรงกันมาเลย โอกาสที่จะมีอัลลีลเป็น aa เท่ากับ Pa2) |
(1/4)Pa2 |
|
1a1 กรณี share อัลลีล a1 |
(1/2)Pr(Z1a1) |
1/4 |
Pa (ได้รับการถ่ายทอด 1 อัลลีล คือ a1 โอกาสที่จะได้อัลลีล a อีกเท่ากับ Pa) |
(1/4)Pa |
|
1a2 กรณี share อัลลีล a2 |
(1/2)Pr(Z1a2) |
1/4 |
Pa (ได้รับถ่ายทอดอัลลีล a2 มาแล้วโอกาสที่จะได้อัลลีล a อีกเท่ากับ Pa) |
(1/4)Pa |
|
2 กรณี share 2 อัลลีล |
Pr(Z2) |
1/4 |
1 (ได้รับถ่ายทอดอัลลีล aa มาแล้วโอกาสที่จะเป็น aa เท่ากับ 1) |
1/4 |
|
|
|
|
|
(1/4).(1+Pa)2 |
1.3 กรณี Full sibling ของ Pr(ab/ab)
|
Pairs of IBD alleles |
Pr(Zi) |
Pr(Zi) sibling โอกาสเกิด IBDกรณีต่างๆ |
Pr(ab/ab, Zi) |
Product
|
|
0 กรณีไม่ share อัลลีล |
Pr(Z0) |
1/4 |
2PaPb (ไม่ได้รับการถ่ายทอดอัลลีล โอกาสที่จะมีอัลลีลเป็น ab เท่ากับ 2PaPb) |
(1/4).2PaPb |
|
1a กรณี share อัลลีล a |
(1/2)Pr(Z1a) |
1/4 |
Pb (ได้รับการถ่ายทอด 1 อัลลีล คือ a โอกาสที่จะได้อัลลีล b เท่ากับ Pb) |
(1/4)Pb |
|
1b กรณี share อัลลีล b |
(1/2)Pr(Z1b) |
1/4 |
Pa (ได้รับถ่ายทอดอัลลีล b มาแล้วโอกาสที่จะได้อัลลีล a เท่ากับ Pa) |
(1/4)Pa |
|
2 กรณี share 2 อัลลีล |
Pr(Z2) |
1/4 |
1 (ได้รับถ่ายทอดอัลลีล aa มาแล้วโอกาสที่จะเป็น aa เท่ากับ 1) |
1/4 |
|
|
|
|
(1/4).(1+Pa+Pb+2PaPb) |
|
1.4 กรณี Full sibling ของ Pr(ab/aa)
|
Pairs of IBD alleles |
Pr(Zi) |
Pr(Zi) sibling โอกาสเกิด IBDกรณีต่างๆ |
Pr(ab/aa, Zi) |
Product
|
|
0 กรณีไม่ share อัลลีล |
Pr(Z0) |
1/4 |
2PaPb (ไม่ได้รับการถ่ายทอดอัลลีล โอกาสที่จะมีอัลลีลเป็น ab เท่ากับ 2PaPb) |
(1/4).2PaPb |
|
1a1 กรณี share อัลลีล a1 |
(1/2)Pr(Z1a1) |
1/4 |
Pb (ได้รับการถ่ายทอด 1 อัลลีล คือ a โอกาสที่จะได้อัลลีล b เท่ากับ Pb) |
(1/4)Pb |
|
1a2 กรณี share อัลลีล a2 |
(1/2)Pr(Z1a2) |
1/4 |
Pb (ได้รับการถ่ายทอด 1 อัลลีล คือ a โอกาสที่จะได้อัลลีล b เท่ากับ Pb) |
(1/4)Pb |
|
2 กรณี share 2 อัลลีล |
Pr(Z2) |
1/4 |
0 (ได้รับถ่ายทอดอัลลีล aa มาแล้วโอกาสที่จะเป็น ab เท่ากับ 0) |
0 |
|
|
|
|
|
(1/2).Pb(1+Pa) |
2. กรณี Half sibling ของ Pr(ab/ab)
|
Pairs of IBD alleles |
Pr(Zi) |
Pr(Zi) sibling โอกาสเกิด IBDกรณีต่างๆ |
Pr(ab/ab, Zi) |
Product
|
|
0 กรณีไม่ share อัลลีล |
Pr(Z0) |
1/2 |
2PaPb (ไม่ได้รับการถ่ายทอดอัลลีล โอกาสที่จะมีอัลลีลเป็น ab เท่ากับ 2PaPb) |
(1/2).2PaPb |
|
1a กรณี share อัลลีล a |
(1/2)Pr(Z1a) |
1/4 |
Pb (ได้รับการถ่ายทอด 1 อัลลีล คือ a โอกาสที่จะได้อัลลีล b เท่ากับ Pb) |
(1/4)Pb |
|
1b กรณี share อัลลีล b |
(1/2)Pr(Z1b) |
1/4 |
Pa (ได้รับการถ่ายทอด 1 อัลลีล คือ b โอกาสที่จะได้อัลลีล a เท่ากับ Pa) |
(1/4)Pa |
|
2 กรณี share 2 อัลลีล |
Pr(Z2) |
0 |
1 (ได้รับถ่ายทอดอัลลีล ab มาแล้วโอกาสที่จะเป็น ab เท่ากับ 1) |
0 |
|
|
|
|
|
(1/4).(Pa+Pb+4PaPb) |
3. กรณี First cousin ของ Pr(ab/aa)
|
Pairs of IBD alleles |
Pr(Zi) |
Pr(Zi) sibling โอกาสเกิด IBDกรณีต่างๆ |
Pr(ab/ab, Zi) |
Product
|
|
0 กรณีไม่ share อัลลีล |
Pr(Z0) |
3/4 |
2PaPb (ไม่ได้รับการถ่ายทอดอัลลีล โอกาสที่จะมีอัลลีลเป็น abเท่ากับ 2PaPb) |
(1/4).6PaPb |
|
1a กรณี share อัลลีล a |
(1/2)Pr(Z1a) |
1/8 |
Pb (ได้รับการถ่ายทอด 1 อัลลีล คือ a โอกาสที่จะได้อัลลีล b เท่ากับ Pb) |
(1/8)Pb |
|
1b กรณี share อัลลีล b |
(1/2)Pr(Z1b) |
1/8 |
Pa (ได้รับการถ่ายทอด 1 อัลลีล คือ b โอกาสที่จะได้อัลลีล a เท่ากับ Pa) |
(1/8)Pa |
|
2 กรณี share 2 อัลลีล |
Pr(Z2) |
0 |
1 (ได้รับถ่ายทอดอัลลีล ab มาแล้วโอกาสที่จะเป็น ab เท่ากับ 1) |
0 |
|
|
|
|
|
(1/8).(Pa+Pb+12PaPb) |
ข้อสังเกตค่า Probability
Pr(G1/G2,Z0) จะมีค่าเท่ากับ Prob(G1)
Pr(G1/G2, Z1) จะมีค่าเท่ากัน เมื่อ G2 เป็น homozygous
Pr(G1/G2, Z2) จะมีค่าเท่ากับ 0 หรือ 1 ขึ้นอยู่กับลักษณะ G1 เหมือนกับ G2 หรือไม่
การหาค่า likelihood ratio
ให้เปรียบเทียบโอกาสที่จะเกิดขึ้นเทียบกับ unrelated (ตัวอย่างข้างล่างนี้มีค่าเป็น 1)
|
คนที่ 1 |
คนที่ 2 |
พ่อ-ลูก |
Full sib |
Half sib |
Cousin |
Unrelated |
|
ab |
cd |
0 |
1/4 |
1/2 |
3/4 |
1 |
|
LR เอาโอกาสในแต่ละกรณี หารด้วย unrelated ซึ่งในตัวอย่างนี้มีค่าเท่ากับ 1 |
0 |
1/4 |
1/2 |
3/4 |
1 |
|
กลับไป ตอนที่ 1
ไปต่อ ตอนที่ 3
ความเห็น (0)
ไม่มีความเห็น