该命令x=[ones(size ) T1 ]、t1、T2、T3、T4 )； [b，bint，r，rint，stats]=regress(y，x )，Q=r'*r，Q2=q/) 15-2 )，stepwise(X ) x，y，[ 1，2，3 ]
ifnargin 2
错误(stats :注册： toofewinputs )、
' regressrequiresatleasttwoinputarguments.'；
elseif nargin==2
阿尔法=0.05；
结束
if~isvector(y )
error(stats:Regress:追寻的小蜜蜂) Y must be a vector.)；
elseifnumel(y ) (=n
error(Stats:Regress:追寻的小蜜蜂) .
' thenumberofrowsinymustequalthenumberofrowsinx.'；
结束
Havenans=any(wasNan )；
if havenans
y(wasNan )=[]；
x(wasnan， )=[]；
n=length(y；
结束
p=sum(ABS(diag(r ) ) max (n，ncolX ) ) EPS ) r )1)；
if p ncolX
警告(stats : regress : rankdefdesignmat )、
' xisrankdeficienttowithinmachineprecision.'；
r=r(1:p、1:p )；
q=q65:1:p；
perm=perm(1:p；
结束
b=Zeros(ncolx，1 )；
b(perm )=r () q'*y )；
if nargout=2
nu=max(0，n-p )； % Residual degrees of freedom
yhat=X*b； % predictedresponsesateachdatapoint。
r=y-yhat； % Residuals。
jqdbgr=jqdbg(r；
if nu ~=0
RMSE=jqdbgr/sqrt(nu； % Root mean square error。
tval=tinv ((1- alpha/2 )，nu )；
else
rmse=NaN；
tval=0；
结束
if nargout=4
hatdiag=sum(ABS(q ).^ 2，2 )；
ok=() (1-hatdiag ) sqrt (EPS ) class (hat diag ) )
hatdiag(~ok )=1；
国际空间站1
Enom=(nu-1 ).* ) * (1-hatdiag )；
sigmai=Zeros(Length(Denom )，1 )；
sigmai(ok )=sqrt (max ) 0，) nu * S2/(nu-1 ) )-(r ) ok ).^2./denom (ok ) ) )
elseif nu==1
else % if nu==0
结束
结束
if nargout==5
r2=1 - SSE/TSS； % r-square静态。
国际足球锦标赛
f=(RSS/) p-1 ) )/s2； % F statistic for regression
else
F=NaN；
结束
prob=1-fcdf(f，p-1，nu )； % significanceprobabilityforregression
stats=[r2 F prob s2]；
if至any (all (x==1，1 ) )
B0=R () Q'*Ones(N，1 )；
if ) sum(ABS(1-x ) :perm (B0 ) ) ) n*sqrt ) EPS ) class ) x ) )
warning (stats : regress :否const )、
' unlessxhasacolumnofones. ntype ' help ' .
' regress' for more information.']；
结束
结束
结束
if havenans
if nargout=3
TMP=repmat(Nan，Length ) wasNan )，1 )；
tmp(wasnan )=r；
r=tmp；
if nargout=4
TMP=repmat(Nan，Length ) wasNan )，2 )；
TMP(wasnan， )=rint；
结束
结束
结束
end % nargout=2
错误由警告和错误提示。 请高手帮忙解决，本人的毕业设计很着急！ 谢谢你！