condense.sanger.snps: Create an HDF5 file with the unique SNP patterns between each...
In DOQTL: Genotyping and QTL Mapping in DO Mice

Description Usage Arguments Details Value Author(s) See Also Examples

Given a set of markers, a SNP VCF and a set of founder strains that are in the Sanger Mouse Genomes, create an HDF5 file that contains one object for each pair of markers.

1	condense.sanger.snps = function(markers, snp.file, strains, hdf.file, ncl = 1)

`markers`	Data.frame containing the markers locations. There must be at least three columns containing the marker name, the chromosome and the Mb position. Other columns are ignored.
`snp.file`	Character string containing the full path to the Sanger SNP VCF file. Get the file from ftp://ftp-mouse.sanger.ac.uk/.
`strains`	Character vector containing the names of strains to extract from the Sanger VCF file. Obtain the names using `get.vcf.strains`.
`hdf.file`	Character string containing the name of the HDF5 file. NOTE: this function will not overwrite an existing file.
`ncl`	Integer that is the number of nodes to use for parallel execution. Default = 1.

The markers are usually MUGA series markers, but they can be any set of genomic positions. The markers are broken up by chromosome and written out to and HDF5 file in which each data object is named using the proximal and distal marker named separated by and "_". At the beginning of the chromosome, the object is named using "start" and the first marker separated by an "_". At the end of the chromosome, the object is named using the last marker and "end" separated by an "_". Each object is a list containing three elements: sdps: Numeric matrix containing the unique strain distribution patterns in the founder strains. map: Numeric vector indicating the index of the SNP at which each SDP occurs with the current pair of markers. snps: Data.frame containing SNP names, positions, reference and alternate alleles. NOTE: We assume (incorrectly) that all SNPs are bimorphic. SNPs are coded as 0 for "equal to reference" and 1 for "not equal to reference".

Writes out and HDF5 file.

Daniel Gatti

read.vcf, read.vcf

  ## Not run:  
    load(url("ftp://ftp.jax.org/MUGA/muga_snps.Rdata"))
	vcf.file = "mgp.v4.snps.dbSNP.vcf.gz"
    condense.sanger.snps(markers = muga_snps, snp.file = vcf.file, strains = get.vcf.strains(vcf.file), hdf.file = "output.h5") 
  
## End(Not run)