boxc_trans.m

计量工具箱

function bdata = boxc_trans(x,lam)
% PURPOSE: compute box-cox transformation
%----------------------------------------------------
% USAGE: bdata = boxc_trans(data,lam)
% where:    lam  = scalar transformation parameter
%           data = matrix nobs x k
%----------------------------------------------------
% RETURNS: bdata = data matrix box-cox transformed

% written by:
% James P. LeSage, Dept of Economics
% University of Toledo
% 2801 W. Bancroft St,
% Toledo, OH 43606
% [e-mail]
 
[n k] = size(x);
z = zeros(n,k);
iota = ones(n,1);

 for i=1:k;
   if lam ~= 0
   z(:,i) = (x(:,i).^lam - iota)/lam;
   else
   z(:,i) = log(abs(x(:,i)));
   end;
 end;

bdata = z;