matlab線性回歸程式

2021-06-21 04:27:32 字數 564 閱讀 6169

最近用matlab做線性回歸,用到命令regress,可是發現他沒辦法做回歸係數的t檢驗,因此,自己就寫了乙個。參考書籍:《計量經濟學》龐皓 科學出版社

function [b,beta,t,t0,sigma,rss,tss,r,r1,f,f0]=pureregress(x,y,con)

%x是因子,n*p,不包括場數列,y是被解釋變數,n*1.con是置信度,預設:0.95,一列代表乙個因子

%[n,k]=size(x);

eyes=eye(k+1,k+1);%對角單位矩陣

x=[ones(n,1),x];%構建結構陣x,

a=x'*x; %求算資訊陣a,

c=eyes/a; %求算資訊陣的逆陣

b=x\y; % 求算回歸統計數向量,其中第一行為回歸截距

y_1=x*b;%y的估計值序列

rss=sum((y-y_1).^2);%殘差平方和

tss=sum((y-mean(y)).^2);%總變差

ess=tss-rss;%回歸變差 

r=1-(rss/tss);%未修正的r方或者說可決係數

r1=1-

用python編寫線性回歸程式

在網上本來找到了程式模板以為新建python專案然後導過去就能用了,可是我還是想的太簡單,導過去後程式有很多庫並沒有,就像這樣。本來還有sklearn庫,不過我用pip install sklearn 安裝了,如果你沒有pip,請先安裝pip。然後pandas也是要安裝的,也是命令列直接pip in...

MATLAB 線性回歸

一 一元線性回歸 1 1 命令 polyfit最小二乘多項式擬合 p,s polyfit x,y,m 多項式y a1xm a2xm 1 amx am 1 其中x x1,x2,xm x1 xm為 n 1 的矩陣 y為 n 1 的矩陣 p a1,a2,am 1 是多項式y a1xm a2xm 1 amx...

線性回歸(Matlab)

單變數線性回歸 假設方程 構造矩陣x時,使x 1 1 代價函式 y是訓練集中實際值,h x 是 值 利用gradient descent 梯度下降 一步步迭代求出使代價函式最小的 注意 對於單變數,是二維的 data load ex1data1.txt x data 1 y data 2 m len...