reduceGeno: Reduction Dimension of Genotype

In FRGEpistasis: Epistasis Analysis for Quantitative Traits by Functional Regression Model

Description

Reduce Dimension of Genotype using Functional Regression Model

Usage

reduceGeno(wDir, pheno, gnoFiles, mapFiles, gLst, rng)

Arguments

wDir	The dataset directory. If the dataset is in the working directory, wDir is ".".
pheno	It is a matrix with two columns. One column is the individual ID and the other is the phenotype. The phenotype can be quantitative trait or binary trait.
gnoFiles	The vector of genotype file names. It contains the genotype file names indicating where to read the genotype files.
mapFiles	The vector of SNP genetic map file names. It contains the map file names indicating where to read the genetic map files.
gLst	Gene annotation which includes gene name, chromosome, start position and end position.
rng	A numeric value which represents gene region extensible scope.

Details

This function reduces the dimension of genotypes of all the genes with Fourier expansion. In real dataset the genotypes on different chromosome are always organized into different files. And each genotype file is very large. This function processes the genotype files in turns. The reduced genotype data (the expansion data) of all the chromosomes are combined together. The expansion data and other information are organized into a list. During Fourier expansion, the physical position information of the SNPs are used. This is one of merits of our method.

Value

Return a list that includes reduced genotype data, gene names,chromosome information, start index and end index of each gene.

Author(s)

Futao Zhang

Examples

wDir <-paste(system.file("extdata", package="FRGEpistasis"),"/",sep="")
gnoFiles<-read.table(system.file("extdata", "list_geno.txt", package="FRGEpistasis"))
mapFiles<-read.table(system.file("extdata", "list_map.txt", package="FRGEpistasis"))
phenoInfo <- read.csv(system.file("extdata", "phenotype.csv", package="FRGEpistasis"),header=TRUE)
gLst<-read.csv(system.file("extdata", "gene.list.csv", package="FRGEpistasis"))
rng=0
gInfo=reduceGeno(wDir,phenoInfo,gnoFiles,mapFiles,gLst,rng)

FRGEpistasis documentation built on Nov. 8, 2020, 5:51 p.m.