lin3

Information about lin3

Published on January 15, 2008

Author: Maurizio

Source: authorstream.com

Content

Demonstration of Linkage Analysis:  Demonstration of Linkage Analysis 林明薇 Ming-Wei Lin, PhD 陽明大學醫學系家庭醫學科 台北榮民總醫院教學研究部 Linkage Data Files:  Linkage Data Files Pedigree file Parameter file Slide3:  Pedigree file format (demo.pre) 1 1 0 0 1 1 5 5 5 5 2 2 2 3 4 7 3 4 5 7 2 7 4 6 3 3 1 5 1 10 13 2 2 1 4 5 3 4 2 4 3 4 4 8 1 1 6 6 2 6 4 7 3 3 1 1 1 12 0 0 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 13 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 1 12 2 2 4 5 4 5 2 2 2 3 4 7 1 3 5 6 2 6 6 7 3 3 1 1 1 3 1 12 1 1 3 5 4 5 2 2 3 3 4 4 1 4 6 7 6 7 4 7 3 3 1 5 1 4 1 12 2 2 4 5 4 5 2 2 2 3 4 4 1 3 5 6 2 6 6 7 3 3 1 1 1 9 13 2 1 1 4 5 3 5 2 4 2 4 7 8 1 3 5 6 1 2 4 6 3 3 1 1 2 1 0 0 2 2 4 5 4 5 2 4 3 3 4 7 1 3 1 6 1 6 5 5 3 3 3 5 2 10 22 6 2 1 4 5 5 6 2 2 2 3 4 7 1 2 6 6 1 6 7 9 3 4 3 5 A B C D E F G H I J K L M N O P Q A: pedigree name B: individual Id C: father's ID D: mother's ID E: sex (1= male, 2= female) F: affection status (1 = unaffected, 2 = affected; 0 = unknown) G-Q: marker genotypes Slide4:  Parameter file format (demo.par) 12 0 0 5 << no loci, risk locus, sexlinked(if 1) 0 0.0 0.0 0 << mut locus, mut rate, haplotype freq(if 1) 1 2 3 4 5 6 7 8 9 10 11 12 << order of loci 1 2 << affection, #alleles [amyloidosis] 0.9999 0.0001 << gene freqs 1 << number of liability classes 0.00002 0.9 0.9 << penetrence 3 7 << numbered alleles, #alleles [d1s498] 0.02 0.29 0.07 0.36 0.21 0.04 0.01 3 7 << numbered alleles, #alleles [d1s2635] 0.02 0.05 0.1 0.38 0.38 0.05 0.02 3 4 << numbered alleles, #alleles [d1s2771] 0.02 0.79 0.07 0.12 ABO Genotype and Phenotype:  ABO Genotype and Phenotype Penetrances = Prob(phenotype given genotype) Incomplete Penetrance - Dominant Disease -:  Incomplete Penetrance - Dominant Disease - Incomplete Penetrance - Recessive Disease -:  Incomplete Penetrance - Recessive Disease - Slide8:  Parameter file format (demo.par) cont.1 3 4 << numbered alleles, #alleles [d1s484] 0.33 0.12 0.5 0.05 3 9 << numbered alleles, #alleles [d1s2675] 0.02 0.05 0.31 0.15 0.02 0.4 0.03 0.01 0.01 3 4 << numbered alleles, #alleles [d1s2768] 0.5 0.12 0.33 0.05 3 9 << numbered alleles, #alleles [d1s2844] 0.02 0.12 0.17 0.05 0.07 0.21 0.1 0.14 0.12 3 7 << numbered alleles, #alleles [d1s2878] 0.29 0.21 0.17 0.02 0.29 0.01 0.01 3 11 << numbered alleles, #alleles [d1s2681] 0.02 0.02 0.05 0.14 0.19 0.4 0.02 0.03 0.05 0.05 0.03 3 4 << numbered alleles, #alleles [d1s194] 0.14 0.02 0.72 0.12 Slide9:  Parameter file format (demo.par) cont.2 3 5 << numbered alleles, #alleles [d1s196] 0.64 0.02 0.05 0.05 0.24 0 0 << sex difference, interference (if 1 or 2) 0 0.09 0.028 0.011 0.01 0.017 0.025 0.025 0.011 0.005 0.03 << recombination values 1 0.05 0.45 << rec varied, increment, finishing value Chromosome 1 Map:  Chromosome 1 Map Markers cM D1S498 D1S2635 D1S2771 D1S484 D1S2675 D1s2768 9.0 2.8 1.7 2.5 1.1 qter 3.0 1.1 1.0 cen 2.5 0.4 D1S2844 D1S2878 D1S2681 D1S194 D1S196 Preplink:  Preplink Use the Preplink to make the parameter file Slide12:  1. 鍵入 Preplink。 Slide13:  2. 開啟 Preplink ,得到下列的畫面。 Slide14:  3. 鍵入Enter ,及選擇選項 (k) see or modify loci description。 Slide15:  4. 選取(e) change locus type ,鍵入locus to change 及選取 new locus type後,按 Enter。 Slide16:  5. 選取(a) see or modify a locus ,並鍵入locus number to see or modify。 Slide17:  6. 選取(d) gene frequencies 並鍵入新 gene frequencies 值。 Slide18:  7. 選取(c) penetrances , 並鍵入三組新penetrance值。 Slide19:  8. 選取(e) EXIT。 Slide20:  9. 選取(a) See or modify a locus及鍵入 Locus number。 Slide21:  10. 選取(a) number of alleles 及鍵入正確之 number of alleles。 Slide22:  11. 選取(b) gene frequencies 及鍵入新的 gene frequencies。 Slide23:  12. 選取(c) EXIT。 Slide24:  13. 選取(f) RETURN TO MAIN MENU。 Slide25:  14. 選取(n) Write datafile。 Slide26:  15. 鍵入 Output file name。 Slide27:  16. 選取(o) Exit。 Makeped:  Makeped The pedigree file (*.pre) has to run the Makeped program first before applying to the LINKAGE program. Slide29:  1. 開啟 Makeped,得到下列的畫面。 Slide30:  2. 鍵入Pedigree file、Output file名稱,及選擇 家族是否含有Loop,Proband 是否自動設定。 Makeped:  Makeped makeped <input pedigree filename> <output pedigree filename> n (no loop) Ex: makeped demo.pre demo.ped n Slide32:  LINKAGE I. Analysis Program of LINKAGE: 1. ILINK: estimation of recombination rates and calculation of the maximum lod score. 2. MLINK: lod score tables and risk analysis. 3. LINKMAP: location scores. 4. LODSCORE: estimation of recombination rate and the maximum lod score from two-locus data. It is almost identical to ILINK. Slide33:  II. Steps for running the LINKAGE program: 1. Input pedigree and genotypic data (pedigree file). 2. Description of loci (parameter file). Program code: 1 = CILINK 2 = CMAP 3 = ILINK 4 = LINKMAP 5 = MLINK 6 = LODSCORE 7 = CLODSCORE Slide34:  Locus description: 0 = Quantitative trait 1 = Affected status 2 = Binary factors 3 = Numbered alleles Recombination information: Sex difference: 0 = no sex-difference 1 = constant sex-difference 2 = variable sex-difference Interference: 0 = no interference 1 = interference without a mapping function 2 = user-specific mapping function Slide35:  Analysis: use LINKAGE CONTROL PROGRAM (LCP) to generate command files to control the analysis programs. Linkage Control Program (LCP):  Linkage Control Program (LCP) Slide37:  Window Control Characters Key Ctrl Chr Action F10 ^U Delete the line. ^ ^I Move <up> to the previous entry. V ^K Move <down> to the next entry. Return ^M Move <down> to the next entry. (on entry screens) Prev Scr ^P Move <up> to the previous screen. Next Scr ^N Move <down> to the next screen. F8 ^Z Exit the current session. F9 ^A Abort the current session.. Linkage Control Program (LCP) - MLINK -:  Linkage Control Program (LCP) - MLINK - Slide39:  1. 開啟LCP,出現下列畫面: Slide40:  2. 輸入COMMAND file name, LOG file name, PEDIGREE file name, PARAMETER file name的名稱。 Slide41:  3. 按Ctrl-N,出現下列畫面,選取General pedigrees。 Slide42:  4. 按Ctrl-N,出現下列畫面,選取MLINK。 Slide43:  5. 按Ctrl-N,出現下列畫面,選取Lod score table。 Slide44:  6. 按Ctrl-N,出現下列畫面,選取 No sex difference。 Slide45:  7. 按Ctrl-N,出現下列畫面,輸入Locus order, Other recomb,再按Ctrl-N及Ctrl-Z,即完成LCP之設定。 Slide46:  8. 在Unix 命令列鍵入 demo,即可執行MLINK的分析。 Slide47:  9. MLINK的分析畫面。 Linkage Control Program (LCP) - LINKMAP -:  Linkage Control Program (LCP) - LINKMAP - Slide49:  1. 開啟LCP,得到下列畫面: Slide50:  2. 輸入COMMAND file name, LOG file name, PEDIGREE file name, PARAMETER file name的名稱。 Slide51:  3. 按Ctrl-N,出現下列畫面,選取General pedigrees。 Slide52:  4. 按Ctrl-N,出現下列畫面,選取LINKMAP。 Slide53:  5. 按Ctrl-N,出現下列畫面,選取All intervals。 Slide54:  6. 按Ctrl-N,出現下列畫面,選取No sex difference。 Slide55:  7. 按Ctrl-N,出現下列畫面,輸入Test loci、Order of fixed loci、Recombination fractions,再按 Ctrl-N Ctrl-Z,即完成LCP之設定。 Slide56:  8. 在Unix 命令列鍵入 demo1,即可執行LINKMAP的 分析。 Slide57:  9. LINKMAP的分析畫面。 Slide58:  10. 在DOS下,執行 Table.exe 的功能。 Microsoft(R) Windows 98 (C)Copyright Microsoft Corp 1981-1998. C:\WINDOWS>cd\ C:\>cd ming-wei C:\Ming-Wei>cd vghworkshop C:\Ming-Wei\vghworkshop>cd 2002 C:\Ming-Wei\vghworkshop\2002>cd demo C:\Ming-Wei\vghworkshop\2002\demo>table demo.res /a C:\Ming-Wei\vghworkshop\2002\demo>table demo1.res /a C:\Ming-Wei\vghworkshop\2002\demo> Slide59:  11. 用 WORD 開啟 demo.tab,得到MLINK的結果。 Slide60:  12. 用 WORD 開啟 demo1.tab,得到LINKMAP的結果。 Slide61:  13. 解讀LINKMAP的結果。 GENEHUNTER:  GENEHUNTER Reference:  Reference Kruglyak L, Daly MJ, Reeve-Daly MR, Lander ES (1996) Parametric and non-parametric analysis: a unified multipoint approach. American Journal of Human Genetics 58: 1347-1363. Slide64:  Running GENEHUNTER 1. type gh.sun 2. photo demo.out 3. skip large off 4. load demo.par 5. scan demo.pre 6. total stat het 7. quit Slide65:  1. 鍵入gh.sun ,在npl下鍵入 photo demo.out 、skip large off 、load demo.par 、 scan demo.pre等指令。 Slide66:  2. 所有家族分析完成後 ,鍵入total stat het指令。 Slide67:  3. 所有家族分析完成後,鍵入total stat het指令, 得到 total lod score 、 NPL score 、 information score。 Slide68:  4. 鍵入 quit 離開 genehunter。

#alleles presentations

Related presentations


Other presentations created by Maurizio

Lecture 13
08. 01. 2008
0 views

Lecture 13

Rona Grimmer Grampian Police
09. 01. 2008
0 views

Rona Grimmer Grampian Police

experimental prob
10. 01. 2008
0 views

experimental prob

The Marble Champ
10. 01. 2008
0 views

The Marble Champ

3 RogerBanks ITC
10. 01. 2008
0 views

3 RogerBanks ITC

Passive Voice and Adverbs
14. 01. 2008
0 views

Passive Voice and Adverbs

powerpoint sample
15. 01. 2008
0 views

powerpoint sample

chapter8
16. 01. 2008
0 views

chapter8

photometry
16. 01. 2008
0 views

photometry

macbethact2
21. 01. 2008
0 views

macbethact2

textile
22. 01. 2008
0 views

textile

BiofuelsJuly2006shor tv2
22. 01. 2008
0 views

BiofuelsJuly2006shor tv2

Stroke 06 001
04. 02. 2008
0 views

Stroke 06 001

Achieving Wireless
04. 02. 2008
0 views

Achieving Wireless

sargam poster
15. 01. 2008
0 views

sargam poster

IST 200701
06. 02. 2008
0 views

IST 200701

PRMs 02 11
18. 02. 2008
0 views

