Displaying 1 result from an estimated 1 matches for "maxabsr".
Did you mean:
imaxabs
2011 Jun 16
1
lines(..., lwd=3) inaccuracy
...r(mar=c(1,2,0,0), xpd=NA)
set.seed(1)
r <- matrix(runif(16), nrow=4, dimnames=list(NULL, paste('x', 1:4, sep='')))
p <- nrow(r)
v <- colnames(r)
plot(c(-.35,p+.5),c(.5,p+.25), type='n', axes=FALSE,
xlab='',ylab='')
text(rep(.5,p), 1:p, v, adj=1)
maxabsr <- max(abs(r[row(r) != col(r)]))
for(i in 1:p) {
for(j in 1:p) {
if(i >= j) next
lines(c(i,i),c(j,j+r[i,j]/maxabsr/2), lwd=3)
lines(c(i-.2,i+.2),c(j,j), col=gray(.7))
}
text(i, i, v[i], srt=-45, adj=0)
}
-----
Frank Harrell
Department of Biostatistics, Vanderbilt Universit...