Skip to content

Biocircos Link track automation using loop #28

@khush876

Description

@khush876

I would like to automate BioCircos plot in R and here is my code. However, I am not able to link the tracks when I am using the full file.

Input file

All_C1.txt

query pos1 sub pos2
AUSP_1 15147 AUSA_1 93903
AUSP_3 15147 AUSA_2 93903
AUSP_9 15204 AUSA_4 93960
AUSP_2 15242 AUSA_3 93998
AUSP_3 15264 AUSA_5 94021
AUSP_4 15287 AUSA_6 94044
AUSP_5 15309 AUSA_7 94066
AUSP_6 15313 AUSA_8 94070
AUSP_7 15319 AUSA_9 94076
AUSP_8 15320 AUSA_1 94077

library(BioCircos)
setwd("/PATH/TO/FILE")
data <- read.table("All_C1.txt")
head(data)

myGenomeAus0 = list("AUSP_1" = 31248581,
"AUSP_2" = 31645086,
"AUSP_3" = 47108296,
"AUSP_4" = 29969644,
"AUSP_5" = 36443716,
"AUSP_6" = 26560673,
"AUSP_7" = 28486927,
"AUSP_8" = 29326496,
"AUSP_9" = 29774650,
"AUSA_1" = 31122193,
"AUSA_2" = 32101233,
"AUSA_3" = 46783639,
"AUSA_4" = 29229108,
"AUSA_5" = 36358360,
"AUSA_6" = 26077416,
"AUSA_7" = 27897604,
"AUSA_8" = 29233132,
"AUSA_9" = 30348868)

links_chromosomes_1 <- c(data)
V1 <- links_chromosomes_1[1]
V1

links_pos_1 <- c(data)
V2 <- links_pos_1[2]
V2

links_chromosomes_2 <- c(data)
V3 <- links_chromosomes_2[3]
V3

links_pos_2 <- c(data)
V4 <- links_pos_2[4]
V4

for (i in data$V1[2]){
print (i)
for (j in data$V2[2]){
print (j)
for (k in data$V3[2]){
print (k)
for (l in data$V4[2]){
print (l)

tracklist = BioCircosLinkTrack('myLinkTrack', i, j, j,
k, l, l, genome = myGenomeAus0,
maxRadius = 1, color = "black", width = "1.0em",
displayAxis = TRUE, axisColor = "black", axisWidth = 0.5,
axisPadding = 0, displayLabel = TRUE, labelColor = "black",
labelSize = "1em")
}
}
}
}

BioCircos(tracklist, genome = myGenomeAus0, genomeFillColor = c("Spectral"),
chrPad = 0.05, displayGenomeBorder = FALSE,
genomeTicksDisplay = FALSE, genomeLabelTextSize = "20pt", genomeLabelDy = 0)

Thanks

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions