guanjue / IDEAS_2018

Jointly characterizing epigenetic dynamics across multiple cell types
MIT License
6 stars 6 forks source link

error report #4

Open KailiBio opened 6 years ago

KailiBio commented 6 years ago

Hi. I got a new error repot show as below.

Error in apply(statemean, 1, range) : dim(X) must have a positive length Calls: createHeatmap -> array -> stateColor -> apply Execution halted

I found 'apply(statemean, 1, range)' comes from createGenomeTracks.R. Do you have an idea how can I solve this? Thanks!

Kaili

guanjue commented 6 years ago

First, you need to run this command outside of your output folder: In your case, your output folder is: /data/zusers/fankaili/ideas/test/run_IDEAS_result_chr1/ You need to run it in this folder: cd /data/zusers/fankaili/ideas/test/ time Rscript bin/get_heatmap.R run_IDEAS_result_chr1/run_IDEAS.para0 FALSE bin/createGenomeTracks.R

If this is not the issue, I wonder if you can check what is in the 'run_IDEAS.para0' file? If it is empty?

KailiBio commented 6 years ago

This is another issue. For the previous one, I gave it the bed file with at least 2 chrs, and it runs well in the first '[IDEAS Analysis]' step. It exhaulted at the '[ Create Custom Tracks ]' step this time, and I still trying to find what's statemean in createGenomeTracks.R

guanjue commented 6 years ago

Does the .inv file recreated?

KailiBio commented 6 years ago

Yeah, I do have .inv file.

Below is all the file I got when they exhaulted. screen shot 2018-07-20 at 4 35 28 pm

And this is the log file: screen shot 2018-07-20 at 4 36 26 pm

KailiBio commented 6 years ago

The files listed are in the result folder. I have .inv file in the work directory.

guanjue commented 6 years ago

Can you check if the .inv file to see if the number of bins in the chromosome (third column) matches with the number of bins in each chromosome in your bed file?

KailiBio commented 6 years ago

The .inv file is correct.

guanjue commented 6 years ago

I wonder if you can send me the para0, so I can check the bug tonight? my email address: gzx103@psu.edu

KailiBio commented 6 years ago

Already send it. Many Thanks!