estimate.cluster.params {DOQTL} | R Documentation

## 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.

### 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.

### 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.

### Author(s)

Daniel Gatti

### See Also

`hmm.intensity`

### Examples

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

[Package DOQTL version 1.18.0 Index]