matlab iPLSToolbox.zip

matlab ipls少有。全面，适合化学计量学方向专业人士使用。

function biModel=bipls（XYno_of_lvprepro_methodintervalsxaxislabelsval_methodsegments）

%  bipls: backwards elimination of non-informative intervals the Leardi way...
%
%  Input:
%  X is the independent variables
%  Y is the dependent variable（s） NOTE: Y is allways autoscaled
%  no_of_lv is the maximum number of PLS components
%  prepro_method （for X only） is ‘mean‘ ‘auto‘ ‘mscmean‘ or ‘mscauto‘
%     Note: msc is performed in each interval
%  intervals is the number of intervals
%     if intervals is a row vector divisions are made based on the elements
%     [startint1 endint1 startint2 endint2 startint3 endint3] see an example in manint
%  xaxislabels （self explainable） if not available type []
%  val_method is ‘test‘ ‘full‘ ‘syst111‘ ‘syst123‘ ‘random‘ or
%     ‘manual‘; the last five are cross validation based methods
%  segments （segments = number of samples corresponds to full cv）
%     if intervals is a cell array cross validation is performed according
%     to this array see the script makeManualSegments
%
%  Output:
%  biModel is a structured array containing all model information
%
%  Copyright Chemometrics Group - KVL Copenhagen Denmark
%  Lars N鴕gaard July 2004
%
%  biModel=bipls（XYno_of_lvprepro_methodintervalsxaxislabelsval_methodsegments）;

%  Functions used: iPLSreverse

if nargin==0
    disp（‘ ‘）
    disp（‘ biModel=bipls（XYno_of_lvprepro_methodintervalsxaxislabelsval_methodsegments）;‘）
    disp（‘ ‘）
    disp（‘ Example:‘）
    disp（‘ biModel=bipls（XY10‘‘mean‘‘20[]‘‘syst123‘‘5）;‘）
    disp（‘ ‘）
    return
end

% Error checks
if ~ismember（val_method{‘test‘ ‘full‘ ‘syst123‘ ‘syst111‘ ‘random‘ ‘manual‘}）
    disp（‘ Not allowed validation method‘）
    biModel=[];
    return
end

if ~ismember（prepro_method{‘mean‘ ‘auto‘ ‘mscmean‘ ‘mscauto‘ ‘none‘}）
    disp（‘ Not allowed preprocessing method‘）
    biModel=[];
    return
end
% End error checks

if strcmp（val_method‘full‘）
    segments=size（X1）;
end

ModelReverse=sub_iPLSreverse（XYno_of_lvprepro_methodintervalsxaxislabelsval_methodsegments）;
[temp1temp2minRMSEwithout（1）ix_for_iterative（1）minRMSEglobal（1）]=sort_ipls（ModelReverse）;  % Subfunction of this file

newX=X;
RevVars（1）=size（newX2）;
newX（:ModelReverse.allint（ix_for_iterative（1）2）:ModelReverse.allint（ix_for_iterative（1）3））=[];

keeptrackofinterval=[（1:intervals）‘ （1:intervals）‘];
RevIntInfo（1:）=keeptrackofinterval（ix_for_iterative（1）:）;
keeptrackofinterval（ix_for_iterative（1）:）=[];
keeptrackofinterval（:1）=（1:（intervals-1））‘;
%IntervalInformation=[（1:intervals）‘ ModelReverse.allint（1:intervals:）];
for i=2:（intervals-1）
    %plot（newX‘）pause
    RevVars（i）=size（newX2）;
    ModelReverse=sub_iPLSreverse（newXYno_of_lvprepro_methodintervals-（i-1）xaxislabelsval_methodsegments）;
    [temp1temp2minRMSEwithout（i）ix_for_iterative（i）minRMSEglobal（i）]=sort_ipls（ModelReverse）;  % Subfun

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         122  2003-11-16 21:29  makeManualIntervals.m
     文件         308  2004-07-16 08:03  makeManualSegments.m
     文件          29  2000-03-16 14:33  makeSampleNames.m
     文件        3241  2005-03-22 11:46  mwpls.m
     文件        1169  2004-07-20 20:58  mwplsdemo.m
     文件        3547  2005-03-21 16:46  mwplsplot.m
     文件        3118  2005-03-21 16:46  mwplsplot_comp.m
     文件        3985  2005-03-22 11:45  plsmodel.m
     文件        4449  2005-03-22 13:05  plspredict.m
     文件        2789  2004-07-19 09:31  plspress.m
     文件       14372  2005-03-22 11:00  plspvsm.m
     文件        4879  2004-07-20 21:03  plsrmse.m
     文件        6442  2005-03-22 11:39  sipls.m
     文件        2696  2004-07-23 11:43  siplsdemo.m
     文件        1559  2004-07-20 21:55  siplstable.m
     文件        4457  2005-03-22 13:39  sub_bipls_limit.m
     文件        1409  2005-03-22 11:51  sub_bipls_vector.m
     文件        1674  2005-03-22 13:51  sub_bipls_vector_limit.m
     文件        3698  2005-03-22 11:48  sub_iplsreverse.m
     文件        2016  2005-03-22 11:29  sub_pca.m
     文件        1842  2005-03-22 11:36  sub_pls.m
     文件         549  2005-03-22 12:00  sub_pls_pre.m
     文件        8939  2005-03-22 11:38  sub_pls_val.m
     文件      544686  2005-03-22 07:52  iToolbox_Manual.pdf
     文件      447248  2004-06-27 13:51  nirbeer.mat
     文件        6216  2005-03-22 11:54  bipls.m
     文件        3916  2005-03-22 13:35  bipls_dyn.m
     文件        2232  2004-07-14 08:43  biplsdemo.m
     文件         949  2004-07-14 08:43  biplstable.m
     文件        2799  2005-03-22 11:03  contents.m
     文件         862  2004-07-11 13:14  intervals.m
............此处省略11个文件信息