[Bioc-devel] Error with InteractionSet package

[Aaron Lun]

The randomness of the errors is concerning, especially given that none 
of these functions should be calling compiled code yet. Can you produce 
a minimal working example that reproducibly triggers at least one of the 
errors? And post your session information.

-Aaron

On 08/11/16 17:34, Ioannis Vardaxis wrote:
> Hi,
>
> I am trying to combine two GRanges from two anchors of same length in to a Gitnteractions object and I get the following error:
>
>> PairedData=InteractionSet::GInteractions(anchor1=Anchor1,anchor2=Anchor2)
> Error in .new_GInteractions(anchor1 = anchor1, anchor2 = anchor2, regions = regions,  :
>   first and second anchor vectors have different lengths
>> length(Anchor1)
> [1] 33705839
>> length(Anchor2)
> [1] 33705839
>> traceback()
> 5: stop(msg)
> 4: .new_GInteractions(anchor1 = anchor1, anchor2 = anchor2, regions = regions,
>        metadata = metadata, mode = mode)
> 3: .local(anchor1, anchor2, regions, ...)
> 2: InteractionSet::GInteractions(anchor1 = Anchor1, anchor2 = Anchor2)
> 1: InteractionSet::GInteractions(anchor1 = Anchor1, anchor2 = Anchor2)
>
> While they clearly have the same length.
>
> If I then run it again it works fine.
>
> If I then run again one more time I get:
> Error in .new_GInteractions(anchor1 = anchor1, anchor2 = anchor2, regions = regions,  :
>   all anchor indices must be finite integers
> In addition: Warning message:
> In new.pos[o] <- new.pos :
>   number of items to replace is not a multiple of replacement length
>
> I run it again and I get the correct results.
>
> It seems like it is quite random if I get en error or not.
>
> --
> Ioannis Vardaxis
> Stipendiat IMF
> NTNU
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> Bioc-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/bioc-devel
>