[R] Overlapping subject-specific histograms

Frank S.
2016-01-14T11:47:37Z

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]]