1. 安装rgl包、evd包和MASS包:
- > install.packages("rgl")
- Installing package(s) into ‘D:/program/R/library’
- (as ‘lib’ is unspecified)
- 试开URL’http://cran.dataguru.cn/bin/windows/contrib/2.15/rgl_0.92.892.zip'
- Content type 'application/zip' length 1937097 bytes (1.8 Mb)
- 打开了URL
- downloaded 1.8 Mb
- package ‘rgl’ successfully unpacked and MD5 sums checked
- The downloaded binary packages are in
- C:\Users\LiQiang\AppData\Local\Temp\RtmpIjwxP2\downloaded_packages
- > install.packages("evd")
- Installing package(s) into ‘D:/program/R/library’
- (as ‘lib’ is unspecified)
- 试开URL’http://cran.dataguru.cn/bin/windows/contrib/2.15/evd_2.3-0.zip'
- Content type 'application/zip' length 1176009 bytes (1.1 Mb)
- 打开了URL
- downloaded 1.1 Mb
- package ‘evd’ successfully unpacked and MD5 sums checked
- The downloaded binary packages are in
- C:\Users\LiQiang\AppData\Local\Temp\RtmpIjwxP2\downloaded_packages
- > install.packages("MASS")
- Installing package(s) into ‘D:/program/R/library’
- (as ‘lib’ is unspecified)
- 试开URL’http://cran.dataguru.cn/bin/windows/contrib/2.15/MASS_7.3-21.zip'
- Content type 'application/zip' length 872270 bytes (851 Kb)
- 打开了URL
- downloaded 851 Kb
- package ‘MASS’ successfully unpacked and MD5 sums checked
- The downloaded binary packages are in
- C:\Users\LiQiang\AppData\Local\Temp\RtmpIjwxP2\downloaded_packages
2. 输入程序
- library(rgl)
- library(MASS)
- library(evd);
- data(lossalae)
- U=rand(lossalae[,1]+rnorm(nrow(lossalae),mean=0,sd=.001))/(nrow(lossalae)+1)
- V=rank(lossalae[,2])/(nrow(lossalae)+1)
- M=kde2d(qnorm(U), qnorm(V), n=35)
- persp3d(M$x, M$y, M$z, col='green',xlab="loss", ylab="alae", zlab="")
3. 输出3D图