search for: chr9

I have a dataframe in the general format: chr1 0.5 chr1 0 chr1 0.75 chr2 0 chr2 0 chr3 1 chr3 1 chr3 0.5 chr7 0.75 chr9 1 chr9 1 chr22 0.5 chr22 0.5 where the first column is the chromosome location and the second column is some value. What I'd like to do is have a histogram created for each chr location (i.e. a separate histogram for chr1, chr2, chr3, chr7, chr9, and chr22). I am just having a hard time gettin...

...just not have any y-value plotted as a point in that window. I'm attaching what I've been able to plot with ggplot - sorry for the file size. I'd like to be able to plot just the points shaded in blue, for example. Any help is appreciated. Kelly V. Code for the ggplot: > ggplot(chr9, aes(x=Window, y=ExplantInv, colour=ExplantInv))+ + geom_point() + + scale_colour_gradientn(colours=rainbow(4)) Example Data: > head(chr9, n=50) Window Explant ExplantInv 1 168030 0.036590781 0.9634092 2 168031 1.000000000 0.0000000 3 168032 1.000000000 0.0000000 4 168033 1.000...

...ot;chr11", "chr12", "chr13", "chr14", "chr15", "chr16", "chr17", "chr18", "chr19", "chr2", "chr3", "chr4", "chr5", "chr6", "chr7", "chr8", "chr9", "chrX"), class = "factor"), start = c(5291000L, 10988025L, 11767950L, 11840900L, 12267450L, 12276675L), end = c(5291926L, 10988526L, 11768676L, 11841851L, 12268076L, 12277051L), peak = c(8L, 7L, 8L, 8L, 12L, 7L)), .Names = c("chrom", "start", "end...

...also trying to make it look "better" using ggplot2. I have data.frame that contains chrom length 1 chr1 249250621 2 chr2 243199373 3 chr3 198022430 4 chr4 191154276 5 chr5 180915260 6 chr6 171115067 7 chr7 159138663 8 chrX 155270560 9 chr8 146364022 10 chr9 141213431 11 chr10 135534747 12 chr11 135006516 13 chr12 133851895 14 chr13 115169878 15 chr14 107349540 16 chr15 102531392 17 chr16 90354753 18 chr17 81195210 19 chr18 78077248 20 chr20 63025520 21 chrY 59373566 22 chr19 59128983 23 chr22 51304566 24 chr21 48129895 I want...

...DDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD") # 60bp C free sequence query.minus <- reverseComplement(query.plus) chrList <- c("chr1", "chr2", "chr3", "chr4", "chr5", "chr6", "chr7", "chr8", "chr9", "chr10", "chr11", "chr12", "chr13", "chr14", "chr15", "chr16", "chr17", "chr18", "chr19", "chr20", "chr21", "chr22", "chrX", "chrY") #access/...

Hello, I have an input file: http://r.789695.n4.nabble.com/file/n3700031/testOut.txt testOut.txt where col 1 is chromosome, column2 is start of region, column 3 is end of region, column 4 and 5 is base position, column 6 is total reads, column 7 is methylation data, and column 8 is the strand. I would like a summary output file such as: