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 when a gene is always zero-expressed #38

Open
YazdanZ opened this issue Jul 22, 2022 · 0 comments
Open

Error when a gene is always zero-expressed #38

YazdanZ opened this issue Jul 22, 2022 · 0 comments

Comments

@YazdanZ
Copy link

YazdanZ commented Jul 22, 2022

I am running PIDC using the BEELINE pipeline (Julia v1.1.0).
I get the following error when a gene is zero-expressed for every cell:

6d6612190f7e:/data grnbeeline/pidc:base /bin/sh -c "time -v -o data/outputs/defpos/GSD/PIDC/time.txt julia runPIDC.jl data/inputs/defpos/GSD/PIDC/ExpressionData.csv data/outputs/defpos/GSD/PIDC/outFile.txt "
ERROR: LoadError: MethodError: no method matching setindex_shape_check(::Int64, ::Int64)
Closest candidates are:
 setindex_shape_check(!Matched::AbstractArray{#s72,1} where #s72, ::Integer) at indices.jl:218
 setindex_shape_check(!Matched::AbstractArray{#s72,1} where #s72, ::Integer, !Matched::Integer) at indices.jl:221
 setindex_shape_check(!Matched::AbstractArray{#s72,2} where #s72, ::Integer, !Matched::Integer) at indices.jl:225
 ...
Stacktrace:
[1] macro expansion at ./multidimensional.jl:694 [inlined]
[2] _unsafe_setindex!(::IndexLinear, ::Array{Int64,1}, ::Int64, ::UnitRange{Int64}) at ./multidimensional.jl:689
[3] _setindex! at ./multidimensional.jl:684 [inlined]
[4] setindex! at ./abstractarray.jl:1020 [inlined]
[5] get_bin_ids!(::Array{Float64,1}, ::String, ::Int64, ::Array{Int64,1}) at /root/.julia/packages/InformationMeasures/fdfJk/src/Discretization.jl:111
[6] Node(::Array{Any,2}, ::String, ::String, ::Int64) at /root/.julia/packages/NetworkInference/z8pnG/src/common.jl:32
[7] #get_nodes#15(::Bool, ::String, ::String, ::Int64, ::Function, ::String) at /root/.julia/packages/NetworkInference/z8pnG/src/infer_network.jl:35
[8] get_nodes(::String) at /root/.julia/packages/NetworkInference/z8pnG/src/infer_network.jl:26
[9] top-level scope at util.jl:156
[10] include at ./boot.jl:326 [inlined]
[11] include_relative(::Module, ::String) at ./loading.jl:1038
[12] include(::Module, ::String) at ./sysimg.jl:29
[13] exec_options(::Base.JLOptions) at ./client.jl:267
[14] _start() at ./client.jl:436

I don't think it's a problem with BEELINE but with PIDC since this is also what has been reported in #29.
As soon as I change the expression of the zero-expressed gene to a value other than 0 in a single cell, the problem goes away.

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