function [Aprederror,MSE,pcnterr,Apred] = modelvalidate(P,A)% modelvalidate.m Compare wafer profile predictions from polymodel%                 to a set of experimental data. Called as%%                 [Aprederror,MSE,pcnterr,Apred] = modelvalidate(P,A)%% INPUT PARAMETERS%%                P : a polymodel object%                A : vector of wafermap objects%% OUTPUT PARAMETERS%%       Aprederror : vector of wafermap objects%              MSE : vector of mean squared errors%          pcnterr : vector of prediction errors, %%            Apred : vector of predicted wafermap objects% Copyright (c) by Raymond A. Adomaitis, 1998-2006% the polymodel parameterspolymodelparam = get(P,'param');Apred = A;Aprederror = A;for i = 1:length(A)% get parameters from each wafermap object    Ap   = get(A(i),'param');    Akey = keys(Ap);% now reset the parameter values in the polymodel according% to the wafer param values    for j = 1:length(Akey)        P = set(P,'param',getval(Ap,Akey{j}),Akey{j});    end    [P,Apred(i)] = residual(P);    Aprederror(i) = A(i) - Apred(i);    MSE(i) = sqrt(mean(Aprederror(i)^2));    pcnterr(i) = mean(abs(Aprederror(i)))/mean(A(i);end