uhlerlab/causaldag

ValueError: On entry to DLASCL parameter number 4 had an illegal value

Opened this issue · 0 comments

Thanks for this very promising package. I am trying to apply this tool on single cell RNA seq data, of once cell type that has been sorted out. I have treatment vs control samples. However on following this tutorial, I encounter the following error. I do not have any NAs in the data set. Please how could I fix this?

!pip install --user causaldag
from causaldag import dci, dci_stability_selection
from causaldag import get_directed_and_undirected_edges
import os
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import networkx as nx
from sklearn.preprocessing import StandardScaler

# Matplotlib command for Jupyter notebooks only
%matplotlib inline


path_in = '/home/power/docker_scratch'
gene_exp_diff_1 = pd.read_csv(os.path.join(path_in, 'log_norm_plus1_counts_control.csv'), index_col=0)
gene_exp_diff_2 = pd.read_csv(os.path.join(path_in, 'log_norm_plus1_counts_treated.csv'), index_col=0)
gene_names = gene_exp_diff_1.columns.values

scaler = StandardScaler(with_std=False)
X1 = scaler.fit_transform(gene_exp_diff_1.values)
X2 = scaler.fit_transform(gene_exp_diff_2.values)

d_dag = dci(X1, X2, alpha_ug = 0.001, alpha_skeleton=0.5, alpha_orient=0.001, max_set_size=3, verbose=1)

ValueError: On entry to DLASCL parameter number 4 had an illegal value