直线矩阵数学之路(2)-四大神器-R(18)
最近研究直線矩陣，稍微總結(jié)一下，以后繼續(xù)補(bǔ)充：
????目標(biāo)是這些點(diǎn)到這條直線的距離的平方和最小，可運(yùn)用最小二乘法，最小二乘法擬合的進(jìn)程就是回歸，這條直線就是回歸線。
????Lsfit()函數(shù)實(shí)現(xiàn)最小二乘法擬合，其主要參數(shù)為：
????X：一個(gè)矩陣的行對應(yīng)的情況和其列對應(yīng)為變量。
????Y：結(jié)果，可所以一個(gè)矩陣，如果你想，以適應(yīng)多種左手側(cè)。
????Wt:可選參數(shù)，加權(quán)最小二乘法的執(zhí)行權(quán)重向量。
????Intercept:是否應(yīng)使用截距項(xiàng)。
????Tolerance:公差將用于在矩陣分解
????Yname:用于響應(yīng)變量的名稱。
????我們以x=(1,2,3,4),y=(2,4,6,8)，可得到回歸線方程為
????Y=2x
????>?y<-c(2,4,6,8)
????>?x<-c(1,2,3,4)
????>?lsfit(x,y)
????$coefficients
????Intercept?????????X?
????????0?????????2?
????........
????........
????上述結(jié)果中,Intercept項(xiàng)表現(xiàn)截距，x項(xiàng)表現(xiàn)方程的x的常數(shù)項(xiàng)。
????我們先假設(shè)回歸線為
????Y=2x+3
????然后，根據(jù)回歸線結(jié)構(gòu)x和y值。
????>?y<-c(5,7,9,11)
????>?x<-c(1,2,3,4)
????執(zhí)行l(wèi)sfit()函數(shù)
????>?lsfit(x,y)
????$coefficients
????Intercept?????????X?
????????3?????????2?
????要正確得出方程的截距為3，x的常數(shù)項(xiàng)為2。現(xiàn)實(shí)生活中，很難有如此精確的模型，我們再多結(jié)構(gòu)一些點(diǎn)：
????>?y<-c(5,7,9,11,16,20)
????>?x<-c(1,2,3,4,7,9)
????>?lsfit(x,y)
????>?x<-c(1,2,3,4,7,9)
????>?y<-c(5,7,9,11,16,20)
????我們通過plot(x,y)來繪制這些點(diǎn)在直角坐標(biāo)系中的位置，這個(gè)圖也被稱為散點(diǎn)圖。
????>?plot(x,y)
????
????>?lsfit(x,y)
????$coefficients
????Intercept?????????X?
?3.338028??1.845070?
????$residuals
每日一道理只有啟程，才會(huì)到達(dá)理想和目的地，只有拼搏，才會(huì)獲得輝煌的成功，只有播種，才會(huì)有收獲。只有追求，才會(huì)品味堂堂正正的人。
????[1]?-0.18309859?-0.02816901??0.12676056??0.28169014?-0.25352113??0.05633803
????Coefficients為系數(shù)，包含截距和x的系數(shù)，residuals表現(xiàn)殘差，殘差分別反響了這些點(diǎn)與直線的差異，殘差越小越好，我們將回歸線也畫上
????>?abline(lsfit(x,y))
????
????可以看到擬合效果還是不錯(cuò)的，我們也可以使用lm（）函數(shù)，來建立線性模型停止回歸分析：
????畫x,y的散點(diǎn)圖：?plot(x,y)
????做相關(guān)回歸分析，結(jié)果存放在xy中：?lm(y~x)->xy
????顯示xy的相關(guān)回歸分析結(jié)果:summary(xy)
????畫回歸線:>??abline(?lm(y~x))
????
本博客全部內(nèi)容是原創(chuàng)，未經(jīng)書面許可，嚴(yán)禁任何情勢的轉(zhuǎn)載。
????
http://blog.csdn.net/u010255642
文章結(jié)束給大家分享下程序員的一些笑話語錄： 打賭
飛機(jī)上，一位工程師和一位程序員坐在一起。程序員問工程師是否樂意和他一起玩一種有趣的游戲。工程師想睡覺，于是他很有禮貌地拒絕了，轉(zhuǎn)身要睡覺。程序員堅(jiān)持要玩并解釋說這是一個(gè)非常有趣的游戲："我問你一個(gè)問題，如果你不知道答案，我付你5美元。然后你問我一個(gè)問題，如果我答不上來，我付你5美元。"然而，工程師又很有禮貌地拒絕了，又要去睡覺。 程序員這時(shí)有些著急了，他說："好吧，如果你不知道答案，你付5美元；如果我不知道答案，我付50美元。"果然，這的確起了作用，工程師答應(yīng)了。程序員就問："從地球到月球有多遠(yuǎn)？"工程師一句話也沒有說，給了程序員5美元。 現(xiàn)在輪到工程師了，他問程序員："什么上山時(shí)有三條腿，下山卻有四條腿？"程序員很吃驚地看著工程師，拿出他的便攜式電腦，查找里面的資料，過了半個(gè)小時(shí)，他叫醒工程師并給了工程師50美元。工程師很禮貌地接過錢又要去睡覺。程序員有些惱怒，問："那么答案是什么呢？"工程師什么也沒有說，掏出錢包，拿出5美元給程序員，轉(zhuǎn)身就去睡覺了。
轉(zhuǎn)載于:https://www.cnblogs.com/jiangu66/archive/2013/05/12/3074153.html
總結(jié)
以上是生活随笔為你收集整理的直线矩阵数学之路(2)-四大神器-R(18)的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。
- 上一篇： 解决GIT代码仓库不同步问题
- 下一篇： jQuery获取CSS样式中的颜色值的问