par5 <- '' par4 <- '' par3 <- '' par2 <- 'TRUE' par1 <- 'grey' bitmap(file='test1.png') par2<- as.logical(par2) # notches par3<-as.numeric(par3) # percentage trim if(par3>45){par3<-45;warning('trim limited to 45%')} if(par3<0){par3<-0;warning('negative trim makes no sense. Trim is zero.')} par4 <- as.numeric(par4) #factor column par5 <- as.numeric(par5) # response column x <- t(x) x1<-as.numeric(x[,par5]) # response f1<-as.character(x[,par4]) # factor x2<-x1[f1=='no'] f2 <- f1[f1=='no'] lotrm<-as.integer(length(x2)*par3/100) hitrm<-as.integer(length(x2)*(100-par3)/100) srt<-order(x2,f2) trmx1<-x2[srt[lotrm:hitrm]] trmf1<-f2[srt[lotrm:hitrm]] x3<-x1[f1=='yes'] f3 <- f1[f1=='yes'] lotrm<-as.integer(length(x3)*par3/100) hitrm<-as.integer(length(x3)*(100-par3)/100) srt<-order(x3,f3) trmx2<-x3[srt[lotrm:hitrm]] trmf2<-f3[srt[lotrm:hitrm]] xtrm<-c(trmx1,trmx2) ftrm<-c(trmf1,trmf2) xtrm[1:6] ftrm[1:6] r<-boxplot(xtrm~as.factor(ftrm), col=par1, notch=par2, main=main, xlab=xlab, ylab=ylab) dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,hyperlink('overview.htm','Boxplot statistics','Boxplot overview'),6,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Placement',1,TRUE) a<-table.element(a,hyperlink('lower_whisker.htm','lower whisker','definition of lower whisker'),1,TRUE) a<-table.element(a,hyperlink('lower_hinge.htm','lower hinge','definition of lower hinge'),1,TRUE) a<-table.element(a,hyperlink('central_tendency.htm','median','definitions about measures of central tendency'),1,TRUE) a<-table.element(a,hyperlink('upper_hinge.htm','upper hinge','definition of upper hinge'),1,TRUE) a<-table.element(a,hyperlink('upper_whisker.htm','upper whisker','definition of upper whisker'),1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'no',1,TRUE) for (j in 1:5) { a<-table.element(a,r$stats[j,1]) } a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'yes',1,TRUE) for (j in 1:5) { a<-table.element(a,r$stats[j,2]) } a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab') tr.mns<-tapply(x1,f1,mean, trim=par3/100) a<-table.start() a<-table.row.start(a) a<-table.element(a,hyperlink('trimmed_mean.htm','Trimmed Mean Equation','Trimmed Mean'),2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'No Placement') a<-table.element(a,'Placement') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,tr.mns[1]) a<-table.element(a,tr.mns[2]) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable1.tab')
|