recreate regulatory matrix BMDC 0hr

[Wybe]

Dear Pertub-seq team,

I'm trying to reproduce some of your results using the data provided on GEO. 

Specifically I'm trying to recreate the regulatory matrix provided in the supplementary material.

This one to be exact: Regulatory matrices for BMDC without modelling cells states, BMDC 0hr. 

I have created the expression matrix Y from  GSM2396857_dc_0hr.mtx.txt.gz.

I have created the covariate matrix from GSM2396857_dc_0hr_cbc_gbc_dict.csv.gz, using the function provided in mimosca.py. X = dict2X(dict, cellNamesList), where dict is a dictionary object created from the given file.

However when running run_model_bycol() with all column names of the covariate matrix as input

I get a different regulatory matrix than the one you provided in the supplementary material.

Any ideas on what might be going wrong? 

Also where to columns BF-BN come from in the regulatory matrix for BMDC 0hr without cell states? 

If you could provide me with more data that would be great. 

Cheers,

Wybe