Displaying 1 result from an estimated 1 matches for "curr_assoc".
2009 Feb 02
2
Defining plot colors based on a variable
...5602540 139 0.6182 6.66 6.289e-05
0.09283 0.9262 TCTX 4560767 0.004973939
etc
And this is the code that I have:
assoc_data <- read.table("master.out", header =TRUE)
par(fig=c(0, 10, 0, 10 )/10, mar=c(10,8,2,8),xpd=NA, cex.axis=2)
attach(assoc_data)
curr_assoc <- assoc_data[CHR == 1 & BP > 500000 & BP < 1000000, ] #these
criteria change based on input from another file
#count the number of transcripts
transcripts <- length(unique(curr_assoc$FILE))
#generate that number of unique ³FILE² entries in my subset of data
my_colors <- ra...