clear
set obs 900
gen x = int((_n - mod(_n-1,30) -1 ) /30 )
gen y = mod(_n-1,30)
gen z1 = -0.491-0.116*x-0.133*y+0.16*x*y
surface x y z1
gr save mygraph1,replace
clear
set obs 900
gen x = int((_n - mod(_n-1,30) -1 ) /30 )
gen y = mod(_n-1,30)
gen z2 = 0.491-0.102*x+0.369*y-0.544*x*y
surface x y z2
gr save mygraph2,replace
clear
set obs 900
gen x = int((_n - mod(_n-1,30) -1 ) /30 )
gen y = mod(_n-1,30)
gen z3 = -0.491-0.186*x-0.133*y-0.077*x*y
surface x y z3
gr save mygraph3,replace
clear
set obs 900
gen x = int((_n - mod(_n-1,30) -1 ) /30 )
gen y = mod(_n-1,30)
gen z4 = 0.491-0.074*x+0.369*y+0.599*x*y
surface x y z4
gr save mygraph4,replace
gr combine mygraph1.gph mygraph2.gph mygraph3.gph mygraph4.gph,replace
*试一下:
clear all
set obs 900
g x = int((_n - mod(_n-1,30) -1 ) /30 )
g y = mod(_n-1,30)
g z1 = -0.491-0.116*x-0.133*y+0.16*x*y
g z2 = 0.491-0.102*x+0.369*y-0.544*x*y
g z3 = -0.491-0.186*x-0.133*y-0.077*x*y
g z4 = 0.491-0.074*x+0.369*y+0.599*x*y
forv i=1/4{
surface x y z`i'
gr save mygraph`i',replace
}
gr combine mygraph1.gph mygraph2.gph mygraph3.gph mygraph4.gph
sungmoo 发表于 2010-3-2 18:15
*试一下:
clear all
set obs 900
g x = int((_n - mod(_n-1,30) -1 ) /30 )
g y = mod(_n-1,30)
g z1 = -0.491-0.116*x-0.133*y+0.16*x*y
g z2 = 0.491-0.102*x+0.369*y-0.544*x*y
g z3 = -0.491-0.186*x-0.133*y-0.077*x*y
g z4 = 0.491-0.074*x+0.369*y+0.599*x*y
forv i=1/4{
surface x y z`i'
gr save mygraph`i',replace
}
gr combine mygraph1.gph mygraph2.gph mygraph3.gph mygraph4.gph
. gr combine mygraph1.gph mygraph2.gph mygraph3.gph mygraph4.gph
file mygraph1.gph is not a Stata live graph
r(198);
显示这样的错误结果。。。。
如果是2维的应该可以。。
你把surface命令升级一下
应该是2008年的版本
. which surface
D:\programs\Stata10\ado\stbplus\s\surface.ado
*! Date : 14 July 2008
*! Version : 1.01
*! Author : Adrian Mander
*! Email : adrian.mander@mrc-hnr.cam.ac.uk