generate_genodata: A Function That Generates a Genome.

Description Usage Arguments Examples

Description

This function generates a genome around significant parental alleles.

Usage

1
generate_genodata(alleles, half_chr_width = 100, teta = 0.05)

Arguments

alleles

A vector that describes significant parental alleles (those that will control your phenotype

half_chr_width

An integer that defines the length of the (half (chromosome size minus 1)).

teta

An interger that specify recombination factor to use.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
alleles = sample(c("BY","RM"), 25, replace=TRUE)
print(alleles)
genodata = generate_genodata(alleles, half_chr_width=50)
print(genodata)
tetas = compute_teta(genodata)
# matlab::imagesc(t(as.matrix(genodata)), xlab="markers", 
#   ylab="segrgants", main="Parental origin of segregant genomes")
plot(density(tetas), main=paste("Recombination Fraction Distribution"), 
  xlab=paste("mean(teta)=", signif(mean(tetas),3), sep=""))
abline(v=mean(tetas), lty=2)

fchuffar/ptldata documentation built on May 16, 2019, 12:07 p.m.