% Newer breast cancer data set from KFSYSCC: 158 samples with ER, Nodes, 
% 	and survival information


% startup if cd'd here

addpath('../data')

load SYSBC.158

whos


[ind,sER]=select_genes_cor(x,C_ER,100); 
		 	 		% vector sERG holds correlations
[a,j]=sort(C_ER); 
imagecf(x(ind,j),'gr')
b=sum(C_ER);  line(0.5+[b b],[0 101])			   % separates +/- arrays
c=0.5+sum(sER(ind)>0); line([0 n+1],[c c]+.5);      % separates pos/neg genes
cmap jet
cmap hot


% ----
% explore ER gene  predictors ... 
%
[beta,v,prob,haty,e,varind]=mregauto(y,x(ind(2:20),:),.05);

varind

beta

stepwise(x(ind(2:20),:)',y)  % click out all but varind .... 


Des(ind(varind),:)

				% click out more?