hello,
maybe this code can be useful for you.
cleber
-----------------------
trimage <- function(f){
x = y = seq( 1, 0, l=181 )
t1 = length(x)
im = aux = numeric(0)
for( i in seq( 1, t1, by = 2 ) ){
#idx = seq( t1**2, i*t1, by = -t1 ) - ((t1 - i):0)
idx = seq( i*t1, t1**2, by = t1 ) - (i-1)
im = c(im, aux, idx, aux )
aux = c(aux, NA)
}
z = outer(X=x, Y=y, FUN=f)
return( matrix(z[im],nr=t1) )
}
### for chemical mixtures
### restriction: sum( x[i]==1 ) and 0 < x[i] < 1
### naive example
f <- function(x1,x2) { x3=1-x2-x1; -100*x1 + 0*x2 + 100*x3 }
windows(w=4.5, h=4.5, restoreConsole = TRUE )
par(mar=c(0,0,0,0), pty='s', xaxt='n', yaxt='n', bty='n' )
trimat <- trimage( f )
image( trimat )
contour( trimat, add=T)
Em 16/2/2010 11:25, kajo escreveu:
Hi all,
I am working on a filled contour plot which shows a triangular matrix data
set (as shown below). Is there a possibilty to draw a triangular filled
contour in a equilateral triangle (like a ternary plot)?
Thanks in advance
Johannes
http://n4.nabble.com/file/n1557386/Bild3.png
--
O bom senso é a coisa do mundo mais bem distribuída:
Todos pensamos tê-lo em tal medida que até os mais difíceis
de contentar nas outras coisas não costumam desejar mais bom senso do que
aquele que têm.
[René Descartes]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
