par5 <- '2'
par4 <- '1'
par3 <- '20'
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')
|
