Publication link ———————————————— https://www.nature.com/articles/s41598-017-04070-4
*This file is part of the CCS package for biclustering analysis
Condition-dependent Correlation Subgroup (CCS) is a biclustering algorithm for comprehensive discovery of functionally coherent biclusters from large-scale gene expression data. For the details of the CCS algorithm, see our paper entitled “A GPU-accelerated algorithm for biclustering analysis and detection of condition-dependent coexpression network modules”. The algorithm is implemented in C. See the steps in the section A and B for compilation and execution of the C code. The structure of the CCS algorithm is particularly suitable for parallel computing. A CUDA C based GPGPU computing code is included here as a parallel version of the algorithm. You need a programmable GPU card and CUDA C complier for compilation and execution of our code. Follow the steps in C and D. The performance of CCS was tested on synthetic and real gene expression datasets. We also showed that there is an equivalence between CCS biclusters and condition-dependent co-expression network modules. The related Python codes (Python 2.7 or higher) are available in the "python_utility" directory. Synthetic and real gene expression datasets, and CCS biclustering results are available in the "Results" directory.
- Change your current directory to "src". $cd src
- Type "make" to create executable "ccsbc" in the home directory $make
- Back to parent directory cd ../
Type the following commands in Linux: ./ccs -t [correlation threshold] -i [input file] -o [output file] Parameters: -t [0-1.0]: Specify correlation threshold "theta" between 0 to 1. Recommended value is 0.8. -i [input_data_file] -o [output_data_file] Example: ./ccs -t 0.8 -i ./Results/Synthetic_data_results/Data/Data_Constant_100_1_bicluster.txt -o ./Results/Output.txt Additional parameters: -m [1 - number of gene/rows in the data matrix]: Set the number of base gene that are to be considered for forming biclusters. Default value is 1000 or maximum number of genes when that is less than 1000. Example: ./ccs -t 0.8 -i ./Results/Synthetic_data_results/Data/Data_Constant_100_1_bicluster.txt -o ./Results/Output.txt -m 90 -g [0.0 - 100.0]: Minimum gene set overlap required for merging the overlapped biclusters. Default value is 100.0 for 100% overlap. Example: ./ccs -t 0.8 -i ./Results/Synthetic_data_results/Data/Data_Constant_100_1_bicluster.txt -o ./Results/Output.txt -g 50.0 -p [0/1]: Set the output format. Default is 0. 0 - Print output in 3 rows.
Row 1: Number_of_rows[\t]Number_of_Columns[\t]Score
Row 2: Gene_name_1[b]Gene_name_2[b] ...
Row 3: Sample_name_1[b]Sample_name_2[b] ...
1 - Print output in 2 rows (Bibench supported format).
Row 1: Row_index_1[b]Row_index_2[b] ...
Row 2: Column_index_1[b]Column_index_2[b] ...
Example: ./ccs -t 0.9 -i ./Results/Synthetic_data_results/Data/Data_Constant_100_1_bicluster.txt -o ./Results/Output_standard.txt -m 50 -p 1 -g 100.0
*Note that a CUDA supported GPU card and CUDA C compiler is required.
- Change your current directory to "CUDA_C". $cd CUDA_C
- Type following in the linux command line $nvcc ./src/ccs.cu -lm -o ccs_cuda
Type following commands in the Linux: ./ccs_cuda -t [correlation threshold] -i [input file] -o [output file] Example: ./ccs_cuda -t 0.9 -i ../Synthetic_data_results/Data/Data_Constant_100_1_bicluster.txt -o ./Output.txt -m 50 -p 1 -g 100.0
• Anindya Bhattacharya, [email protected] • Yan Cui, [email protected]
If you have comments or questions, or if you would like to contribute to the further development of CCS, please send us an email at [email protected] and [email protected]
This projected is licensed under the terms of the GNU General Public License v3.0.