WCQ (version 0.2)

wcq.fdrc: WCQ QTL detection with false discovery rate control

Description

Performs WCQ with optional false discovery rate control.

Usage

wcq.fdrc(marker_data, trait_data, alleles, fdrc.method = "none")

Arguments

marker_data

m x n matrix/array containing marker score data; m is the number of markers, n is the number of samples

trait_data

t x n matrix/array containing trait value data; t is the number of traits, n is the number of samples

alleles

vector containing the marker scores representing the alleles; currently supports diploids only; third entry assumed to be score for heterozygous marker

fdrc.method

string value representing method to be used for false discovery rate control; uses the p.adjust method in the R default library; default value is "none"

Value

adjusted_pval_matrix

m x t array/matrix containing adjusted p-values of each marker as being a potential QTL for a trait; if fdrc.method is set to "none", returns the pval_matrix computed using the wcq function

See Also

p.adjust

Examples

data(sample.markers)
data(sample.traits)
alleles <- c(1,2,3)
wcq.fdrc(sample.markers, sample.traits, alleles, fdrc.method="bonferroni")