c語言版的線性回歸分析函數(shù)(存儲版)
【正文】 n+1,+(bdatamcount。(a0。p(p[n]+cols]}Answer)。//0{bY的平均值SquarePoor[0](mm(i(*p*=if (rows cols )n。=981,1386,1672,1864,2143,2176,1769,2207,2251,2390,2287,2388,2093,2083,2390,doubleinti)。printf(剩余平方和：%dat。pforv*pdouble(MultipleRegression((double*)data,SquarePoor,0。什么？你問我懂不懂？呵呵，不瞞你說，我的主業(yè)就是統(tǒng)計，而且統(tǒng)計師職務已經有20年，也就是文革 后第一批評定的，而且第一批全國自學考試統(tǒng)計大專畢業(yè)，編程序只是我的業(yè)余愛好，不過我退養(yǎng)休息了近10年，也忘得差不多了，但是還是能看懂這些簡單的東 東。其運行結果是一樣的，可能以前我為了DOS下的運行速度，單獨寫了一個函數(shù)，因為畢竟多元回歸分析很少用到，而一元回歸是經常使用的。其實，一元線性回歸是多元線性回歸的一個特例，完全可以使用同一個函數(shù)，如前面的例子：if (MultipleRegression((double*)data1, 12, 2, Answer, SquarePoor) == 0)Display((double*)data,0)}**p,=ipSquarePoor[3],SquarePoor[1],%.5lf*X%d,double1613,1612,1400,1920,1328,2199,2016,1955,1873,2050,1890,1970,1565,1771,1536,X4下面是一個多元回歸的例子：double}else回歸方差////*a++)rows。data。b(nLinearEquations(dat,//*ammn*p。rows。p*++acols)m,i*ndat[(nSum(Xn)(cols//=bn,=to0。(double)rows。*||||n,*Answer,返回值：0求解成功，1錯誤int...BnXn//求多元回歸方程：Ycount)。=1]。dat[count1的元=0。dat[n][j]++)=/nmmemcpy(dat[n],sizeof(double))。amp。m++)sizeof(double))。*for=*dn。double解線性方程。i(i同一元線性回歸相比，多元線性回歸分析代碼可就復雜多了，必須求解線性方程，因此本代碼中包含一個可獨立使用的線性方程求解函數(shù)：0。SquarePoor,amp。SquarePoor[4]。system(pause)。v,*p++,=iforprintf(剩余分析: )。/驗：%printf(離差平方和：%回歸方差：% ,cols。i,double*SquarePoor,{,},{,回歸平方和：剩余平方和：實例計算：double}///Lxy))。SquarePoor[1]U(回歸平方和)+p0。0。b*XayaLxy=//ya))。((*p*(p++。mmb(data,Lxxint*a,SquarePoor[4]：返回方差分析指標:求線性回歸方程：Y+LinearRegression(doubledoubleLxy||||*p/=Y平均值mm{Lxxxa))。=*Sum((X}
數(shù)學相關推薦
鄂ICP備17016276號-1