How to find the coefficient of determination (R2) between 2 variables in netcdf?

Question

Keegan Carvalho 2018-12-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/433715-how-to-find-the-coefficient-of-determination-r2-between-2-variables-in-netcdf

编辑： Keegan Carvalho 2018-12-4

ssh.zip

I have attached a zipped netcdf file ('ssh.nc').

I wanted to do a corellation test (OR to find the R2 value between two variables). The two variables are 'zos' and 'bottomT' which have the following dimensions of longitude,latitude,time. After this I wanted to plot a map so I can find out in which regions the correlation is high/low.

I tried to use ncread and extract grid-by-grid values but it is time- consuming and I might be prone to making errors. Is there any code that can find the correlation between two variables?

Would be grateful for assistance in this.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

KSSV 2018-12-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/433715-how-to-find-the-coefficient-of-determination-r2-between-2-variables-in-netcdf#answer_350442

在 MATLAB Online 中打开

It depends whether you want to find it spatially or temporally. I feel temporal regression is the one you are looking for. Check this code:

ncfile = 'ssh.nc' ;
A = ncread(ncfile,'zos') ;
B = ncread(ncfile,'bottomT') ;
[nx,ny,nt] = size(A) ;
R = zeros(nx,ny) ;
for i = 1:nx
    for j = 1:ny
        Ai = squeeze(A(i,j,:)) ;
        Bi = squeeze(B(i,j,:)) ;
        R(i,j) = regression(Ai',Bi')^2 ;
    end
end