作者carlos0987 (carol)
看板Statistics
标题[程式] 关於R的3D绘图问题
时间Tue Apr 16 01:20:36 2013
各位大大好，我卡了快一个星期，才厚着脸皮来到这项各位高手求救，因为我是统计软体
新手，之前是英文系的同学，所以完全没用过统计软体，但我自己有摸索R软体快三个月
，我还是不大懂如何画3D立体图形，拜托请各位大大教我如何解决我现在遇到的问题，
小弟感激不尽!
问题1.
请问X/Y在CSV档案中，是表示X轴为1-~1，而Y轴为1-~1 ,请问我输入指令之後，R软体也
会是认为X轴为-1~1而Y轴为-1~1吗?下图的80表为X轴为-1，Y轴为-1时的值。
如下所示:
X/Y -1 -0.9 . . . . . 1
-1 80 . . . . . . .
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
1 . . . . . . . 90
问题2.
请问用这样的档案如何画成3D图档呢?小弟已经卡了快一星期了...请问可以求助各位高手
帮我解答吗?谢谢各位高手。
我想画成X轴为-1~1，而Y轴为-1~1，Z轴为X跟Y轴的值。
这是我的指令，如下:
na.df<-read.table("fig3.csv", header=T, row.name=NULL, sep=",", dec=".")
attach(na.df)
x<-y<-seq(-1,1,length=202)
f<-function(x,y){(exp(-0.5*x^2-0.5*y^2))/(2*pi)}
z<-outer(x,y,f)
persp(na.df,
theta=35,phi=25,col='lightblue',ylab = expression(paste(phi[2])))
title(main = "Exchange rare volatility",
ylab = expression(paste(phi[2])),
xlab = expression(paste(phi[1])),
zlab = "Variance")
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 134.208.61.165
※ 编辑: carlos0987 来自: 134.208.61.165 (04/16 01:22)
1F：→ andrew43:persp()的资料引入错了. 要persp(x,y,z,...) 04/16 01:48
2F：→ andrew43:自己看一下help(persp)的例子. 04/16 01:49
3F：→ carlos0987:是的，我改了，图形能跑，但是我这样的CODE没有问题吗? 04/16 01:56
4F：→ carlos0987:不过还是谢谢你热心答覆我XD 04/16 02:07
5F：→ andrew43:你既不知道你的问题是什麽, 我当然也不知道. 04/16 02:33
6F：→ andrew43:请明确指出问题是什麽, 别人才方便帮助你. 04/16 02:34
7F：→ carlos0987:我的问题是CSV档这样的资料，r读进去也是跟csv档所表示 04/16 08:56
8F：→ carlos0987:的这样吗?还有我用ATTACH的指令是对的吗? 04/16 08:57
9F：→ andrew43:不知道你的csv内容所以不能确定. 但这不是自己试试就知了 04/17 21:06
10F：→ andrew43:吗? 04/17 21:06
11F：→ andrew43:如果不知道某function的功能就用help()了解. 04/17 21:11
12F：→ andrew43:再不行可以买本r的工具书, 不然你这样学下去会很辛苦. 04/17 21:12