I have a function in 3D - let's say De Jong function:
fdejong <- function (x, y) {
return (x^2 + y^2)
}
How can I draw it's plot in 3D? I want to achieve an effect similar to this from wikipedia:
Asked
Active
Viewed 1.2k times
5
-
see the `curve3d()` function from the `emdbook` package. – Ben Bolker Jun 05 '14 at 16:42
-
stackoverflow.com/questions/11875941/3d-equivalent-of-the-curve-function-in-r/ – Ben Bolker Jun 05 '14 at 16:43
2 Answers
9
try this :
fdejong <- function (x, y) {
return (x^2 + y^2)
}
x <- seq(-10, 10, length= 30)
y <- x
z <- outer(x, y, fdejong)
z[is.na(z)] <- 1
op <- par(bg = "white")
persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "lightblue")
MrFlick
- 195,160
- 17
- 277
- 295
Vincent Guyader
- 2,927
- 1
- 26
- 43
-
It's almost what I want, but is it any way to color it like that plot in wikipedia (color depend on the function value)? – SathOkh Jun 05 '14 at 16:56
-
8
You can also use the Lattice wireframe function. (Using @user1020027's data)
fdejong <- function (x, y) {
return (x^2 + y^2)
}
x <- seq(-10, 10, length= 30)
y <- x
z <- outer(x, y, fdejong)
z[is.na(z)] <- 1
require(lattice)
wireframe(z, drape=T, col.regions=rainbow(100))
MrFlick
- 195,160
- 17
- 277
- 295