Ich mache Simulationsstudien am linearen Modell y = X * beta + eps mit der Größe (X) = [n d]. Ich betrachte den Effekt der Dimensionalität d basierend auf zwei Methoden. Ich führe 10 simulierte Daten und bekomme entsprechende Beta-Schätzungen, und dann Ich möchte den Mittelwert der Beta über die 10 simulierten Daten berechnen.Zell-Array-Mittelwert in Matlab
Mein MATLAB-Code Spielzeug ist wie folgt:
nsim=10; %iteration number
dd=[4 6]; %two dimension cases,\beta=(\beta_1,\cdots,\beta_d)^T
ddlen=length(dd);
nmethod=2; %two methods
seednum=0;
BH = cell(nsim,ddlen,nmethod); %estimation of beta w.r.t two dimension cases and two methods
for di = 1:ddlen
d = dd(di);
for ni = 1:nsim
seednum = seednum + di*ni;
randn('seed', seednum);
betahat=randn(d,1);
for method = 1:nmethod
if method==1
BH{ni,di,method} = betahat;
else
BH{ni,di,method} = 10*betahat;
end
end
end
end
Dann können wir
BH(:,:,1) =
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
BH(:,:,2) =
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
[4x1 double] [6x1 double]
erhalten möchte ich den Mittelwert über 10 Zeilen (NSIM = 10) und erhalten so etwas wie
mean(BH(:,:,1))=
[4x1 double] [6x1 double]
mean(BH(:,:,2)) =
[4x1 double] [6x1 double]
Irgendwelche Ideen? Vielen Dank!
Thanks @ EBH. Aber deine Antwort ist nicht, dass ich will. Die Rückkehr sollte zwei Vektoren sein, einer ist [4x1 double] und der andere ist [6x1 double] und mit anderen Worten, der Mittelwert von 10 [4x1 double] bzw. der Mittelwert von 10 [6x1 double]. –
Warum verwenden Sie Zellenarrays, wenn alle Ihre Arrays die gleiche Größe haben? – beaker
Warum legen Sie einen neuen Zufalls-Seed in jeder Schleife fest? – EBH