Displaying 2 results from an estimated 2 matches for "genedensity".
2012 Jul 09
1
boxplot with "cut"
Dear UseRs,
I'm making box plots from a data set that looks like this:
Chr Start End GeneDensity ReadCount_Explant ReadCount_Callus ReadCount_Regen
1 1 1 10000 107.82 1.243 1.047 1.496
2 1 10001 20000 202.50 0.835 0.869 0.456
3 1 20001 30000 158.80 1.813 1.529 1.131
4 1...
2012 Jul 10
0
Thanks! RE: boxplot with "cut"
...nge(unlist(x)), xaxt = "n", ...)
for(i in seq_len(n)[-1])
boxplot(x[[i]] ~ by, at = 0:len*n1 + i, xaxt = "n", add = TRUE, ...)
axis(1, at = 0:len*n1 + n1/2, labels = names(sp), tick = TRUE) }
cols <- grep("ReadCount", names(GD))
multi.boxplot(GD[, cols], cut(GD$GeneDensity, breaks=10))
If this is it and you don't like those x-axis tick lables, use as.integer(cut(...etc...)).
Hope this helps,
Rui Barradas
Em 09-07-2012 20:51, Vining, Kelly escreveu:
> Dear UseRs,
> I'm making box plots from a data set that looks like this:
>
>
> Chr Star...