PRMs 02 11

Main presentation
20. 02. 2008
0 views

Main presentation

imyanitov
28. 02. 2008
0 views

imyanitov

0 Introduction
16. 01. 2008
0 views

0 Introduction

leFriant Chris
20. 01. 2008
0 views

leFriant Chris

Lecture 29 SPT DEC
21. 01. 2008
0 views

Lecture 29 SPT DEC

ActuarialProfession0 607
14. 03. 2008
0 views

ActuarialProfession0 607

20 3 Vertebrates form function
19. 03. 2008
0 views

20 3 Vertebrates form function

maynews
03. 04. 2008
0 views

maynews

Lecture22 LargeSatsII
13. 02. 2008
0 views

Lecture22 LargeSatsII

Worldmapper presentation small
28. 03. 2008
0 views

Worldmapper presentation small

safety presentation
05. 02. 2008
0 views

safety presentation

rcdt 061127 04
16. 04. 2008
0 views

rcdt 061127 04

Kickoff Poster Brazil
22. 01. 2008
0 views

Kickoff Poster Brazil

PHED190sportpsycholo gy
17. 04. 2008
0 views

PHED190sportpsycholo gy

speech final Macau
23. 04. 2008
0 views

speech final Macau

The Big Hope march 08
24. 04. 2008
0 views

The Big Hope march 08

PRE DOHA REV
07. 05. 2008
0 views

PRE DOHA REV

2007531101543139
02. 05. 2008
0 views

2007531101543139

BP Idol
10. 01. 2008
0 views

BP Idol

AC2007KhanLecture
05. 03. 2008
0 views

AC2007KhanLecture

LandTour05
25. 02. 2008
0 views

LandTour05

RO CLIL
17. 01. 2008
0 views

RO CLIL

David Brock NASA MSFC
13. 01. 2008
0 views

David Brock NASA MSFC

cfg07 quiz
04. 02. 2008
0 views

cfg07 quiz

133 2005ss elisagaleotti2
11. 01. 2008
0 views

133 2005ss elisagaleotti2