用ltrend 命令，快速求出数据的趋势线！或者给数据去掉倾向！

展望2014

运行完了，q define 会出现趋势系数

centaurus

感谢楼主！！！但是我Q DEFINE后出现x=-1,bb=-xxxxx,aa=-xxxxxxx，我画的是时间序列图，按照y=kx+b，我得到的哪个是k，哪个是b啊？按照图，k应该是正值的还是感谢楼主，看了GS，我大致知道了，只要编GS的时候不像楼主那样‘ltrend bb'而是'ltrend bb k s'就能输出系数k了，但那个s应该是个类似残差的东西，值没对应上。不大明白，得回去再琢磨琢磨

守望之月/;-D

学习~感谢LZ~~

hikingbaby

相关系数还是没出来，我再看看gs。

hikingbaby

centaurus 发表于 2014-4-13 15:07
感谢楼主！！！但是我Q DEFINE后出现x=-1,bb=-xxxxx,aa=-xxxxxxx，我画的是时间序列图，按照y=kx+b，我得到 ...

k输出是个负值？这个是相关系数r吗？还是斜率？

centaurus

hikingbaby 发表于 2014-4-19 10:02
k输出是个负值？这个是相关系数r吗？还是斜率？

不好意思我是分了两次回复结果中间忘加标点符号了。你编脚本输入‘ltrend bb k s'的话k是斜率，s估计是残差平方和。你找几个点算一下应该是这样的，反正我算的是的

Pynomid

请问一下楼主，我画图的时候有的数据能输出趋势线和斜率，
但是有的就不行，提示[qdim info] calendar mode not set yet,standard calendar is assumed
这是什么情况呢？谢谢

liuzf406

有知道如何检验趋势显著性的吗？

/xin村儿/

楼主是否验证过举例的正确性？
个人觉得这两行貌似有点问题，求趋势的时候，litrend是如何循环的1947-2009的Jan,Feb,.....Dce？
'set time   Jan1947   Dec2009'
'ltrend  ssta  bb'  
盼指点，谢谢！

/xin村儿/

兰溪之水 发表于 2013-7-8 12:34
使用这个还需要qdims.gsf这个文件，可以到这里下：
http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts ...

兰溪师兄，我感觉楼主例子中这两行貌似有点小问题，担忧具体不知道怎么修改！
就是求趋势的时候，litrend是如何循环的1947-2009的Jan,Feb,.....Dce？
'set time   Jan1947   Dec2009'
'ltrend  ssta  bb'  
盼指点，谢谢！
比如，我用您的grads2ascii读写如下，结果txtl里面竟然不是12个趋势值，而是56行*12个，即56年*12月
'sdfopen E:/sst.mnmean.nc'
'set lon 315'
'set lat 55'
'set time jan1958 dec2013'
'ltrend  sst  trend'
'grads2ascii trend E:/sst-trend-5801-1312.txt %13.5f 12 0'