meanIntensityByScanChrom: Calculate Means & Standard Deviations of Intensities
In amstilp/GWASTools: Tools for Genome Wide Association Studies

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

Function to calculate the mean and standard deviation of the intensity for each chromosome for each scan.

1 2	meanIntensityByScanChrom(intenData, vars = c("X", "Y"), snp.exclude = NULL, verbose = TRUE)

`intenData`	`IntensityData` object. Chromosomes are expected to be in contiguous blocks.
`vars`	Character vector with the names of one or two intensity variables.
`snp.exclude`	An integer vector containing SNPs to be excluded.
`verbose`	Logical value specifying whether to show progress information.

The names of two intensity variables in intenData may be supplied. If two variables are given, the mean of their sum is computed as well. The default is to compute the mean and standard deviation for X and Y intensity.

A list with two components for each variable in "vars": 'mean.var' and 'sd.var'. If two variables are given, the first two elements of the list will be mean and sd for the sum of the intensity variables:

`mean.intensity`	A matrix with one row per scan and one column per chromosome containing the means of the summed intensity values for each scan and chromosome.
`sd.intensity`	A matrix with one row per scan and one column per chromosome containing the standard deviations of the summed intensity values for each scan and chromosome.
`mean.var`	A matrix with one row per scan and one column per chromosome containing the means of the intensity values for each scan and chromosome.
`sd.var`	A matrix with one row per scan and one column per chromosome containing the standard deviations of the intensity values for each scan and chromosome.

Cathy Laurie

IntensityData, mean, sd

file <- system.file("extdata", "illumina_qxy.gds", package="GWASdata")
gds <- GdsIntensityReader(file)
intenData <- IntensityData(gds)

meanInten <- meanIntensityByScanChrom(intenData)
close(intenData)