基于Matlab的线性回归函数-经管之家官网！

%自变量矩阵：X

%因变量矩阵：Y

%线性回归系数：RegCoff

%可决系数：R2

%修正的可决系数：R21

%复相关系数：R

%方程显著性检验F值：F

%各因子显著性检验F值：FX

%各因子显著性检验t值：tX

format long;

n = size(X,1);

k = size(X,2);

RegCoff = zeros(k+1,1);%回归系数

Z = mean(X);

yp = mean(Y);

A = transpose(X)*X-n*transpose(Z)*Z;

C = transpose(X)*Y-n*transpose(Z)*yp;

RegCoff(2:k+1) = A\C;

RegCoff(1) = yp-Z*RegCoff(2:k+1);%回归系数的常数项

%显著性检验

TSS = norm(Y)^2-n*yp^2;%总离差平方和 TSS

YR = X*RegCoff(2:k+1) + RegCoff(1)*ones(n,1);

ESS = transpose(RegCoff(2:k+1))*C;%回归平方和 ESS

RSS = TSS-ESS;%剩余平方和 RSS

R2 = ESS/TSS;

R21 = 1-(1-R2)*((n-1)/(n-k-1));

R0 = sqrt(R2);

UR = ESS/(length(RegCoff)-1);

QR = RSS/(n-length(RegCoff));

s = sqrt(QR);

inA = inv(A);

F = UR/QR;

for i=1:length(RegCoff)-1%各因子显著性检验

FX(i)=RegCoff(i+1)^2/inA(i,i)/QR;

tX(i)=RegCoff(i+1)/sqrt(inA(i,i))/s;

end

format short;

%%%%%%%%%%

FX(i)、tX(i)只能检验各因子显著性，本人还想编出回归方程常数项显著性的F检验、t检验的代码，但是尝试了没有成功。

同时用regress也没有试成功。

在这里向大神求助，先谢谢大神！拜谢！！！

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP