Error in Variables Regression

Question

Fabrice Lambert 2023-9-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2015981-error-in-variables-regression

编辑： Fabrice Lambert 2023-9-4

Hi. I'd like to perform a multilinear Error-in-Variables Regression, i.e. I have uncertainty in the predictors as well as in the dependent variable. This is sometimes also named Total Least Squares Regression. To my surprise, this is not part of the standard matlab functions as far as I can tell. Does anyone know if I missed something or if this (rather basic functionality I would say) is indeed not included with Matlab?

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

请先登录，再进行评论。

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

Answer 1

Torsten 2023-9-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2015981-error-in-variables-regression#answer_1299791

编辑：Torsten 2023-9-2

Linear Total Least Squares Regression means minimizing the sum of distances of projections of your data on a linear subspace. This is usually accomplished by using "svd".

If you don't want to adapt "svd" according to your need, you can try the followig package:

https://uk.mathworks.com/matlabcentral/fileexchange/31109-total-least-squares-method

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

Torsten 2023-9-2

https://www.youtube.com/watch?v=u2f4LG6OE5o

might help.

Fabrice Lambert 2023-9-4

编辑：Fabrice Lambert 2023-9-4

在 MATLAB Online 中打开

Thanks.

I modified his code to make it a bit more flexible in terms of dimensions in case anyone's interested.

function [Crit_TLS, ThetaEstimTLS] = tls(Pred,Crit,Predstd,Critstd)
%tls Total Least Squares Regression
%   Pred: Predictors in column form
%   Crit: Criterions in column form
%   Predstd: standard deviation of Predictors, horizontal vector
%   Critstd: standard deviation of Criterions, horizontal vector.
%   Modified from Antonio Sala Piqueras, Universitat Politècnica de València.
if size(Pred,1) ~= size(Crit,1)
    error('X and Y must have the same number of rows')
end
Crit_N = size(Crit,2);
Pred_N = size(Pred,2);
Crit_std = diag(Critstd);
Pred_std = diag(Predstd);
Crit_scaled = (Crit-mean(Crit)) / Crit_std; % scale to mean zero (data) and unit variance (noise)
Pred_scaled = (Pred-mean(Pred)) / Pred_std; % scale to mean zero (data) and unit variance (noise)
Data_scaled=[Crit_scaled Pred_scaled];
[N,~]=size(Data_scaled);
[~,~,V]=svd(Data_scaled / sqrt(N-1),'econ'); %TLS and SVD are the same with the proposed scaling.
Model_scaled = V(:,Pred_N+1:end);
ModPred_sc = Model_scaled(end-Pred_N+1:end,:); 
ModCrit_sc = Model_scaled(1:Crit_N,:);
ModCrit = Crit_std \ ModCrit_sc;
ModPred = Pred_std \ ModPred_sc;
ThetaEstimTLS = -ModPred / ModCrit;
Crit_TLS = Pred * ThetaEstimTLS;
end