svr_functionally_coupled

Get functionally_coupled neighbors (neighbors that tend to co-occur).

------

Example:

    svr_all_features 3702.1 peg | svr_functionall_coupled

would produce a 3-column table. The first column would contain PEG IDs for genes occurring in genome 3702.1, the second would give the functional coupling score, and the third would give the PEG that is functionally coupled.

------

The standard input should be a tab-separated table (i.e., each line is a tab-separated set of fields). Normally, the last field in each line would contain the PEG for which functions are being requested. If some other column contains the PEGs, use

    -c N

where N is the column (from 1) that contains the PEG in each case.

This is a pipe command. The input is taken from the standard input, and the output is to the standard output.

Command-Line Options

-c Column

This is used only if the column containing PEGs is not the last.

-max MaxReturned [default 1000000]

This is used to restrict the number of results displayed for a single incoming line (not a restriction on the number of tatal lines)

-n MinSc [default 15]

Output Format

The standard output is a tab-delimited file. It consists of lines from the input file that are for PEGs that are functionally coupled. The lines will have two appended columns: the functional coupling score (number of distinct OTUs in which the pair co-occur), and the functionally coupled PEG.