search for: colsumz

Hi This question is far less simple than the title suggests, please read carefully, thanks. I have 2 sets of data, both read into R >data1<-read.table ("1.txt", header=T, sep="\t") >data2<-read.table ("2.txt", header=T, sep="\t") >data1 Taxon stage1 stage2 stage3 stage4 T1 0 0 1 1 T2 0

Is there a way to create a 'bubble plot' in R? For example, if we define the following data frame containing the level of y observed for 5 patients at three time points: time<-c(rep('time 1',5),rep('time 2',5),rep('time 3',5))

...xt="n", # no y axis lables bty="n", # no box around the plot xaxs = "i", yaxs = "i", xlim=xlim, ylim=ylim, ... ) ny <- nlevels(ztab$y) nx <- nlevels(ztab$x) sumz <- sum(ztab$z, na.rm=TRUE) colsumz <- sapply(split( ztab$z, ztab$y), sum, na.rm=TRUE) # works rowsumz <- sapply(split( ztab$z, ztab$x), sum, na.rm=TRUE) # broken if(show.margins) { ## column totals text( x=(1:nx) + nlabels.y*rowmar + 0.25 -1, y=0.25, labels=format(c(sumz,...