Nothing
# Version: 30-11-2012
estPlotPair <- function(x,col,highlight,hlCol,pch,zoom,...){
ifelse(is.matrix(x$probs),Nprobs <- dim(x$probs)[1], Nprobs <- length(x$probs))
ifelse(is.matrix(x$probs),Nval <- dim(x$probs)[2], Nval <- 1)
markThese <- rep(col,Nval)
markThese[highlight] <- hlCol
pairComb <- getPairComb(x$goi,x$order)
oldpar <- par(no.readonly = TRUE) # code line i
on.exit(par(oldpar)) # code line i + 1
par(mfrow=c(Nprobs-1,Nprobs-1),
pty="s",
bty="c",
oma=c(0,0,3,0),
mar=c(0,0,0,0))
run <- 1
for(xRun in 1:(Nprobs-1))
{
# First fill up with empty fields
if(xRun>1)
{
for(i in 1:(xRun-1))
{
plot(-10,-10,xlim=c(0,1),ylim=c(0,1),xaxt="n",yaxt="n",pch=pch,xlab="",ylab="",xaxt="n",yaxt="n")
if(i==1)
{
name1 <- paste("P(", pairComb[run,1],"<",pairComb[run,2],")",sep="")
if(is.matrix(x$probs)){
mtext(paste(rownames(x$probs)[rownames(x$probs)==name1],sep=""), WEST<-2, line=0.3, adj=0.5, cex=0.8, col="black", outer=FALSE)
} else {
mtext(paste(names(x$probs)[names(x$probs)==name1],sep=""), WEST<-2, line=0.3, adj=0.5, cex=0.8, col="black", outer=FALSE)
}
axis(2,at=seq(0,1,0.2),labels=seq(0,1,0.2),tck=0.01,mgp=c(0,-1.5,0))
}
}
}
yValues <- getPairProb(x,pairComb[run,1],pairComb[run,2])
yMin <- ifelse(zoom,min(yValues),0)
yMax <- ifelse(zoom,max(yValues),1)
for(yRun in (xRun+1):Nprobs)
{
if(xRun==1 & yRun == 2) par(oma=c(0,3,3,0), mar=c(0,0,0,0))
#if(xRun==(Nprobs-1) & yRun == Nprobs) par(bty="o")
ifelse(yRun==(Nprobs), par(bty="o"),par(bty="c"))
xValues <- getPairProb(x,pairComb[run,3],pairComb[run,4])
xMin <- ifelse(zoom,min(xValues),0)
xMax <- ifelse(zoom,max(xValues),1)
plot(xValues,yValues,xlim=c(max(xMin-0.05,0),min(xMax+0.05,1)),ylim=c(max(yMin-0.05,0),min(yMax+0.05,1)),xaxt="n",yaxt="n",col=markThese,pch=pch,xlab="",ylab="",xaxt="n",yaxt="n")
lines(c(0.5,0.5),c(-1,2),lty="dotted")
lines(c(-1,2),c(0.5,0.5),lty="dotted")
if(!is.null(highlight)) {
points(xValues[highlight],yValues[highlight],xaxt="n",yaxt="n",col=hlCol,pch=pch,xlab="",ylab="",xaxt="n",yaxt="n")
}
name1 <- paste("P(", pairComb[run,1],"<",pairComb[run,2],")",sep="")
name2 <- paste("P(", pairComb[run,3],"<",pairComb[run,4],")",sep="")
if(xRun==1){
if(is.matrix(x$probs)){
mtext(paste(rownames(x$probs)[rownames(x$probs)==name2],sep=""), NORTH<-3, line=0.3, adj=0.5, cex=0.8, col="black", outer=FALSE)
} else {
mtext(paste(names(x$probs)[names(x$probs)==name2],sep=""), NORTH<-3, line=0.3, adj=0.5, cex=0.8, col="black", outer=FALSE)
}
}
if(xRun==1 & yRun == 2){
if(is.matrix(x$probs)){
mtext(paste(rownames(x$probs)[rownames(x$probs)==name1],sep=""), WEST<-2, line=0.3, adj=0.5, cex=0.8, col="black", outer=FALSE)
} else {
mtext(paste(names(x$probs)[names(x$probs)==name1],sep=""), WEST<-2, line=0.3, adj=0.5, cex=0.8, col="black", outer=FALSE)
}
}
if(xRun==1 & yRun == 2) axis(2,at=seq(0,1,0.2),labels=seq(0,1,0.2),tck=0.01,mgp=c(0,-1.5,0))
axis(1,at=seq(0,1,0.2),labels=seq(0,1,0.2))
run <- run + 1
}
}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.