rforscience

• rforscienceのサイトは、海洋データなどの扱いを中心に、微分方程式の扱い、球面プロット・３次元プロットなどに関していろいろ提供している
• ３次元プロットにはrglパッケージがあるが、少なくとも、ある面では、rforscienceが提供しているplot3Dパッケージは便利そうにみえる
• また、plot3Dパッケージは静止表示だが、それをrgl的に動かすパッケージplot3Drglもあって、これもよい

install.packages(c("plot3D","plot3Drgl"))
library(plot3D)
library(plot3Drgl)
X       <- seq(0, pi, length.out = 50)
 Y       <- seq(0, 2*pi, length.out = 50)
 M       <- mesh(X, Y)
 phi     <- M$x
 theta   <- M$y

# x, y and z grids
 r <- sin(4*phi)^3 + cos(2*phi)^3 + sin(6*theta)^2 + cos(6*theta)^4
 x <- r * sin(phi) * cos(theta)
 y <- r * cos(phi)
 z <- r * sin(phi) * sin(theta)

# full colored image
 surf3D(x, y, z, colvar = y, colkey = FALSE, shade = 0.5,
        box = FALSE, theta = 60)

# rgl化
surf3Drgl(x,y,z,colvar=y,colkey=FALSE,shade=0.5,box=FALSE)