Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error in the clustering step of G1/2 cells in scRT.infer() #12

Open
annahyan opened this issue Nov 13, 2024 · 0 comments
Open

Error in the clustering step of G1/2 cells in scRT.infer() #12

annahyan opened this issue Nov 13, 2024 · 0 comments

Comments

@annahyan
Copy link

Hi,

I am trying PERT inference on my data without providing the clone_id and I get the following error when running the scrt.infer function:

# run inference using PERT
... cn_s_with_scrt, supp_s_output, cn_g_with_scrt, supp_g_output = filtered_scrt.infer(level='pyro')
 20122010 trying with max k=20
 20122011 trying with k=2
 20122122 trying with k=3
 20122213 trying with k=4
 20122317 trying with k=5
 20122433 trying with k=6
 20122551 trying with k=7
 20122673 trying with k=8
 20122830 trying with k=9
 20123017 trying with k=10
 20123183 trying with k=11
 20123364 trying with k=12
 20123543 trying with k=13
 20123754 trying with k=14
 20123965 trying with k=15
 20124194 trying with k=16
 20124428 trying with k=17
 20124688 trying with k=18
 20124929 trying with k=19
 20125170 trying with k=20
 20125424 selected k=17
Traceback (most recent call last):
  File "<stdin>", line 2, in <module>
  File "/app/scdna_replication_tools/infer_scRT.py", line 122, in infer
    self.cn_s, supp_s_out_df, cn_g1_out, supp_g1_out_df = self.infer_pert_model()
  File "/app/scdna_replication_tools/infer_scRT.py", line 148, in infer_pert_model
    cell_col=self.cell_col, chr_col=self.chr_col, start_col=self.start_col)
  File "/app/scdna_replication_tools/assign_s_to_clones.py", line 73, in assign_s_to_clones
    copy_corrs = clone_correlations(clone_df, temp_cell_cn, col_name)
  File "/app/scdna_replication_tools/assign_s_to_clones.py", line 43, in clone_correlations
    r, pval = pearsonr(merged_df[col_name], clone_cn)
  File "/usr/local/lib/python3.7/site-packages/scipy/stats/stats.py", line 3906, in pearsonr
    raise ValueError('x and y must have length at least 2.')
ValueError: x and y must have length at least 2.

I'd appreciate your help.

Best

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant