[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