motifclstest: Find motifs enriched in clustering

Description Usage Arguments Details Value Author(s) Examples

View source: R/DynaMO.R

Description

Calculate enrichment levels of motifs in clusters generated by K-means clustering based on read counts at predicted motif sites across different samples.

Usage

1
motifclstest(motifcls, motifname)

Arguments

motifcls

A matrix object with each row representing a motif site and each column representing read counts from one sample. The last three columns are the ids of clusters, the ids of motifs and the ids of motif sites.

motifname

A vector of ids of motifs. For example, if x motifs are examined, motifname is 1:x

Details

This function counts the numbers of motif sites in each cluster for each motif and calculates the fold enrichments and p values for each motif in each cluster and outputs in a matrix where each row represents a motif.

Value

A table is written in the "motif_cluster_enrichment.txt".

Author(s)

Zheng Kuang

Examples

1
2
3
motiffdrread=read.table(system.file("extdata","motif_filter_reads_cluster.txt",
package="DynaMO",mustWork=TRUE),header=TRUE,sep="\t")
motiffdrclstest=motifclstest(motiffdrread,1:5)

spo111/DynaMO documentation built on May 30, 2019, 7:59 a.m.