Aimin, please provide self-contained code. This is not "data" but a printed table. Next time, provide real or simulated data in a way that it can be copied directly into R or an R editor. Also, it is quite impossible to tell from your question how you want to plot it. The grouping is not the problem, but what do you mean by "plot Mean with SD(+/-)." How should it look like? Also, some SDs are NAs. What to do about them?
Here is my guess for what I think you might want to do. Pay special attention to the tapply commands that compute mean MEANS and mean SDs for each LACxDIM bin. Also pay special attention to how you have to recode the DIM and LAC variables (here recoded into DIM2 and LAC2) to plot the MEAN and SDs with the right DIM and LAC indicator. #Simulate data DIM=rep(1:4,each=10) LAC=rep(1:5,8) MEAN=rnorm(40,45,10) SD=rgamma(40,15,2) #Summary for DIMxLAC bins tapply(MEAN,list(DIM,LAC),mean) #the same in one vector #pay attention how they are ordered c(tapply(MEAN,list(DIM,LAC),mean)) MEAN2=c(tapply(MEAN,list(DIM,LAC),mean)) SD2=c(tapply(SD,list(DIM,LAC),mean)) INT.plus=MEAN2+SD2 INT.minus=MEAN2-SD2 #recode LAC2 and DIM2 accordingly #to match the order in MEAN2 and SD2 LAC2=rep(1:5,each=4) DIM2=rep(1:4,5) #plot #may require to install package Hmisc library(Hmisc) Dotplot(LAC2~Cbind(MEAN2,INT.minus,INT.plus)|DIM2,method="bars") hth Daniel ------------------------- cuncta stricte discussurus ------------------------- -----Ursprüngliche Nachricht----- Von: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Im Auftrag von jlfmssm Gesendet: Tuesday, June 23, 2009 6:07 PM An: r-help@r-project.org Betreff: [R] ask for help xyplot I want to make such plot using the following data set plot Mean with Sd(+/-) by Dim group by Lac here is data Does anyone know how to plot using xyploy Thanks, Aimin Mean Sd Var Min Max Dim Lac 704 44.00000 NA NA 44.0 44.0 -30 3 703 45.92000 9.5484030 91.1720000 60.0 34.4 -30 4 702 57.40000 NA NA 57.4 57.4 -30 5 707 49.00000 NA NA 49.0 49.0 -29 3 706 43.12000 16.5762481 274.7720000 61.0 16.8 -29 4 705 60.40000 NA NA 60.4 60.4 -29 5 710 32.65000 21.2987480 453.6366667 52.0 8.4 -28 3 709 36.30000 13.2717746 176.1400000 51.2 22.2 -28 4 708 39.60000 NA NA 39.6 39.6 -28 6 713 43.84000 16.3226223 266.4280000 57.4 15.8 -27 3 712 39.46667 10.7930842 116.4906667 50.2 24.8 -27 4 711 23.40000 NA NA 23.4 23.4 -27 6 717 36.86667 16.5373113 273.4826667 52.0 11.0 -26 3 716 51.77143 21.5293377 463.5123810 93.0 20.2 -26 4 715 56.20000 NA NA 56.2 56.2 -26 5 714 38.00000 NA NA 38.0 38.0 -26 6 628 33.33333 22.0927741 488.0906667 60.6 0.6 -25 3 627 46.20000 11.9002521 141.6160000 55.0 22.8 -25 4 626 45.00000 16.1220346 259.9200000 56.4 33.6 -25 5 738 36.76000 16.8014285 282.2880000 49.8 9.4 -24 3 737 45.05000 17.6313168 310.8633333 61.2 21.6 -24 4 736 53.70000 2.1213203 4.5000000 55.2 52.2 -24 5 122 38.14286 17.1455395 293.9695238 58.0 8.0 -23 3 121 49.36667 18.7456306 351.3986667 69.2 15.6 -23 4 120 50.25000 9.2885234 86.2766667 58.8 37.8 -23 5 119 80.00000 NA NA 80.0 80.0 -23 6 126 47.37778 5.8072751 33.7244444 55.2 37.8 -22 3 125 59.10000 18.5004865 342.2680000 87.6 40.8 -22 4 124 57.44000 24.3065423 590.8080000 99.0 38.6 -22 5 123 30.40000 NA NA 30.4 30.4 -22 7 129 50.00000 13.8258557 191.1542857 78.0 30.0 -21 3 128 40.92000 20.8502278 434.7320000 62.8 7.0 -21 4 127 49.80000 1.0392305 1.0800000 50.4 48.6 -21 5 721 52.20000 6.2369865 38.9000000 61.2 43.2 -20 3 720 49.36000 6.7740682 45.8880000 56.8 42.2 -20 4 [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org 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. ______________________________________________ R-help@r-project.org 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.
