作者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