makeoutput {hopach} | R Documentation

## function to write a text file with hopach output

### Description

The function `makeoutput` takes a data matrix and corresponding `hopach` clustering output, plus possibly bootstrap resampling output, and makes a table summarizing the clustering results. The table is written to a tab delimited text file.

### Usage

    makeoutput(data, hopachobj, bootobj = NULL, file = "HOPACH.out",
    gene.names = NULL)

### Arguments

`data` | data matrix, data frame or exprSet of gene expression measurements. Typically, each column corresponds to an array, and each row corresponds to a gene. For clustering arrays, the arrays appear in the rows and the genes in the columns. All values must be numeric. Missing values are ignored.
`hopachobj` | output of the `hopach` function.
`bootobj` | optional output of `boothopach` or `bootmedoids` - a matrix of bootstrap estimated cluster membership probabilities, with a row for each row in `data` and a column for each cluster.
`file` | filename for the table produced.
`gene.names` | optional names or annotations for the genes (arrays), which can be different from the row names of `data`

### Details

The output table contains information about the rows of `data`: Index, UID and Name; the main cluster results: Cluster Number, Cluster Label (from `hopach`), and Cluster Level Order (the ordering of the elements in the level of the `hopach` tree at which the main clusters were identified); and the final level of the tree: Final Label, and Final Order (the ordering of elements in the final level of the tree). Sorting this table on Index results in the rows having the same order as in data. Sorting on Cluster Level Order results in the rows being ordered by cluster, and then within cluster based on the value of the argument `ord` to `hopach` (default is distance to the medoid). Sorting on Final Level Order results in the rows being ordered as in the leaves of the `hopach` tree, where clusters are still ordered and elements near each other in the ordering will have small pairwise distances.

### Value

The function `makeoutput` has no value. It writes a tab delimited text file to the current working directory.

### Note

Thank you to Karen Vranizan <vranizan@uclink.berkeley.edu> for helping to write this function.

### Author(s)

Katherine S. Pollard <kpollard@soe.ucsc.edu>

### References

van der Laan, M.J. and Pollard, K.S. A new algorithm for hybrid hierarchical clustering with visualization and the bootstrap. Journal of Statistical Planning and Inference, 2003, 117, pp. 275-303.

http://www.stat.berkeley.edu/~laan/Research/Research_subpages/Papers/hopach.pdf

### See Also

`hopach`, `boothopach`, `bootmedoids`

### Examples

    #25 variables from two groups with 3 observations per variable
    mydata<-rbind(cbind(rnorm(10,0,0.5),rnorm(10,0,0.5),rnorm(10,0,0.5)),cbind(rnorm(15,5,0.5),rnorm(15,5,0.5),rnorm(15,5,0.5)))
    dimnames(mydata)<-list(paste("Var",1:25,sep=""),paste("Exp",1:3,sep=""))
    mydist<-distancematrix(mydata,d="cosangle") #compute the distance matrix.

    #clusters and final tree
    clustresult<-hopach(mydata,dmat=mydist)

    #bootstrap resampling
    myobj<-boothopach(mydata,clustresult)

    #write output file
    makeoutput(mydata,clustresult,myobj)

[Package hopach version 1.6.0 Index]