空間と時間におけるプロセス(Y)の進化を説明する次の確率モデルがあります。DsとDtは、空間(x軸とy軸の2D)と時間(t軸の1D)のドメインです。このモデルは通常、混合効果モデルまたは変動成分モデルとして知られています。
私は現在、次のようにYを開発しています。
%# Time parameters
T=1:1:20; % input
nT=numel(T);
%# Grid and model parameters
nRow=100;
nCol=100;
[Grid.Nx,Grid.Ny,Grid.Nt] = meshgrid(1:1:nCol,1:1:nRow,T);
xPower=0.1;
tPower=1;
noisePower=1;
detConstant=1;
deterministic_mu = detConstant.*(((Grid.Nt).^tPower)./((Grid.Nx).^xPower));
beta_s = randn(nRow,nCol); % mean-zero random effect representing location specific variability common to all times
gammaTemp = randn(nT,1);
for t = 1:nT
gamma_t(:,:,t) = repmat(gammaTemp(t),nRow,nCol); % mean-zero random effect representing time specific variability common to all locations
end
var=0.1;% noise has variance = 0.1
for t=1:nT
kappa_st(:,:,t) = sqrt(var)*randn(nRow,nCol);
end
for t=1:nT
Y(:,:,t) = deterministic_mu(:,:,t) + beta_s + gamma_t(:,:,t) + kappa_st(:,:,t);
end
私の質問は次のとおりです。
- Yの式でデルタを生成し、カッパとデルタの差を生成するにはどうすればよいですか?
- 私がYを正しく生成しているかどうかを、Matlabを使用したいくつかの図を通して説明するのを手伝ってください。
さらに情報や説明が必要な場合はお知らせください。ありがとう。