DAWN analysis

The users can investigate the relationship between categories and identify the specific type of categories that are clustered within the network of categories of interest.

  • -e, –eig_vector: Eigen vector file. This is the output file from calculation of effective number of tests. The file name must have pattern *eig_vecs*.zarr.

  • -c, –corr_mat: Category correlation matrix file. This is the output file from categorization. The file name must have pattern *correlation_matrix*.pkl.

  • -P, –permut_test: Permutation test file. This is the output file from burden test. The file name must have pattern *permutation_test*.txt.gz.

  • -o_dir, –output_directory: Path to the directory where the output files will be saved. By default, outputs will be saved at $CWAS_WORKSPACE.

  • -r, –range: Range (i.e., (start,end)) to find optimal K for k-means clustering. It must contain two integers that are comma-separated. The first integer refers to the start number and must be above 1. The second integer refers to the end.

  • -k, –k_val: K for K-means clustering. With this argument, users can determine K manually. -r and -k arguments are mutually exclusive. If -k is given, -r will be ignored.

  • -s, –seed: Seed value for t-SNE. Same seed will generate same results for the same inputs.

  • -T, –tsen_method: Gradient calculation algorithm for t-SNE, which is used in TSNE of sklearn. If the dataset is large, ‘barnes_hut’ is recommended. By default, exact.

  • -t, –tag: Tag used for the name of the output files. By default, None.

  • -c_count, –cat_count: Path of the categories counts file from burden test.

  • -C, –count_threshold: The treshold of variant (or sample) counts. The least amount of variants a category should have.

  • -R, –corr_threshold: The threshold of correlation values between clusters. Computed by the mean value of correlation values of categories within a cluster.

  • -S, –size_threshold: The threshold of the number of categories per cluster. The least amount of categories a cluster should have.

  • -p, –num_proc: Number of worker processes that will be used for the DAWN analysis. By default, 1.

cwas dawn -e INPUT_EIG_VEC \
-c INPUT_CORR_MATRIX \
-P INPUT_PERMUATION_RESULT \
-o_dir OUTPUT_DIR \
-r 2,500 \
-s 123 \
-t test \
-c_count CATEGORY_COUNTS.txt \
-C 20 \
-R 0.7 \
-S 2 \
-p 8