medoQtc: Find medoids of QTCs

In QTCAT/qtcat: Quantitative Trait Cluster Association Test

Description

Find a medoid of for each quantitative trait cluster (QTC).

Usage

medoQtc(object, geno, alpha = 0.05, min.absCor = 0.05)

Arguments

object	an object of class qtcatHit.
geno	an object of class qtcatGeno.
alpha	an alpha level for significance estimation.
min.absCor	minimum absolute value of correlation considered.

Examples

# If you want to run the examples, use:
# example(medoQtc, run.dontrun = TRUE)
## Not run: 
# files containing example data for SNP data and the phenotype
pfile <- system.file("extdata/phenodata.csv", package = "qtcat")
gfile <- system.file("extdata/snpdata.csv", package = "qtcat")
pdat <- read.csv(pfile, header = TRUE)
snp <- read.snpData(gfile, sep = ",")
clust <- qtcatClust(snp)
geno <- qtcatGeno(snp, clust)
pheno <- qtcatPheno(names = pdat[, 1],
                    pheno = pdat[, 2],
                    covariates = model.matrix(~ pdat[, 3]))
fitted <- qtcatHit(pheno, geno)

# QTC medoids
medo <- medoQtc(fitted, geno)

## End(Not run)

QTCAT/qtcat documentation built on April 20, 2021, 11:20 p.m.