再窺R(包括一個R示例)
一、將輸出作為輸入——結果的重用
R有個特點,分析的結果可以保存下來,並可以作為進一步分析的輸入使用。下面我們通過R中的預先安裝好的數據集作為實例(這裏只講述原理,涉及到的統計知識,以後在另外的統計專欄中描述):
1、利用汽車數據mtcars執行一次簡單的線性回歸,通過車身重量(wt)預測每加侖行駛的英裏數(mpg),語句如下:
lm(mpg~wt,data= mtcars)
運行結果(結果只是展示在屏幕上,並沒有保存下來):
2、再次執行回歸。這次將其結果保存起來,語句如下:
myfit <- lm(mpg~wt,data=mtcars)
上述語句創建了一個名為myfit的列表對象,其包含了預測值、殘差、回歸系數等信息。雖然我們在屏幕上沒有看到任何信息,但分析結果已經被保存下來,我們就可以在接下來的分析中使用這些分析信息。
3、執行summary(myfit) ,顯示分析結果的統計概要。效果如下:
執行plot(myfit) ,生成回歸診斷圖。 根據提示,按下enter鍵,最後效果如下:
運行語句:
cook <-cooks.distance(myfit)
plot(cook)
第一句首先計算影響度量統計量,第二句是對結果繪圖。效果如下:
執行語句:predict(myfit,mydata)
可以根據你給定的新的車身重量對每加侖行駛的英裏數進行預測。
註:可以根據前文提到的help(對象名)或者 ?對象名 方法來查閱對象的幫助文檔。其中,value部分對應著函數的返回值。
二、R實踐示例
任務描述:
①打開幫助文檔首頁,查閱其中的“Introduction to R”;
②安裝vcd包(一個用於可視化類別數據的包)
③列出次包中可用的函數和數據集
④載入這個包,並閱讀數據集Arthritis的描述
⑤顯示數據集Arthritis的內容;
⑥運行數據集Arthritis自帶的示例。(基本顯示接收治療的關節炎患者較接受安慰劑的患者在病情上有了更多改善)
最後退出。
代碼清單和演示效果如果:
help.start()
install.packages("vcd")
help(package = "vcd")
library(grid)
library(vcd)
help(Arthritis)
效果:
小結
加上上篇介紹,我們從R的安裝開始,討論了如果通過下載附加包來增強R功能;探索了R的基本界面,以交互和批處理的兩種方式運行R程序;並嘗試繪制一些圖形;學習了如何將工作保存到文本和圖形文件中;並花了一些時間來了解如何訪問大量現成的幫助文檔。
接下來的環節,我們將學習R處理的各種數據類型,以及如何從文本文件、其他程序和數據庫系統中導入數據。
再窺R(包括一個R示例)