Esempi di statistica descrittiva e inferenziale/Analisi del proprio genoma: differenze tra le versioni

Contenuto cancellato Contenuto aggiunto
Riga 113:
filter(rsid=='rs1799971')
</syntaxhighlight>
 
 
rsid chromosome position genotype
Line 132 ⟶ 133:
 
DISEASE.TRAIT STRONGEST.SNP.RISK.ALLELE genotype
1 Skin reflectance (Melanin index) rs2272756-A AG
2 Pancreatic cancer rs13303010-G AA
3 Body mass index rs3934834-G CT
4 Urate levels rs9442380-T CC
5 Epithelial ovarian cancer rs9442387-? CT
6 Blood protein levels rs3766186-C AC
STUDY
1 A genome-wide association study of skin and iris pigmentation among individuals of South Asian ancestry.
Line 162 ⟶ 163:
head()
</syntaxhighlight>
 
DISEASE.TRAIT STRONGEST.SNP.RISK.ALLELE genotype
1 Skin reflectance (Melanin index) rs2272756-A AG
6 Blood protein levels rs3766186-C AC
11 Serum alkaline phosphatase levels rs1123571-A AG
16 Body mass index rs7535528-G GG
18 Autoimmune thyroid disease rs2234167-A AG
19 Medication use (thyroid preparations) rs2234167-A AG
 
 
Si estraggono dal dataset le righe relative agli alleli di rischio contenuti nel genotipo del signore 2 volte. Ad esempio se l'allele è A e nel genotipo si trova AA.
Line 176 ⟶ 186:
}
 
v<-df1 %>%
select(DISEASE.TRAIT, STRONGEST.SNP.RISK.ALLELE, genotype,STUDY)
</syntaxhighlight>
 
DISEASE.TRAIT STRONGEST.SNP.RISK.ALLELE genotype
16 Body mass index rs7535528-G GG
23 Body mass index rs6667605-T TT
35 Serum total cholesterol levels rs1456465-G GG
36 Mean corpuscular volume rs1569419-C CC
37 Plateletcrit rs1569419-C CC
38 Platelet distribution width rs1569419-C CC
 
<syntaxhighlight lang="rsplus">