[R] Overlapping subject-specific histograms
Frank S.
f_j_rod at hotmail.com
Thu Jan 14 12:47:37 CET 2016
Dear R users,
First of all, excuse me if my doubt is very trivial, but so far I haven't been able to solve it.
My question is this: I have a data frame which contains repeated measurements on 4 subjects coded
as "id", and I want to plot, for each subject, not only the corresponding "counts" variable histogram,
but also overlapping to the right side the corresponding results of "sim" variable (I want to do it in basic
R code, i.e., without any specific R package). I have almost the right code (see the example code below),
but I can not overlap the "sim" variable.
Thanks in advance for suggestions!!
Frank
data <- data.frame(id = rep(c(1,3,4,7), c(9,5,3,3)),
count = c(0, 10, 15, 0, 16, 7, 14, 11, 12, 1, 8, 17, 19, 0, 9, 10, 14, 2, 3, 10),
sims = c(1, 9, 15, 1, 14, 5, 12, 10, 12, 2, 6, 15, 18, 1, 9, 9, 12, 5, 3, 9))
# The actual code I have
# ------------------------------
windows(height = 5, width = 5)
par(mfrow = c(2, 2), oma = c(1, 2, 2, 1), mar=c(3, 2, 1, 1), las = 1)
for(i in 1:length(unique(data$id))){
kat <- factor(data$id, labels = 1:length(unique(data$id)))
plot(data$count[kat == i],
type = "h", col = 1, lwd = 3, xaxt = "n", xlab ="", main = "",
xlim = c(1, max(table(data$id))), ylim = c(0, 20))
axis(1, at = 1:max(table(data$id)))
mtext( bquote(paste("id = ", .(unique(data$id)[i]))), side = 3, cex = 0.9, line = 0.5)
tab <- table( as.matrix( data$id ) )
dist.overlap <- 0.4 # Distance of right overlapping of the "sim" variable
# points( factor(names(tab)) + dist.overlap, data$sim[kat == i] , type="h", col=2, lw =4) ## =======> Line I can not solve
}
[[alternative HTML version deleted]]
More information about the R-help
mailing list