i. proc allele (sas/genetics) ii. single snp analysis iii. tests for multi-allelic markers iv....
TRANSCRIPT
![Page 1: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/1.jpg)
I. Proc Allele (SAS/Genetics)
II. Single SNP analysis
III. Tests for Multi-allelic Markers
IV. Haplotype testsi. Macro %HAPPY (http://www.hsph.harvard.edu/faculty/kraft/soft.htm)
ii. Hapstat (http://www.bios.unc.edu/~lin/hapstat/)
Lab 2
Yu-Chun Jean [email protected] Rm. 200
![Page 2: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/2.jpg)
Sample DatasetYahoo.dat: a comma delimited file Case-control with 1680 subjects, 16 SNPs
![Page 3: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/3.jpg)
![Page 4: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/4.jpg)
![Page 5: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/5.jpg)
![Page 6: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/6.jpg)
![Page 7: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/7.jpg)
![Page 8: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/8.jpg)
![Page 9: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/9.jpg)
Make sure you're modeling correct outcome
![Page 10: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/10.jpg)
Dominant Model
![Page 11: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/11.jpg)
Codominant Model
![Page 12: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/12.jpg)
Additive Model
![Page 13: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/13.jpg)
![Page 14: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/14.jpg)
proc multtest pdata=outpeas bon hom fdr;run;
BY threshold is
.05/3.38~.015
![Page 15: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/15.jpg)
![Page 16: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/16.jpg)
![Page 17: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/17.jpg)
![Page 18: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/18.jpg)
Most General Model
![Page 19: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/19.jpg)
Model -2 log L
Additivez1 z2 z3 / ref=z0 711.218
z0 z2 z3 / ref=z1 711.218
Dominantd1 d2 d3 / ref=d0 711.797d0 d2 d3 / ref=d1 710.719
Additive Model and Dominant Model
![Page 20: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/20.jpg)
![Page 21: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/21.jpg)
%happy(indsn=yahoo, keep=a7 a8 a19 a20, style=MAR, outadd=yahadd, range=.05);
SNPs (two allele vars each)
Tells happy input format
Where should I save additive scores?
How rare is rare? Rare haplotypes are lumped in "other" category.
![Page 22: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/22.jpg)
![Page 23: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/23.jpg)
proc print data=origfreq;run;
Dataset containing ordered hap freqs and reference numbers
![Page 24: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/24.jpg)
![Page 25: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/25.jpg)
/* set up analysis data set */proc sort data=yahadd; by id;proc sort data=yahoo; by id;data UseMe; merge yahoo yahadd; by id;run;
/* do global haplotype test */proc logistic data=UseMe; model d = z2 z3;run;
/* test haplotype 2-4 specifically */proc logistic data=UseMe; model d = z2;run;
![Page 26: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/26.jpg)
![Page 27: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/27.jpg)
![Page 28: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/28.jpg)
![Page 29: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/29.jpg)
![Page 30: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/30.jpg)
![Page 31: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/31.jpg)
![Page 32: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/32.jpg)
![Page 33: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/33.jpg)
![Page 34: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/34.jpg)
![Page 35: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/35.jpg)
● Haplotype 11(z2) & Haplotype 01(z3)Log-Likelihood -2167.7410
● Haplotype 11(z2)Log-Likelihood -2168.0977
● Null modelLog-Likelihood -2174.8674
LRT = 14.25
LRT =13.54
![Page 36: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/36.jpg)
![Page 37: I. Proc Allele (SAS/Genetics) II. Single SNP analysis III. Tests for Multi-allelic Markers IV. Haplotype tests i. Macro %HAPPY (](https://reader035.vdocument.in/reader035/viewer/2022062600/5a4d1bc27f8b9ab0599d3a4e/html5/thumbnails/37.jpg)