grads能否画出线性回归的图

shirly

发现用grads画图时，总会有很常用但是无法实现的，具体如下：

1、之前用grads作图，做出过某一点随时间的变化率，即一元线性回归的系数b（y=bx+a），可以画出这个要素随时间的变化曲线图，但是要怎么画那条拟合线？

2、有某一点的两个要素a1、a2（a1和a2接近，存在线性关系），想要用grads画出一张横坐标为a1，纵坐标为a2的点图（如附件中图1），再画出一条拟合线，和问题1类似，要怎么画出那条拟合线？

如果哪位知道怎么解决（尤其问题2），还望指教，十分感谢~！

bjj948

有个Bin Guan's GrADS Script Library
http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html

可能有相关gs文件。

river

第一个比较好实现，使用scorr或者tcorr 来实现，第二个问题和第一个基本一样。还想说的一点不是grads实现不了，只能说用grads不能一步到位，我个人觉得用其他软件也方便不到哪去。楼主说实现不了，只能说明楼主对于grads的掌握程度不够，也许基础的东西掌握很好，但是不会灵活运用······

shirly

river 发表于 2014-3-22 10:36

                               
登录/注册后可看大图

第一个比较好实现，使用scorr或者tcorr 来实现，第二个问题和第一个基本一样。还想说的一点不是grads实现不 ...

谢谢river的回复，我确实用grads没那么灵活，所以来这里求指导，感谢你的指导，我先试试~出问题的地方还请多指教~

展望2014

我知道现在可以加拟合线了，去这个网站http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html，下载ltrend.gs和qdims.gsf这个脚本，就可以实现了。你先把这两个脚本放到存脚本的文件夹里。在自己画图的gs文件里加上'ltrend <input> [<output> [<slope> [<rmse>]' i，nput 是要画的变量(y轴是变量，x轴是时间)，output改成display；lope和rmse分别是回归系数和均方差，可以加上也可以不加上，需要的话就写上slope和rmse就可以，运行结束后可以用q define命令，grads会在命令框里告诉你这两个参数的值。

shirly

展望2014 发表于 2014-3-22 21:54

                               
登录/注册后可看大图

我知道现在可以加拟合线了，去这个网站http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html，下载lt ...

特别感谢你能花时间给我回复帖子，我也是刚刚看到家园有个帖子，说了下载ltrend.gs的这个方法，刚下载了，其实现在需要解决的是怎么解决横坐标不是时间的问题，我再查查，再次感谢{:loveliness:}

展望2014

shirly 发表于 2014-3-22 22:03

                               
登录/注册后可看大图

特别感谢你能花时间给我回复帖子，我也是刚刚看到家园有个帖子，说了下载ltrend.gs的这个方 ...

不客气啦，我是新手，大家相互学习嘛。要是你能弄懂了求@

shirly

bjj948 发表于 2014-3-22 08:42

                               
登录/注册后可看大图

有个Bin Guan's GrADS Script Library
http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html

感谢

heyezi

谢谢，查找用grads计算rmse,果真在这找到链接了。

5xx