library(GenomicRanges)
anchor.one = GRanges(c("chr1", "chr1", "chr1", "chr1"), IRanges(c(10, 20, 30, 20), width=5))
anchor.two = GRanges(c("chr1", "chr1", "chr1", "chr2"), IRanges(c(100, 200, 300, 50), width=5))
interaction_counts = sample(1:10, 4)
test <- GenomicInteractions(anchor.one, anchor.two, experiment_name="test",
description="this is a test", counts=interaction_counts)
name(test) <- "Mouse test"
name(test)
description(test) <- "This is a test using the mouse genome"
description(test)
interactionCounts(test) <- c(2,3,8,5)
interactionCounts(test)
Run the code above in your browser using DataLab