24小时热门版块排行榜

>论坛更新日志 (3856)
>文献求助 (403)
>虫友互识 (328)
>导师招生 (313)
>博后之家 (155)
>硕博家园 (154)
>休闲灌水 (138)
>考博 (124)
>论文投稿 (110)
>基金申请 (90)
>招聘信息布告栏 (54)
>考研 (54)
>教师之家 (48)
>公派出国 (43)
>绿色求助(高悬赏) (37)
>找工作 (37)

返回列表

赵环

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 6
帖子: 1
在线: 22分钟
虫号: 3865278
注册: 2015-05-12
性别: GG
专业: 数理统计

[交流] 谁知道这里面的data和newdata指的是什么。。谢谢了

function [cdata,me,ctest]=center(data,opt,newdata);
%#
%#  AIM: Centering along columns, rows or double centering
%#
%#  PRINCIPLE:  Removal of the column-, row- or overall mean from
%#             each column, row or both, respectively
%#    If a test data set is available it can ONLY be
%#             column centered using the mean from the calibration
%#             data set.
%#
%#
%#  INPUT: data: (m x n) matrix with m rows and n variables
%# opt: optional
%#       1 = column centering
%#       2 = row centering
%#       3 = double centering
%#       newdata: (mt x n) test matrix with mt rows and n variables
%#
%#
%#  OUTPUT: cdata: (m x n) matrix containing centered data
%# me: mean vector, overall mean (scalar)
%#             newdata: (mt*n) test matrix centered with the mean of data
%#
%#
%#  AUTHOR: Andrea Candolfi
%#         Copyright(c) 1997 for ChemoAc
%#       FABI, Vrije Universiteit Brussel
%#       Laarbeeklaan 103 1090 Jette
%#
%# VERSION: 1.2 (25/02/2002)
%#
%#  TEST:    I. Stanimirova & S. Gourv閚ec & M. Zhang
%#

function [cdata,me,cnewdata]=center(data,opt,newdata);

[m,n]=size(data);

if nargin==1;
  opt=[4];
  while opt>3 | opt<=0
opt=input('column centering(1), row centering(2), double centering(3):');
  end
end

if opt==1 % column centering
me=mean(data);
cdata=data-ones(m,1)*me;
end

if opt==2 % row centering
me=mean(data')';
cdata=data-me*ones(1,n);
end

if opt==3 % double centering
me=mean(mean(data));
mej=mean(data');
mei=mean(data);
cdata=data-(ones(m,1)*mei)-(ones(n,1)*mej)'+(ones(m,n)*me);
end

if exist('newdata')==1 % center new data
[mt,n]=size(newdata);

if opt==1 % column centering
      me=mean(data);
      cnewdata=newdata-ones(mt,1)*me;
else
      error('Row centering and double centering are impossible to perform on a test set');
end

end

回复此楼