Displaying 1 result from an estimated 1 matches for "yachse".
Did you mean:
sachse
2002 Mar 18
2
persp(): add second plane (second, long question)
...package...
What I can do is to get the x,y and z-range for the xlim,ylim and zlim so
that the two planes would fit the box (x and y values are the same for both
(1x6vectors)). And I'm able to store the z-values in two different matrices
(6x6).
What I did then is
> xachse<-c(1:6)
> yachse<-xachse
>
> par(bg = "white")
>
> trans3d <- function(x,y,z,pmat) {
+ tmat <- t((cbind(x,y,z,1)%*% pmat))
+ list(x=tmat[1,]/tmat[4,],y=tmat[2,]/tmat[4,])
+ }
> pmat <- persp(xachse, yachse, betaon1.median.plot,zlim=c(2.75,3.30), theta
= 60, phi = 20,
+...