missing=sum(isnan(data)')';                                             % Location of Missing Values                        

sample=zeros(NT,3);                                                     % Sample of Complete Data        

sample(:,1)=ones(length(Y),1);sample(missing>0)=0;                      % Sample of Complete Data        
        
balanced_indi=repmat((1:N),T,1);sample(:,2)=balanced_indi(:);           % Index of Individuals
   
sample(:,3)=repmat((1:T),1,N)';                                         % Index of Time
   
Qadj=Q(sample(:,1)==1);                                                 % Values of Qit on the Adjuted Sample

Xadj=X(sample(:,1)==1,:);                                               % Values of Xit on the Adjuted Sample

sampleAdj=sample(sample(:,1)==1,:);                                     % Adjusted Sample

NT=sum(sample(:,1));                                                    % Total number of observations

% Centered Endogeneous Variable

y_cent=centered(Y,sample);                                   % Vector of Centered Data

% Centered Explicative Variables X(i,t)

for indic_k=1:K                                             % Loop on the number of explicative variable
      
    x_cent(:,indic_k)=centered(X(:,indic_k),sample);         % Centered Vector of X
   
end


if nargin==5
   
    disp('WARNING: Initial Condition Provided by User')
   
else                                               % Begin of the condition on initial condition provided by user

红色部分循环提示“带有下表的赋值维度不匹配”，求教解决方法。