Re: [R] barchart producing incorrect number of barcharts when	columns renamed On 25 Apr 2007, at 19:50, Deepayan Sarkar wrote: > You seem to have missed this footer that appears in every r-help > message. Your code is not reproducible, and not minimal by a long, > long, shot. > Indeed you are correct. My apologies. This more minimal code should demonstrate the problem. I've removed   all references to external files and verified that I have the same   problem without them. When I change the column names I get a 4x4 array of barcharts,   whereas when I do not change the column names I correctly get the 5x5   array. The question is what causes this discrepancy? library('lattice') rm(list=ls()) textFontSize=6; clipping="(NOT Clipped)"; roi.errs=data.frame(    cbind(          ##Group          c("ctrl", "ctrl", "short", "short", "long", "long"),          ##Subject          c("66101_2", "66112_2", "66119", "66121", "66026_3",   "66030_2")),    cbind(          ##Sub-brick          c(0, 0, 0, 0, 0, 0),          ##Mean_1          c(0.546532, 6.057756, 0.241484, 8.861095, 6.712401, 6.714763),          ##Mean_2          c(0.541001, 1.821168, -0.071227, 7.387406, 4.164371, 6.196456),          ##Mean_3          c(-0.211658, 1.053770, 4.772759, 9.780944, 7.530651, 3.357583),          ##Mean_4          c(0.592144, -2.386457, 4.998830, 12.428849, 14.037164,   4.049794),          ##Mean_5          c(1.449212, 2.445002, -2.718933, 8.900946, 7.919560, 8.350304),          ##Mean_6          c(1.919896, 2.466083, -3.906275, 0.822132, 5.910088, 4.678836),          ##Mean_7          c(0.665386, -3.790489, 1.191782, 3.613919, 3.808240, 5.036492),          ##Mean_8          c(0.865090, 1.741351, 2.012702, 10.093662, 0.829254, 3.880265),          ##Mean_9          c(0.342638, -1.200232, 9.458908, 7.756559, 6.308535, 1.497788),          ##Mean_10          c(1.522509, 1.780352, 4.471669, -1.122842, 4.923762, 4.025056),          ##Mean_11          c(0.615833, 1.363775, 0.215138, 5.988049, 7.545206, 2.117970),          ##Mean_12          c(-0.932574, -1.140839, -5.778118, 5.749871, 8.312596,   4.937120),          ##Mean_13          c(-1.403353, 9.795787, 1.129139, 2.383282, 4.998841, 3.721432),          ##Mean_14          c(0.534640, 4.651436, -0.163107, 3.714310, 6.594728, 3.636995),          ##Mean_15          c(0.951416, 2.384751, -2.804867, 4.662075, 4.375787, 5.767994),          ##Mean_16          c(0.333624, -4.837685, -2.951958, -3.335181, 2.915522,   2.617556),          ##Mean_17          c(2.184298, 7.578843, -2.862990, 6.923186, -2.569342,   0.913788),          ##Mean_18          c(0.382605, 6.441820, 0.579722, -2.424237, 6.884646, 3.127198),          ##Mean_19          c(-0.737142, 3.306980, 2.417347, -7.521070, 6.570706,   1.556814),          ##Mean_20          c(0.545103, 8.721134, -1.766060, 21.433794, 5.829050,   4.031095))); colnames(roi.errs)=c("Group","Subject","Sub-  brick","Mean_1","Mean_2","Mean_3","Mean_4","Mean_5","Mean_6","Mean_7", "Mean_8","Mean_9","Mean_10","Mean_11","Mean_12","Mean_13","Mean_14","Mea  n_15","Mean_16","Mean_17","Mean_18", "Mean_19","Mean_20"); roi.errs.names=names(roi.errs); ## ctrl roi.errs.ctrl<-roi.errs[roi.errs[,"Group"]=="ctrl", 4:length(roi.errs)] roi.errs.ctrl.subjects=roi.errs[roi.errs[,"Group"]=="ctrl", 2] ## short roi.errs.short<-roi.errs[roi.errs[,"Group"]=="short", 4:length  (roi.errs)] roi.errs.short.subjects=roi.errs[roi.errs[,"Group"]=="short", 2] ## long roi.errs.long<-roi.errs[roi.errs[,"Group"]=="long", 4:length(roi.errs)] roi.errs.long.subjects=roi.errs[roi.errs[,"Group"]=="long", 2] roi.errs.ctrl.matrix=as.matrix(roi.errs.ctrl) roi.errs.short.matrix=as.matrix(roi.errs.short) roi.errs.long.matrix=as.matrix(roi.errs.long) ##################################### ### WITH column name change ##################################### roi.errs.ctrl.means=colMeans(roi.errs.ctrl.matrix) roi.errs.short.means=colMeans(roi.errs.short.matrix) roi.errs.long.means=colMeans(roi.errs.long.matrix) yylim=c(floor(min(roi.errs[, 4:length(roi.errs)])), ceiling(max  (roi.errs[, 4:length(roi.errs)]))) errs.Means=rbind(roi.errs.ctrl.means, roi.errs.short.means,   roi.errs.long.means) rownames(errs.Means)=c('control', 'short', 'long') cat("errs.Means dimensions before col name change ", dim(errs.Means),   "\\n"); colnames(errs.Means) = list("R Cing Gy", "R Smar Gy", "L Inf Fron Gy",            "R Thalamus", "R Mid Fron Gy", "L Mid Temp Gy", "L Inf Par   Lbl",            "R Mid Fron Gy", "L Cing Gy", "R Thalamus", "R Mid Fron   Gy", "R Lent Nuc",            "R Inf Fron Gy", "R Mid Temp Gy", "R Sup Par Lbl", "R Cere   Ton",            "L Prec", "R Cere Ton", "L Pyramis", "R Prec Gy"); cat("errs.Means dimensions after col name change ", dim(errs.Means),   "\\n"); get(getOption("device"))() fontsize=trellis.par.get("fontsize"); fontsize$text=textFontSize; trellis.par.set("fontsize", fontsize); print(barchart(errs.Means, groups=rownames(errs.Means), xlab='Mean   Intensity',                 main=paste("Mean Cluster Intensity for Errors",   clipping),                 ylab='Group', col=rainbow(length(rownames  (errs.Means))), border=rainbow(length(rownames(errs.Means))))) ##################################### ### WITHOUT column name change ##################################### errs.Means=rbind(roi.errs.ctrl.means, roi.errs.short.means,   roi.errs.long.means) rownames(errs.Means)=c('control', 'short', 'long') get(getOption("device"))() fontsize=trellis.par.get("fontsize"); fontsize$text=textFontSize; trellis.par.set("fontsize", fontsize); print(barchart(errs.Means, groups=rownames(errs.Means), xlab='Mean   Intensity',                 main=paste("Mean Cluster Intensity for Errors",   clipping),                 ylab='Group', col=rainbow(length(rownames  (errs.Means))), border=rainbow(length(rownames(errs.Means))))) Thanks, Regards, -- Dr Colm G. Connolly School of Psychology and Institute of Neuroscience The Lloyd Building University of Dublin Trinity College, Dublin 2, Ãire Tel: +353-1-896-8475 Fax: +353-1-671-3183 ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.