function?F=carspls（XyAfoldmethodnum）?
%+++?CARS:?Competitive?Adaptive?Reweighted?Sampling?method?for?variable?selection.
%+++?X:?The?data?matrix?of?size?m?x?p
%+++?y:?The?reponse?vector?of?size?m?x?1
%+++?A:?the?maximal?principle?to?extract.
%+++?fold:?the?group?number?for?cross?validation.
%+++?num:?the??number?of?Monte?Carlo?Sampling?runs.
%+++?method:?pretreatment?method.
%+++?Hongdong?Li?Dec.15?2008.
%+++?Advisor:?Yizeng?Liang?yizeng_liang@263.net
%+++?lhdcsu@gmail.com
%+++?Ref:??Hongdong?Li?Yizeng?Liang?Qingsong?Xu?Dongsheng?Cao?Key
%????wavelengths?screening?using?competitive?adaptive?reweighted?sampling?
%????method?for?multivariate?calibration?Anal?Chim?Acta?2009?648（1）:77-84


tic;
%+++?Initial?settings.
if?nargin<6;num=50;end;
if?nargin<5;method=‘center‘;end;
if?nargin<4;fold=5;end;
if?nargin<3;A=2;end;

%+++?Initial?settings.
[MxNx]=size（X）;
A=min（[Mx?Nx?A]）;
index=1:Nx;
ratio=0.9;
r0=1;
r1=2/Nx;
Vsel=1:Nx;
Q=floor（Mx*ratio）;
W=zeros（Nxnum）;
Ratio=zeros（1num）;

%+++?Parameter?of?exponentially?decreasing?function.?
b=log（r0/r1）/（num-1）;??a=r0*exp（b）;

%+++?Main?Loop
for?iter=1:num
?????
?????perm=randperm（Mx）;???
?????Xcal=X（perm（1:Q）:）;?ycal=y（perm（1:Q））;???%+++?Monte-Carlo?Sampling.
?????
?????PLS=pls（Xcal（:Vsel）ycalAmethod）;????%+++?PLS?model
?????w=zeros（Nx1）;coef=PLS.coef_origin（1:end-1end）;
?????w（Vsel）=coef;W（:iter）=w;?
?????w=abs（w）;??????????????????????????????????%+++?weights
?????[wsindexw]=sort（-w）;??????????????????????%+++?sort?weights
?????
?????ratio=a*exp（-b*（iter+1））;??????????????????????%+++?Ratio?of?retained?variables.
?????Ratio（iter）=ratio;
?????K=round（Nx*ratio）;??
?????
?????
?????w（indexw（K+1:end））=0;??????????????????????%+++?Eliminate?some?variables?with?small?coefficients.??
?????
?????Vsel=randsample（NxNxtruew）;?????????????????%+++?Reweighted?Sampling?from?the?pool?of?retained?variables.?????????????????
?????Vsel=unique（Vsel）;??????????????
?????fprintf（‘The?%dth?variable?sampling?finished.\n‘iter）;????%+++?Screen?output.
?end

%+++??Cross-Validation?to?choose?an?optimal?subset;
RMSEP=zeros（1num）;
Q2_max=zeros（1num）;
Rpc=zeros（1num）;
for?i=1:num
???vsel=find（W（:i）~=0）;
?
???CV=plscvfold（X（:vsel）yAfoldmethod0）;??
???RMSEP（i）=CV.RMSECV;
???Q2_max（i）=CV.Q2_max;???
???
???Rpc（i）=CV.optPC;
???fprintf（‘The?%d/%dth?subset?finished.\n‘inum）;
end
[RminindexOPT]=min（RMSEP）;
Q2_max=max（Q2_max）;




%+++?save?results;
time=toc;
%+++?output
F.W=W;
F.time=time;
F.cv=RMSEP;
F.Q2_max=Q2_max;
F.minRMSECV=Rmin;
F.iterOPT=indexOPT;
F.optPC=Rpc（indexOPT）;
Ft.ratio=Ratio;
F.vsel=find（W（:indexOPT）~=0）‘;



function?sel=weightsampling_in（w）
%Bootstrap?sampling
%2007.9.6H.D.?Li.

w=w/sum（w）;
N1=length（w）;
min_sec（1）=0;?max_sec（1）=w（1）;
for?j=2:N1
???max_sec（j）=sum（w（1:j））;
???min_sec（j）=sum（w（1:j-1））;
end
%?figure;plot（max_sec‘r‘）;hold?on;plot（min_sec）;
??????
for?i=1:N1
??

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2017-03-30?21:57??CARS\
?????文件???????53760??2010-05-20?16:23??CARS\CARS_manual.doc
?????文件?????????776??2011-01-14?16:44??CARS\CIP2pred.m
?????文件???????46264??2010-04-10?10:35??CARS\LOGO_CARS.JPG
?????文件????????1132??2010-05-20?10:31??CARS\Manne_bi.m
?????文件?????????236??2016-12-15?13:43??CARS\Read?me.txt
?????文件????????3399??2011-05-21?15:53??CARS\carspls.m
?????文件?????????920??2010-08-13?10:51??CARS\classplot2.m
?????文件?????????429??2011-05-21?10:28??CARS\compute_T2DM_data.asv
?????文件??????448880??2010-02-12?14:42??CARS\corn_m51.mat
?????文件?????????278??2010-03-30?14:20??CARS\csvd.m
?????文件?????????540??2011-04-21?09:49??CARS\databin.m
?????文件????????1102??2000-10-13?00:00??CARS\dosc.M
?????文件?????????506??2010-04-16?15:52??CARS\example_nir.m
?????文件?????????465??2008-04-20?16:25??CARS\expred1.m
?????文件?????????736??2008-06-24?11:37??CARS\expred2.m
?????文件?????????812??2009-05-08?15:26??CARS\ks.m
?????文件?????????966??2010-04-12?11:10??CARS\lsreg.m
?????文件????????3031??2011-01-28?04:09??CARS\mcs.m
?????文件????????1217??2010-04-12?16:08??CARS\mcuvepls.m
?????文件?????????976??2008-12-10?08:58??CARS\mwpls.m
?????文件????????3013??2009-04-14?15:38??CARS\oscplscv.m
?????文件?????????835??2010-04-19?08:30??CARS\plotcars.m
?????文件?????????592??2010-03-08?19:07??CARS\plotmcs.m
?????文件????????1953??2010-02-14?23:57??CARS\pls.m
?????文件????????1788??2010-09-07?17:11??CARS\pls_nipals.m
?????文件????????2639??2011-01-08?18:20??CARS\plscvfold.m
?????文件????????1955??2009-12-16?11:14??CARS\plsdcv.m
?????文件????????2977??2010-04-16?15:21??CARS\plsmccv.m
?????文件????????1734??2009-10-27?08:37??CARS\plsnipals.m
?????文件????????2179??2010-06-30?11:14??CARS\plsrdcv.m
............此處省略13個文件信息