用R绘制交互式3D图

1. 安装rgl包、evd包和MASS包：

1. > install.packages("rgl")  
2. Installing package(s) into ‘D:/program/R/library’  
3. (as ‘lib’ is unspecified)  
4. 试开URL’'  
5. Content type 'application/zip' length 1937097 bytes (1.8 Mb)  
6. 打开了URL  
7. downloaded 1.8 Mb  
8.   
9. package ‘rgl’ successfully unpacked and MD5 sums checked  
10.   
11. The downloaded binary packages are in  
12.     C:\Users\LiQiang\AppData\Local\Temp\RtmpIjwxP2\downloaded_packages  
13. > install.packages("evd")  
14. Installing package(s) into ‘D:/program/R/library’  
15. (as ‘lib’ is unspecified)  
16. 试开URL’'  
17. Content type 'application/zip' length 1176009 bytes (1.1 Mb)  
18. 打开了URL  
19. downloaded 1.1 Mb  
20.   
21. package ‘evd’ successfully unpacked and MD5 sums checked  
22.   
23. The downloaded binary packages are in  
24. 　　C:\Users\LiQiang\AppData\Local\Temp\RtmpIjwxP2\downloaded_packages  
25. > install.packages("MASS")  
26. Installing package(s) into ‘D:/program/R/library’  
27. (as ‘lib’ is unspecified)  
28. 试开URL’'  
29. Content type 'application/zip' length 872270 bytes (851 Kb)  
30. 打开了URL  
31. downloaded 851 Kb  
32.   
33. package ‘MASS’ successfully unpacked and MD5 sums checked  
34.   
35. The downloaded binary packages are in  
36.     C:\Users\LiQiang\AppData\Local\Temp\RtmpIjwxP2\downloaded_packages  

2. 输入程序

1. library(rgl)  
2. library(MASS)  
3. library(evd);  
4. data(lossalae)  
5. U=rand(lossalae[,1]+rnorm(nrow(lossalae),mean=0,sd=.001))/(nrow(lossalae)+1)  
6. V=rank(lossalae[,2])/(nrow(lossalae)+1)  
7. M=kde2d(qnorm(U), qnorm(V), n=35)  
8. persp3d(M$x, M$y, M$z, col='green',xlab="loss", ylab="alae", zlab="")  

3. 输出3D图