Correlation coefficient calculation error

7 次查看（过去 30 天）

显示更早的评论

Josh 2022-6-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1746395-correlation-coefficient-calculation-error

评论： Mathieu NOE 2022-6-23

采纳的回答： Mathieu NOE

DQS.mat

在 MATLAB Online 中打开

As I try to find correlation coefficient for the respective variables i.e. 'td' & 'cumulat', in for loop, some error is coming in the process of execution. Please help.

clear;clc;close all
S = load('DQS.mat');
F = fieldnames(S);
for n = 1:numel(F)
    F1 = fieldnames(S.(F{n}));
    F11{n} = F1;
    % allocate the variables
    nn = cell(numel(F1),1);
    PArea = cell(numel(F1),1);
    t_max = cell(numel(F1),1);
    for n1 = 1:numel(F1)
        t = S.(F{n}).(F1{n1}).cc.t;
        t_max{n1} = max(t);
        nn{n1} = F1{n1};
        output{n1} = cumtrapz(t);
        NF = @(p,q) max(output{n1}(t<=q)) - min(output{n1}(t>=p));
        PArea{n1} = NF(5, 6);   
    end  
    T_max{n} = t_max;
    tt = cell2mat(T_max{1,n}(:,end));
    td{n} = tt./max(tt);
    cumulat{n} = PArea;
    [R{n},P{n}] = corrcoef(td,cumulat); %correlation coeff. calculation
end
Error using corrcoef
Cannot compute p-values for complex inputs.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Jonas 2022-6-23

can you please edit your question, format your code as code and press the compile button (green rectangle), this way we can see what is happening

请先登录，再进行评论。

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

采纳的回答

Mathieu NOE 2022-6-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1746395-correlation-coefficient-calculation-error#answer_991735

在 MATLAB Online 中打开

hello

seems to me the code could be made simpler, why the need to make all this cell acrobatics ?

clear;clc;close all
S = load('DQS.mat');
F = fieldnames(S);
for n = 1:numel(F)
    F1 = fieldnames(S.(F{n}));
    F11{n} = F1;
    % allocate the variables
    nn = cell(numel(F1),1);
    for n1 = 1:numel(F1)
        t = S.(F{n}).(F1{n1}).cc.t;
        t_max(n1) = max(t);
        nn{n1} = F1{n1};
        output{n1} = cumtrapz(t);
        NF = @(p,q) max(output{n1}(t<=q)) - min(output{n1}(t>=p));
        PArea(n1) = NF(5, 6);   
    end  
td = t_max./max(t_max);
cumulat = PArea;
[R{n},P{n}] = corrcoef(td,cumulat); %correlation coeff. calculation
end