做个参考吧
- id=seq(-5,5,len=100)
- z <- c(dnorm(id),dnorm(id,0.5,1),dnorm(id,-0.5,1),dnorm(id,0,0.5))
- x <- rep(id, 4)
- y <- rep(1:4, rep(100,4))
- sd=scatterplot3d(x=x[1:100], y=y[1:100], z=z[1:100],type = 'l', highlight.3d=TRUE,
- xlim = c(-5, 5), ylim = c(0, 5),zlim=c(min(z),max(z)),box = FALSE,
- xlab = 'x',ylab = 'y',zlab = 'z'
- )
- sd$points3d(x[101:200], y[101:200], z[101:200],type = 'l',col=2)
- sd$points3d(x[201:300], y[201:300], z[201:300],type = 'l',col=3)
- sd$points3d(x[301:400], y[301:400], z[301:400],type = 'l',col=4)
- legend('topright',c('1','2','3','4'),col=1:4,lty=1,bty = 'n')
复制代码
|