First install the cellmix R package following these instructions
Also, in R, do
install.packages("R.matlab")
The input: Download a GSE using get_GSE_data. Instructions will be written later.
This creates a mat file that contains the following two variables. let it be <INPUT_MAT_FILE>
expression - [genes X samples]
gene_symbol - [genes X 1]
Create a text file markers.txt that has the following format .
Verify that the names of the markers exist in the variable "gene_symbol"."GENE NAME" "Class name"
Some marker files already were set in /cortex/code/cellmix/marker_sets/*.txt
"STMN2" "Neurons"
"BRUNOL4" "Neurons"
"BRUNOL5" "Neurons"
"BRUNOL6" "Neurons"
"SYT1" "Neurons"
"GFAP" "Astrocytes"
"ALDH1L1" "Astrocytes"
"OLIG2" "Oligodendrocytes"
"MBP" "Oligodendrocytes"
Run the following matlab code. The file is also located in /cortex/code/cellmix/run_cellmix.m
addpath('/cortex/code/cellmix')
input_mat_file_name = <INPUT_MAT_FILE>;
output_mat_file_name = <OUTPUT_MAT_FILE>;
marker_file_name = 'markers.txt';
algorithm = 'DECONF'; %'NMF'; 'NMFforb'; 'meanProfile'; 'DSA'; 'DECONF';
run_R_script('do_cellmix', input_mat_file_name, output_mat_file_name, marker_file_name, algorithm);