Learn R Programming

DOQTL (version 1.8.0)

estimate.cluster.params: Estimate genotype cluster means and variances

Description

Given the X and Y intensity data, perform model based clustering and estimate the genotype state cluster means and variances.

Usage

estimate.cluster.params(founders, data, chr) keep.homozygotes(founders)

Arguments

founders
List containing founder information for non-DO or CC crosses. Not required for DO. When method = intensity: x: Numeric matrix, num.samples x num.snps, with X intensities for all founders and FALSE1s (if available). Sample IDs and SNP IDs must be in rownames and colnames. y: Numeric matrix, num.samples x num.snps, with Y intensities for all founders and FALSE1s (if available). Sample IDs and SNP IDs must be in rownames and colnames. sex: Character vector, containing "M" or F indicating sex. Sample IDs must be in names. code: Character vector containing two letter genotype codes for each founder sample. Sample IDs must be in names.
data
A list with named elements containing the information needed to reconstruct genomes. When method = intensity: x: Numeric matrix, num.samples x num.snps, with X intensities for all samples. Sample IDs and SNP IDs must be in rownames and colnames. y: Numeric matrix, num.samples x num.snps, with Y intensities for all samples. Sample IDs and SNP IDs must be in rownames and colnames. sex: Character vector, containing "M" or F indicating sex. Sample IDs must be in names. gen: Character matrix containing the generation of DO outbreeding for each sample. FALSEor the DO, this should be "DO" followed by a number with no space between them. FALSEor CC mice, this should be CC. Sample IDs must be in names.
chr
Character containing the current chromosome.

Value

List containing two elements:
r.t.means
Three dimensional array containing rho and theta genotype cluster means.
r.t.covars
Three dimensional array containing rho and theta genotype cluster variances.

Details

At each marker, use mclust to perform model based clustering on all of the data and get estimates of the means and variances for each cluster. Then assign each of the 36 genotype states to the nearest founder cluster. keep.homozygotes is an internal helper function.

See Also

hmm.intensity

Examples

Run this code
  ## Not run:  estimate.cluster.params(founders, data, chr) 

Run the code above in your browser using DataLab