R中如何画出这种三维网格图

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请高手帮忙，这种图R肯定可以画，但是该用什么包什么函数呢？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏5 回帖

关键词：如何

1. require(grDevices)
   
2. x <- seq(-10, 10, length= 30)
   
3. y <- x
   
4. f <- function(x, y) { r <- sqrt(x^2+y^2); 10 * sin(r)/r }
   
5. z <- outer(x, y, f)
   
6. z[is.na(z)] <- 1
   
7. op <- par(bg = "white")
   
8. 
   
9. persp(x, y, z, theta = 30, phi = 30,
   
10.       expand = 0.5, col = drapecol(z))
    
11. persp(x, y, z, theta = 45, phi = 20,
    
12.       expand = 0.5, col = drapecol(z),
    
13.       r=180,
    
14.       ltheta = 120,
    
15.       shade = 0.75,
    
16.       ticktype = "detailed",
    
17.       xlab = "X", ylab = "Y", zlab = "Sinc( r )" ,
    
18.       #border=30
    
19. )

复制代码

我也新手，网格线和投影没弄出来，哪位高手补充？？

ryoeng 发表于 2015-2-4 15:18
rlg程序包
http://cran.r-project.org/web/packages/rgl/index.html

谢谢指点！！！！

ryoeng 发表于 2015-2-4 15:18
rlg程序包
http://cran.r-project.org/web/packages/rgl/index.html

rlg? 打错了吧，不是rgl 吗

同求解决方案，楼主解决了么

niu9146 发表于 2015-2-4 13:34
我也新手，网格线和投影没弄出来，哪位高手补充？？

大神，我这有两组时间序列数据（X,Y），只能知道他们各自的经验分布函数（分布函数未知），怎么画出类似上面的图，考察他们的尾部相依结构

require(grDevices)
> x <- seq(-10, 10, length= 30)
> y <- x
> f <- function(x, y) { r <- sqrt(x^2+y^2); 10 * sin(r)/r }
> z <- outer(x, y, f)
> z[is.na(z)] <- 1
> op <- par(bg = "white")
>
> persp(x, y, z, theta = 30, phi = 30,
+       expand = 0.5, col = drapecol(z))
Error in persp.default(x, y, z, theta = 30, phi = 30, expand = 0.5, col = drapecol(z)) :
  没有"drapecol"这个函数


请问这是为什么呢

╰﹀ヤ埖瓣雨 发表于 2016-3-4 16:48
require(grDevices)
> x  y  f  persp(x, y, z, theta = 30, phi = 30,
+       expand = 0.5, col = d ...

col = drapecol(z)改成col = rainbow(1000)试试

niu9146 发表于 2015-2-4 13:34
我也新手，网格线和投影没弄出来，哪位高手补充？？

您这个是z=x^2+y^2吧？也就是说z轴上的数据是根据x轴和y轴算出来的，如果z轴的数据不是通过x,y轴算出来的，而是它本身变量的呢？这个代码怎么写？