r - R中的3d圆锥和圆柱

Question

你如何在 R 中绘制 3d 圆锥和圆柱？我认为您正在使用线框，但我在任何地方都找不到任何示例。我正在尝试做 Matlab 和 Mathematica 通常用来做的事情。如果绘制 3d 对象可以像那些程序一样简单，我真的很喜欢 R。

Answer 1

这就是 Jaakko Nevalainen 和 Deepayan Sarkar 使用 lattice 的方式来自： https ://stat.ethz.ch/pipermail/r-help/2009-April/196715.html

library(lattice)

A<-matrix(ncol=2, nrow=64)

for(i in 0:63) {  A[i+1,1]<-sin(i/10)  
A[i+1,2]<-cos(i/10) }

Sigma<-matrix(c(0.5,0.1,0.1,0.25),byrow=TRUE,nrow=2) 
G<-eigen(Sigma)

E1<-t(G$vector%*%t(A)) 
E2<-t(diag(sqrt(G$values))%*%t(E1)) 
mu<-c(0.1,0.2) 
E3<-sweep(E2,2,-mu)

a<-sqrt(max(rowSums(sweep(E3,2,mu)**2))) 
b<-sqrt(min(rowSums(sweep(E3,2,mu)**2)))

astar<-as.numeric(a+abs(mu[1])) 
bstar<-as.numeric(b+abs(mu[2]))

xstar<-seq(-astar,astar,len=50) 
ystar<-seq(-bstar,bstar,len=50)

g<-expand.grid(x=xstar,y=ystar)

p1<-2*g$x*mu[1]/a**2+2*g$y*mu[2]/b**2 
p2<-(g$x**2/a**2+g$y**2/b**2) 
p3<-mu[1]**2/a**2+mu[2]**2/b**2-1

q<-(p1+sqrt(p1**2-4*p2*p3))/(2*p2) 
z<-sqrt(1-(q*g$x)**2-(q*g$y)**2) 
zstar<-(z/q) 
ind0<-!(q<1) 
g$z<-zstar 
sc<-matrix(c(rep(c(-1,-1,-1),sum(ind0))),nrow=sum(ind0),byrow=TRUE) 
gstar<-rbind(g[ind0,],sc*g[ind0,])

group<-c(rep(1,nrow(gstar)/2),rep(2,nrow(gstar)/2)) 
gstar$group<-group
gstar <- g 
gstar$z[!ind0] <- NA 
gstar$z2 <- -gstar$z 
wireframe(z + z2 ~ x * y, gstar, colorkey=TRUE,drape=TRUE, scales=list(arrows=FALSE))