cluster_tissues {EMMIXgene}

Clusters tissues

Description

Clusters tissues

Usage

cluster_tissues(gen, clusters, method = "t", q = 6, G = 2)

Arguments

gen

EMMIXgene object

clusters

mclust object

method

Method for separating tissue classes. Can be either 't' for a univariate mixture of t-distributions on gene cluster means, or 'mfa' for a mixture of factor analyzers.

q

number of factors if using mfa

G

number of components if using mfa

Value

a clustering for each sample (columns) by each group(rows)

Examples

#only run on first 100 genes for speed
alon_sel <- select_genes(alon_data[seq_len(100), ]) 
alon_clust<- cluster_genes(alon_sel,2)
alon_tissue_t<-
   cluster_tissues(alon_sel,alon_clust,method='t')
alon_tissue_mfa<-
   cluster_tissues(alon_sel, alon_clust,method='mfa',q=2,G=2)