phenotype | R Documentation |
The function calculates phenotying for the input haplotype object using the QTN supplied effects.
phenotype(pop,qtn_effect,tgv_only,vr,parental_genotypes)
pop |
Required: a population |
qtn_effect |
Required: qtn_effect data frame |
tgv_only |
Required: if TRUE, only calculated breeding value will be export, without environmental effects, default is FALSE |
vr |
Required: Vr value |
parental_genotypes |
Required: parental_genotypes object |
Vr value
gff = "crossword_processing/peanut.gff"
chr_stat = "crossword_processing/chr_siz.txt"
chr_length = "crossword_processing/chr_len.txt"
gen2phy = physical2genomic(gff,chr_stat,chr_length,200000)
homo = TRUE
input_pheno = "crossword_processing/pheno.txt"
parental_genotypes = get_parental_genotypes(input,gen2phy,homo)
biased_selection = TRUE
min_qtn_freq = 0
dominant=FALSE
effect_distribution = "equal"
pop3 = create_families(pop=pop2,S=15,id="pop3",chr_length=chr_length)
qtn_effect = random_qtn_assign(qtn = 5,gen2phy=gen2phy,biased_selection=biased_selection,parental_genotypes=parental_genotypes,min_qtn_freq=min_qtn_freq,dominant=min_qtn_freq,effect_distribution=effect_distribution)
h2 = 0.6
tgv_only = FALSE
vr1 = get_vr(qtn_effect = qtn_effect,h2 = 0.6 ,parental_genotypes=parental_genotypes,dominant=TRUE,heritability_mode="absolute")
phen1 = get_vr(pop=pop3,qtn_effect=qtn_effect,vr=vr1)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.