On Fri, 18 Feb 2005 10:50:37 -0500 Campo Elías PARDO wrote: > Thank for your help > > I obtained profiles and I found mosaicplot as an interesting > alternative. > > I don't like my solution about legend in profiles graphics: I inserted > empty extra columns in order to avoid tue superimposed of legend. > > #Data > N <- matrix(0,3,6) > N[1,] <- c(7,7,5,0,4,4) > N[2,] <- c(0,0,0,5,5,5) > N[3,] <- c(4,4,0,0,3,0) > rownames(N) <- c("Ple1","Ple2","Ple3") > colnames(N) <- c("A","B","C","D","E","F") > N > # Row profiles > PF <- N/rowSums(N)
This could also be done by prop.table(N, 1) > # Columns profiles > PC <- t(N)/rowSums(t(N)) and this by prop.table(N, 2) > # Graphics > par(mfrow=c(2,2)) > barplot(cbind(cbind(t(PF),0),0),legend=colnames(F),density=100,cex.na Note, that you could also do cbind(t(PF), 0, 0) Instead of adding 0s, you could also set xlim. > mes=0.8) > barplot(cbind(cbind(cbind(t(PC),0),0),0),legend=rownames(F),density=1 > 00,cex.names=0.8) > mosaicplot(N,color=TRUE) Note that the mosaicplot is very similar to the stacked barplot (in fact it is a generalization of stacked barplots). R> mosaicplot(N, color = TRUE) gives you the distribution of A-F conditional on Ple1-3, but it also visualizes the marginal distribution of Ple1-3. If you should want to ignore the latter, you could do R> mosaicplot(prop.table(N, 1), color = TRUE) and the latter is essentially equivalent to your barplot R> barplot(t(PF), legend = TRUE, xlim = c(0.2, 4)) And, vice versa, if you want to visualize the conditional distribution of Ple1-3 given A-F, you can do. R> mosaicplot(t(N), color = TRUE) R> mosaicplot(t(prop.table(N, 2)), color = TRUE) hth, Z > library(ade4) > table.cont(t(N),csize=4.5,col.labels=rownames(N)) > > > > Campo Elías PARDO > [EMAIL PROTECTED] > [EMAIL PROTECTED] > ______________________________________________ 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