par1 = 12greygrey0.010.10,010,010,010,10,10,10,010,010,010,010,112DefaultDefaultDefault486060121212additive120.01 ; par2 = nono0.990.90,990,990,990,990,990,990,990,990,990,990,99grey11111112Triple0.99 ; par3 = 0.010.10,010,010,010,10,10,10,010,10,10,10,1FALSE000001additive0.01 ; par4 = Unknown00000012 ; par5 = 121212121212 ; par6 = White NoiseWhite NoiseWhite NoiseWhite NoiseWhite NoiseWhite Noise ; par7 = 0.950.950.950.950.950.95 ; par8 = 480 ; |
par1 <- as(par1,'numeric') par2 <- as(par2,'numeric') par3 <- as(par3,'numeric') library(Hmisc) myseq <- seq(par1, par2, par3) hd <- hdquantile(x, probs = myseq, se = TRUE, na.rm = FALSE, names = TRUE, weights=FALSE) bitmap(file='test1.png') plot(myseq,hd,col=2,main=main,xlab=xlab,ylab=ylab) grid() dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Harrell-Davis Quantiles',3,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'quantiles',header=TRUE) a<-table.element(a,'value',header=TRUE) a<-table.element(a,'standard error',header=TRUE) a<-table.row.end(a) length(hd) for (i in 1:length(hd)) { a<-table.row.start(a) a<-table.element(a,as(labels(hd)[i],'numeric'),header=TRUE) a<-table.element(a,signif(as.matrix(hd[i])[1,1],6)) a<-table.element(a,signif(as.matrix(attr(hd,'se')[i])[1,1],6)) a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable.tab')
|