线性回归问题

推开世界的门

各位大佬，非常抱歉。叨扰大家，我想请问一下如何求两组数据的回归方程呢？，如下列两组数据a1和a2，时间序列都是1950-1990，求一个线性回归方程y=kx+b,其中k应该为两组数据的相关性，那么b用哪个函数可以确定呢，数组比较多，所以如果ncl处理更方便一点。

a1	27.96774	a2	1.801017
	27.91092		1.419363
	27.90236		1.535143
	27.93173		1.66675
	28.00713		1.961259
	27.83422		1.117052
	27.87683		1.510553
	28.06966		2.067458
	27.92752		1.875387
	28.05965		1.699253
	28.15343		1.324489
	28.04453		2.082381
	28.15018		1.828656
	28.32475		2.767909
	28.1543		1.83427
	28.03699		2.297135
	28.30981		1.901899
	28.27892		2.310006
	28.1358		2.65395
	28.07014		1.857434
	28.06121		0.9075552
	28.0159		1.273672
	27.70146		2.693526
	27.76564		0.6403334
	28.02465		2.289139
	27.70432		0.9882634
	27.6315		2.049191
	27.90511		1.616049
	27.8028		1.408753
	27.88092		1.762757
	27.88425		2.36526
	28.00392		1.591004
	27.73189		3.020699
	27.50834		1.052505
	28.20459		1.736361
	28.03667		1.774717
	27.99497		2.865948
	27.76427		2.032925
	28.04434		1.158222
	28.08507		2.105983
	28.11559		2.762481

topmad

为啥要用ncl呢
python或者R语言 一行就出结果了啊

推开世界的门

topmad 发表于 2021-10-19 23:48
为啥要用ncl呢
python或者R语言 一行就出结果了啊

好的 谢谢您

zhoushiyang123

https://www.ncl.ucar.edu/Applications/regress.shtml

guoguohh

rc    = regline(a1,a2)
print(rc)
得到的xave、yave是平均值；
tval是t统计量；
rstd是回归系数的标准差；
yintercept是y截距，也就是你要的b；
nptxy是样本个数；
最下面的那个数值就是你要的斜率k

推开世界的门

guoguohh 发表于 2021-10-27 09:55
rc    = regline(a1,a2)
print(rc)
得到的xave、yave是平均值；

谢谢大佬，如果对于一个三维数组，怎么才能更方便的求出k和b呀，如我想对一个三维（时间|50，纬度|121，经度|121）的nc文件的每个格点求一个回归一元线性回归方程y=kx+b，x为用其他数据定义的自变量，请问一下怎么才能将每个格点的b算出来呢，下面是我写的代码，如果只算每个格点的k值的话，可以算出来，但是不知道值该怎么才能一次性将所有格点的b值算出来k=new((/120,120/),float)
do i=0,119
  do j=0,119
    k(i,j)   = regline (x,y(:,i,j))
  end do
end do

print(k)

guoguohh

推开世界的门 发表于 2021-10-27 11:46
谢谢大佬，如果对于一个三维数组，怎么才能更方便的求出k和b呀，如我想对一个三维（时间|50，纬度|121， ...

https://www.ncl.ucar.edu/Document/Functions/Built-in/regCoef_n.shtml
你参考一下官网的这个函数，多维的不用写循环，看每个的b值直接print(变量名@yintercept)就得到了

推开世界的门

guoguohh 发表于 2021-10-27 14:47
https://www.ncl.ucar.edu/Document/Functions/Built-in/regCoef_n.shtml
你参考一下官网的这个函数，多 ...

好的 大佬！ 谢谢您！问题已解决

guoguohh

推开世界的门 发表于 2021-10-27 16:28
好的 大佬！ 谢谢您！问题已解决

没事哦，能帮到你就好~