Taesam Lee

function [fcond,Pcond,lambda,yi]=NPcond(y,X,Xc,yi)
%
% Code written by Taesam Lee and Dominique Fasbender 04/03/2009 : INRS-ETE Quebec
% Nonparametric Condtional Distribution
% X : data of condtion variable (numvar*length), y : current variable(length)
% Xc: condtion values
% yi: estimated set for fcond
% example
% t=mvnrnd([0,0],[.9 .4; .4 .3],100);
% X=t(:,1)';y=t(:,2)';yi=-2:0.01:2;Xc=[1];
% [fcond,Pcond,yi]=NPcond(y,X,Xc,yi);
% plot(yi,fcond);plot(yi,Pcond);

Xall=[y;X];
sz=size(Xall);d=sz(1);n=sz(2);

% Estimate the bandwidth ratio
lambda_ref = (4/(d+2))^(1/(d+4))*n^(-1/(d+4));
%min_lambda = 0.25*lambda_ref;max_lambda = 1.1*lambda_ref;
%lambda = fminbnd(@(lambda)LSCV(lambda,Xall',n,d),min_lambda,max_lambda);
lambda=lambda_ref;

% Covariance
Sall=cov(Xall');
Syy=Sall(1,1);Sxy=Sall(2:d,1);Sxx=Sall(2:d,2:d);
Sd=Syy-Sxy'*inv(Sxx)*Sxy;

% Conditional Distribution
difX=repmat(Xc,1,n)-X;% need to be changed
s1=inv(Sxx)
%for in=1:n
% dmat1(in)=(difX(:,in))'*s1*(difX(:,in));
%end
%es=exp(-1/2*lambda_ref^2*dmat1);

es2 = 0;
for i=size(s1,1)
for j=size(s1,2)
es2 = es2 + exp(-1/2*lambda_ref^2*difX(i,:).*difX(i,:)*s1(i,j));
end
end

wi=es2./sum(es2);
bi=Sxy'*inv(Sxx)*difX;

for it=1:length(yi)
fcond(it)=sum(1/sqrt(2*pi*lambda^2*Sd)*wi.*exp(-(yi(it)-bi).^2./(2*lambda^2*Sd)));
Pcond(it)=sum(wi.*normcdf(yi(it),bi,sqrt(lambda^2*Sd)));
end

Hello,
Can you post the data that you used for plotting the graphic UH.jpg?
It would with good to test your code with your data.

I can't get this to work.... the script executes successfully, in that I don't get any error messages, but returns an empty array of Z[]... I have posted a question on matlab central to ask for some advice... any suggestions for debugging would be helpful, thanks.