Palimpsest icon indicating copy to clipboard operation
Palimpsest copied to clipboard

Error running deconvolution_fit_SV with plotting

Open anoronh4 opened this issue 3 years ago • 5 comments

Using the example script with the example Rdata gives me the following error:

Error in RCircos.Validate.Genomic.Data(genomic.data, "link", genomic.columns) : 
  One or more genomic position in plot data is
outside of chromosome length for chr4.
Calls: deconvolution_fit_SV ... RCircos.Get.Paired.Points.Positions -> RCircos.Validate.Genomic.Data

anoronh4 avatar Jul 28 '22 21:07 anoronh4

I have the same problem! Could you solve it? Thanks!

edceeyuchen avatar Oct 19 '22 11:10 edceeyuchen

For anyone with this problem in the future, deconvolution_fit_SV has an argument gv, or 'genome version'. This is set to "hg38" default, and needs to be manually set for the script to "hg19" instead.

jpmcelderry avatar Jan 28 '23 00:01 jpmcelderry

@jpmcelderry Hello! Sorry for lately reply, I want to know, if I used my SV data, there is no need to change the 'genome version', right? Hope your help~

edceeyuchen avatar Mar 30 '23 04:03 edceeyuchen

@jpmcelderry Hello! Sorry for lately reply, I want to know, if I used my SV data, there is no need to change the 'genome version', right? Hope your help~

I imagine it just depends on your SV data genome version. If your SV data is in hg38 coordinates then this matches the default value of gv. If your data is in hg19 coordinates you will need to set gv="hg19".

jpmcelderry avatar Mar 30 '23 18:03 jpmcelderry

@jpmcelderry Hello! Sorry for lately reply, I want to know, if I used my SV data, there is no need to change the 'genome version', right? Hope your help~

I imagi ne it just depends on your SV data genome version. If your SV data is in hg38 coordinates then this matches the default value of gv. If your data is in hg19 coordinates you will need to set gv="hg19".

Thank you! And totally I understand now~

edceeyuchen avatar Mar 31 '23 02:03 edceeyuchen