求助，利用GRADS将满足条件的值用点绘制出来，出现问题，请高手帮忙

风儿飘飘

最近做回归分析，将通过显著性检验的点标记为黑点出错，查了好多资料，也理不清，请高手帮忙！这是我的gs文件！
出图结果附后1
*OLR regression coefficient'reinit''open F:/REGanalysis/Regressdata/OLR/regolr.ctl''open F:/REGanalysis/Regressdata/OLR/Ffc.ctl''enable print F:/REGanalysis/Regressdata/OLR/REGolr.gmf''set z 1''set t 1''set parea 0.75 10. 0.5 7.5''set grads off''set grid off''set mpdset mres'*'set mpdset cnworld''set gxout shaded''set csmooth on''set annot 1 8''set xlopts 1 5 0.15 ''set ylopts 1 5 0.15 ''set xlint 10''set ylint 10''d regxs''cbarn 1 1 10.25 4. 'lon1=60lat1=0while(lon1<=160&lat1<=60)'set lon 'lon1'''set lat 'lat1''if(f.2>=4.10)'q w2xy 'lon1' 'lat1''x1=subwrd(result,3)y1=subwrd(result,6)*say 'x1='x1*say 'y1='y1 'set cmax -100''d f.2''draw mark 3 'x1' 'y1' 3'elsebreak endiflon1=lon1+1lat1=lat1+1endwhile*'set xlpos -20'   *'set ylpos -20'*'set black -0.3 0.3'*'set rbcols 4'*'set ccolor 15''draw title \\Regression of the AOD left coefficient onto OLR ''printim F:/REGanalysis/Regressdata/OLR/REGolr.png X800 Y600 white' 'print''disable print' ;

风儿飘飘

*OLR regression coefficient
'reinit'
'open F:/REGanalysis/Regressdata/OLR/regolr.ctl'
'open F:/REGanalysis/Regressdata/OLR/Ffc.ctl'
'enable print F:/REGanalysis/Regressdata/OLR/REGolr.gmf'
'set z 1'
'set t 1'
'set parea 0.75 10. 0.5 7.5'
'set grads off'
'set grid off'
'set mpdset mres'
*'set mpdset cnworld'
'set gxout shaded'
'set csmooth on'
'set annot 1 8'
'set xlopts 1 5 0.15 '
'set ylopts 1 5 0.15 '
'set xlint 10'
'set ylint 10'
'd regxs'
'cbarn 1 1 10.25 4. '
lon1=60
lat1=0
while(lon1<=160&lat1<=60)
'set lon 'lon1''
'set lat 'lat1''
if(f.2>=4.10)
'q w2xy 'lon1' 'lat1''
x1=subwrd(result,3)
y1=subwrd(result,6)
*say 'x1='x1
*say 'y1='y1
'set cmax -100'
'd f.2'
'draw mark 3 'x1' 'y1' 3'
else
break
endif
lon1=lon1+1
lat1=lat1+1
endwhile
*'set xlpos -20'   
*'set ylpos -20'
*'set black -0.3 0.3'
*'set rbcols 4'
*'set ccolor 15'
'draw title \\Regression of the AOD left coefficient onto OLR '
'printim F:/REGanalysis/Regressdata/OLR/REGolr.png X800 Y600 white'
'print'
'disable print'
;

风儿飘飘

前面发的gs文件太乱，重新发了一下，请哪位高人给指点下，感谢

river

风儿飘飘 发表于 2013-7-29 18:28

                               
登录/注册后可看大图

*OLR regression coefficient
'reinit'
'open F:/REGanalysis/Regressdata/OLR/regolr.ctl'

'set cmax -100'，不是已经有了if(f.2>=4.10)的限制，为什么还要限制不画大于-100的值呢，感觉是矛盾啊

风儿飘飘

原来我对set cmax 的理解不够深，那么只要设定set cmax 4.1行吗，不需要循环吗

风儿飘飘

刚才又试过了不行

风儿飘飘

求助同学解决了此问题！

fairyes

风儿飘飘 发表于 2013-7-30 12:23

                               
登录/注册后可看大图

求助同学解决了此问题！

请问下最后怎么解决的，我现在也需要把显著的地方标记出来