Average and error bars in Matlab for a complex array

Hi,

I have a function M that is an 4-D complex variable of dimensions a*b*c*d.

Now I want to average over the last dimension only, resulting in an array Mean=a*b*c.

I understand this is done with:

 Mean=mean(A,4)

I want to show the error bars for this function, that is, the standard error.

What I have done is:

 R = sqrt(Mean/(d-1));
 Error=R/sqrt(length(d));

Now I have to rescale the parameters for my problem:

 for it=1:length(c)
     for r=0:a-1
         Err(r+1,it)=sum(diag(Error(:,:,it),r))/(a-r);
     end
 end

 for it=1:length(c)
     for r=0:a-1
         Value(r+1,it)=sum(diag(Mean(:,:,it),r))/(a-r);
     end
 end

And I plot:

 errorbar(1:a,Value(1:a,end),Err(1:a,end))

Is this correct?

Thanks in advance!

[Stephen23]

"I have a function M that is an 4-D complex array..."

A function is not an array in memory. M is likely to be a variable.

[Auryn_]

Corrected

Answer 1

[Stephen23]

MA=mean(A,4)

"Is this correct?"

Yes.