myprojects-hg-reborn

Fork

rm(list=ls())

#let us read all the data we need

slice <- 20

edge<-read.table('edge_list_complete.dat', header=FALSE)
edge <- as.matrix(edge)

sel <- which(edge[,2]!=edge[,3])

edge_red <- edge[sel, ] 


edge[,1 ] <- (edge[,1]-min(edge[,1]))/slice

edge_rev<-read.table('edge_list_complete_rev.dat', header=FALSE)
edge_rev <- as.matrix(edge_rev)

sel_rev <- which(edge_rev[,2]!=edge_rev[,3])

edge_red_rev <- edge_rev[sel_rev, ] 



edge_rev[,1 ] <- (edge_rev[,1]-min(edge_rev[,1]))/slice



r1_rev <- scan("reg_1_rev.dat", what="")
r2_rev <- scan("reg_2_rev.dat", what="")
r3_rev <- scan("reg_3_rev.dat", what="")

r1 <- scan("reg_1_.dat", what="")
r2 <- scan("reg_2_.dat", what="")
r3 <- scan("reg_3_.dat", what="")

target <- "7359700484475972"
target2 <- 73597004

sel <- which(r3 == target)

sel_rev <- which(r3_rev==target)

time <- as.numeric(r1[sel])
time_rev <- as.numeric(r1_rev[sel_rev])

sel_edge <- c(which((edge[,2]==target2)), which((edge[,3]==target2)))

sel_edge_rev <- c(which((edge_rev[,2]==target2)), which((edge_rev[,3]==target2)))


new_target <- 78118966

sel <- c(which(edge_red[,2]==new_target), which(edge_red[,3]==new_target))

target_appearance <- edge_red[sel,]


sel_rev <- c(which(edge_red_rev[,2]==new_target), which(edge_red_rev[,3]==new_target))

target_appearance_rev <- edge_red_rev[sel_rev,]



print("So far so good")