spclust {GeneNT}R Documentation

Network constrained "posterior" clustering

Description

This function does network constrained "posterior" clustering based on Floyd-Warshall algorithm, R function allshortestpaths()in R packpage e1071.

Usage

spclust(p, pG2, kG2)

Arguments

p p is the scaling factor that need to be adjusted using biological prior knowledge.
pG2 pG2 is the gene pairs that are screened using the two-stage algorithm based on Pearson correlation statistic.
kG2 kG2 is the gene pairs that are screened using the two-stage algorithm based on Kendall correlation statistic.

Details

This function is written in comparison to the "prior" clustering implemented as priorclust().

Value

This function returns a "posterior" distance matrix that can be used by any distance based clustering software.

Author(s)

Dongxiao Zhu (http://www-personal.umich.edu/~zhud)

References

Zhu, D., Hero, A.O., Qin, Z.S. and Swaroop, A. High throughput screening of co-expressed gene pairs with controlled False Discovery Rate (FDR) and Minimum Acceptable Strength (MAS). Submitted.

See Also

priorclust

Examples

# load GeneNT and GeneTS library
library(GeneTS)
library(GeneNT)
library(e1071)

#EITHER use the internal dataset
data(dat) 
#OR use the following if you want to import external data 
#dat <- read.table("gal.txt", h = T, row.names = 1) 
#Note, data matrix name has to be "dat"
#use (FDR, MAS) criteria (0.2, 0.5) as example to screen gene pairs
#g1 <- corfdrci(0.2, 0.5)  
#pG1 <- g1$pG1
#pG2 contains gene pairs that passed two-stage screening
#pG2 <- g1$pG2  
#use (FDR, MAS) criteria (0.2, 0.5) as example to screen gene pairs 
#g2 <- kendallfdrci(0.2, 0.5) 
#kG1 <- g2$kG1
#kG2 contains gene pairs that passed two-stage screening
#kG2 <- g2$kG2 
#generate Pajek compatible matrix to visualize network
#getBM(pG2, kG2) 
#clustering from network using network constraint clustering
#spclust(3, pG2, kG2) 

[Package GeneNT version 1.2 Index]