options ls=77 ps=58;
data tex1;
infile 'd:\b\s102\Examples\Spatial\texas.d';
input x y z;
run;
;
proc g3d;
scatter x*y=z;
run;
;
proc gplot;
bubble x*y=z;
run;
;
proc variogram outvar=tex2;
compute lagdistance=10 maxlags=20 robust;
coordinates xcoord=x ycoord=y;
directions 45,135;
var z;
run;
;
proc print data=tex2;
run;
;
proc gplot data=tex2;
symbol1 value=circle interpol=join;
plot variog*distance=angle;
plot rvario*distance=angle;
run;
;
proc reg data=tex1;
model z=x y;
output r=res1;
run;
;
proc variogram outvar=tex3;
compute lagdistance=10 maxlags=20 robust;
coordinates xcoord=x ycoord=y;
directions 45,135;
var res1;
run;
;
proc print data=tex3;
run;
;
proc gplot data=tex3;
symbol1 value=circle interpol=join;
plot variog*distance=angle;
plot rvario*distance=angle;
run;
;