Problem 1: QTL Data
Estimate genotypic values \(a\) and
\(d\) and predict breeding values for
all animals using the QTL-data given under:
https://charlotte-ngs.github.io/lbgfs2023/data/p1_qtl_1_loci.csv
Your Solution
- Start by re-formatting the paternal and maternal alleles into a
genotype
- Use a linear regression to estimate genotypic values
- Predict breeding values for all animals
Problem 2: Increase Effects of Genotype on Phenotype
Change the phenotypic records in the above given dataset such that
the QTL explains \(50%\) of the genetic
variation when a heritability of \(0.45\) is assumed. It is assumed that the
QTL acts purely additively, hence the genotypic value of the
heterozygotes can be set to \(d=0\).
Show the results as a scatter plot of all phenotypic values for the
QTL genotypes.
Your Solution
- Compute the genetic variance that can be attributed to the QTL based
on the phenotypic variance on the heritability and on the amount of
variation explained by the QTL.
- Read the data and convert the paternal and the maternal alleles to
QTL-genotypes
- Compute allele frequencies
- Compute the genotypic value \(a\)
from the QTL variance
- Add genotypic value to the phenotypes according to the
QTL-Genotype
- Fit regression of new phenotypes on genotypes
- Show the results with plots
Latest Changes: 2023-10-13 09:09:07 (pvr)
LS0tCnRpdGxlOiBMaXZlc3RvY2sgQnJlZWRpbmcgYW5kIEdlbm9taWNzIC0gTm90ZWJvb2sgNQphdXRob3I6IFBldGVyIHZvbiBSb2hyCmRhdGU6ICcyMDIzLTEwLTEzJwpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKcGFyYW1zOgogIGRvY3R5cGU6CiAgICBsYWJlbDogRG9jdW1lbnQgVHlwZQogICAgdmFsdWU6IHNvbHV0aW9uCiAgICBjaG9pY2VzOgogICAgLSBleGVyY2lzZQogICAgLSBzb2x1dGlvbgogICAgLSBub3RlYm9vawogIGlzb25saW5lOgogICAgbGFiZWw6IE9ubGluZSAoeS9uKQogICAgdmFsdWU6IHRydWUKICAgIGNob2ljZXM6CiAgICAtIHRydWUKICAgIC0gZmFsc2UKLS0tCgoKCmBgYHtyIHNldHVwLCBpbmNsdWRlPUZBTFNFfQprbml0cjo6b3B0c19jaHVuayRzZXQoZWNobyA9IFRSVUUpCmBgYAoKCiMjIFByb2JsZW0gMTogUVRMIERhdGEKYGBge3IgcDAxLXNldHVwLCBlY2hvPUZBTFNFfQpzX2RhdGFfcGF0aCA8LSAiaHR0cHM6Ly9jaGFybG90dGUtbmdzLmdpdGh1Yi5pby9sYmdmczIwMjMvZGF0YS9wMV9xdGxfMV9sb2NpLmNzdiIgIApgYGAKCkVzdGltYXRlIGdlbm90eXBpYyB2YWx1ZXMgJGEkIGFuZCAkZCQgYW5kIHByZWRpY3QgYnJlZWRpbmcgdmFsdWVzIGZvciBhbGwgYW5pbWFscyB1c2luZyB0aGUgUVRMLWRhdGEgZ2l2ZW4gdW5kZXI6CgpgYGB7ciwgZWNobz1GQUxTRSwgcmVzdWx0cz0nYXNpcyd9CmNhdChzX2RhdGFfcGF0aCwgIlxuIikKYGBgCgoKIyMjIFlvdXIgU29sdXRpb24KCiogU3RhcnQgYnkgcmUtZm9ybWF0dGluZyB0aGUgcGF0ZXJuYWwgYW5kIG1hdGVybmFsIGFsbGVsZXMgaW50byBhIGdlbm90eXBlCiogVXNlIGEgbGluZWFyIHJlZ3Jlc3Npb24gdG8gZXN0aW1hdGUgZ2Vub3R5cGljIHZhbHVlcwoqIFByZWRpY3QgYnJlZWRpbmcgdmFsdWVzIGZvciBhbGwgYW5pbWFscwoKCgogCgojIyBQcm9ibGVtIDI6IEluY3JlYXNlIEVmZmVjdHMgb2YgR2Vub3R5cGUgb24gUGhlbm90eXBlCmBgYHtyIHAwMi1zZXR1cCwgZWNobz1GQUxTRX0Kc19kYXRhX3BhdGggPC0gImh0dHBzOi8vY2hhcmxvdHRlLW5ncy5naXRodWIuaW8vbGJnZnMyMDIzL2RhdGEvcDFfcXRsXzFfbG9jaS5jc3YiICAKbl9oMl9hbGwgPC0gMC40NQpuX3JhdGlvX3F0bCA8LSAwLjUgCmBgYAoKQ2hhbmdlIHRoZSBwaGVub3R5cGljIHJlY29yZHMgaW4gdGhlIGFib3ZlIGdpdmVuIGRhdGFzZXQgc3VjaCB0aGF0IHRoZSBRVEwgZXhwbGFpbnMgJGByIG5fcmF0aW9fcXRsICogMTAwYCUkIG9mIHRoZSBnZW5ldGljIHZhcmlhdGlvbiB3aGVuIGEgaGVyaXRhYmlsaXR5IG9mICRgciBuX2gyX2FsbGAkIGlzIGFzc3VtZWQuIEl0IGlzIGFzc3VtZWQgdGhhdCB0aGUgUVRMIGFjdHMgcHVyZWx5IGFkZGl0aXZlbHksIGhlbmNlIHRoZSBnZW5vdHlwaWMgdmFsdWUgb2YgdGhlIGhldGVyb3p5Z290ZXMgY2FuIGJlIHNldCB0byAkZD0wJC4gCgpTaG93IHRoZSByZXN1bHRzIGFzIGEgc2NhdHRlciBwbG90IG9mIGFsbCBwaGVub3R5cGljIHZhbHVlcyBmb3IgdGhlIFFUTCBnZW5vdHlwZXMuCgojIyMgWW91ciBTb2x1dGlvbgoKKiBDb21wdXRlIHRoZSBnZW5ldGljIHZhcmlhbmNlIHRoYXQgY2FuIGJlIGF0dHJpYnV0ZWQgdG8gdGhlIFFUTCBiYXNlZCBvbiB0aGUgcGhlbm90eXBpYyB2YXJpYW5jZSBvbiB0aGUgaGVyaXRhYmlsaXR5IGFuZCBvbiB0aGUgYW1vdW50IG9mIHZhcmlhdGlvbiBleHBsYWluZWQgYnkgdGhlIFFUTC4gCiogUmVhZCB0aGUgZGF0YSBhbmQgY29udmVydCB0aGUgcGF0ZXJuYWwgYW5kIHRoZSBtYXRlcm5hbCBhbGxlbGVzIHRvIFFUTC1nZW5vdHlwZXMKKiBDb21wdXRlIGFsbGVsZSBmcmVxdWVuY2llcwoqIENvbXB1dGUgdGhlIGdlbm90eXBpYyB2YWx1ZSAkYSQgZnJvbSB0aGUgUVRMIHZhcmlhbmNlCiogQWRkIGdlbm90eXBpYyB2YWx1ZSB0byB0aGUgcGhlbm90eXBlcyBhY2NvcmRpbmcgdG8gdGhlIFFUTC1HZW5vdHlwZQoqIEZpdCByZWdyZXNzaW9uIG9mIG5ldyBwaGVub3R5cGVzIG9uIGdlbm90eXBlcwoqIFNob3cgdGhlIHJlc3VsdHMgd2l0aCBwbG90cwoKCgoKIAoKCmBgYHtyLCBlY2hvPUZBTFNFLCByZXN1bHRzPSdhc2lzJ30KY2F0KCdcbi0tLVxuXG4gX0xhdGVzdCBDaGFuZ2VzOiAnLCBmb3JtYXQoU3lzLnRpbWUoKSwgJyVZLSVtLSVkICVIOiVNOiVTJyksICcgKCcsIFN5cy5pbmZvKClbJ3VzZXInXSwgJylfXG4nLCBzZXAgPSAnJykKYGBgCiAK