matlab - matlab を使用したビニング後の RMSE 計算

Question

こんにちは、RMSE を計算したいのですが、論理的な間違いを犯しているようです。これは以下のコードです。

topEdge = 20; % upper limit 
botEdge = 5; % lower limit
numBins = 40;   % define number of bins 
b_vector = speed; a_vector = power;
[N,edges,bins] = histcounts(b_vector,numBins);
Pow_means = [];
speed_means = [];

  For n = 1:numBins
    Pow_means(n,1) = mean(a_vector(bins==n,1));  % for each bins mean value calculation.
   speed_means(n,1) = mean(b_vector(bins==n,1));  % for each bins mean value calculation.
   pow_std(n,1) = std(a_vector(bins==n,1));       % for Standard deviation calculation
  end

residual = power -Pow_means; RMSE = rms(residual); variation = RMSE/mean(power);

RMSE を計算する私のアプローチは正しいですか?? またはRMSEを計算するためにどのような修正を行う必要がありますか