library(GenKern) x <- as.data.frame(t(x)) par1 <- as.numeric(par1) par2 <- as.numeric(par2) par3 <- as.numeric(par3) par4 <- as.numeric(par4) bitmap(file='test1.png') par(mfrow=c(2,2)) op <- KernSur(as.numeric(x[,par3]),as.numeric(x[,par1]), xgridsize=150, ygridsize=150,na.rm=T) image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE, xlab=par7, ylab=par5) contour(op$xords, op$yords, op$zden, add=TRUE) box() lines(stats::lowess(cbind(as.numeric(x[,par3]),as.numeric(x[,par1]))),col='white') op <- KernSur(as.numeric(x[,par4]),as.numeric(x[,par1]), xgridsize=150, ygridsize=150,na.rm=T) image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE, xlab=par8, ylab=par5) contour(op$xords, op$yords, op$zden, add=TRUE) box() lines(stats::lowess(cbind(as.numeric(x[,par4]),as.numeric(x[,par1]))),col='white') op <- KernSur(as.numeric(x[,par3]),as.numeric(x[,par2]), xgridsize=150, ygridsize=150,na.rm=T) image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE, xlab=par7, ylab=par6) contour(op$xords, op$yords, op$zden, add=TRUE) box() lines(stats::lowess(cbind(as.numeric(x[,par3]),as.numeric(x[,par2]))),col='white') op <- KernSur(as.numeric(x[,par4]),as.numeric(x[,par2]), xgridsize=150, ygridsize=150,na.rm=T) image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE, xlab=par8, ylab=par6) contour(op$xords, op$yords, op$zden, add=TRUE) box() lines(stats::lowess(cbind(as.numeric(x[,par4]),as.numeric(x[,par2]))),col='white') dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Kendall Rank Correlations',3,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'',header=TRUE) a<-table.element(a,'tau',header=TRUE) a<-table.element(a,'p-value',header=TRUE) a<-table.row.end(a) r <- cor.test(as.numeric(x[,par3]),as.numeric(x[,par1]),method='kendall') a<-table.row.start(a) dum <- 'tau(' dum <- paste(dum,par7,sep='') dum <- paste(dum,par5,sep=',') dum <- paste(dum,')',sep='') a<-table.element(a,dum,header=TRUE) a<-table.element(a,r$estimate) a<-table.element(a,r$p.value) a<-table.row.end(a) r <- cor.test(as.numeric(x[,par4]),as.numeric(x[,par1]),method='kendall') a<-table.row.start(a) dum <- 'tau(' dum <- paste(dum,par8,sep='') dum <- paste(dum,par5,sep=',') dum <- paste(dum,')',sep='') a<-table.element(a,dum,header=TRUE) a<-table.element(a,r$estimate) a<-table.element(a,r$p.value) a<-table.row.end(a) r <- cor.test(as.numeric(x[,par3]),as.numeric(x[,par2]),method='kendall') a<-table.row.start(a) dum <- 'tau(' dum <- paste(dum,par7,sep='') dum <- paste(dum,par6,sep=',') dum <- paste(dum,')',sep='') a<-table.element(a,dum,header=TRUE) a<-table.element(a,r$estimate) a<-table.element(a,r$p.value) a<-table.row.end(a) r <- cor.test(as.numeric(x[,par4]),as.numeric(x[,par2]),method='kendall') a<-table.row.start(a) dum <- 'tau(' dum <- paste(dum,par8,sep='') dum <- paste(dum,par6,sep=',') dum <- paste(dum,')',sep='') a<-table.element(a,dum,header=TRUE) a<-table.element(a,r$estimate) a<-table.element(a,r$p.value) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
|