R語言筆記分組計(jì)算描述性統(tǒng)計(jì)量aggregate,summaryBy,describe.by

在比較多組個(gè)體或觀測時(shí)，關(guān)注的焦點(diǎn)經(jīng)常是各組的描述性統(tǒng)計(jì)信息，而不是樣本整體的描

創(chuàng)新互聯(lián)專注為客戶提供全方位的互聯(lián)網(wǎng)綜合服務(wù)，包含不限于成都網(wǎng)站建設(shè)、做網(wǎng)站、岫巖網(wǎng)絡(luò)推廣、微信小程序定制開發(fā)、岫巖網(wǎng)絡(luò)營銷、岫巖企業(yè)策劃、岫巖品牌公關(guān)、搜索引擎seo、人物專訪、企業(yè)宣傳片、企業(yè)代運(yùn)營等，從售前售中售后，我們都將竭誠為您服務(wù)，您的肯定，是我們最大的嘉獎(jiǎng)；創(chuàng)新互聯(lián)為所有大學(xué)生創(chuàng)業(yè)者提供岫巖建站搭建服務(wù)，24小時(shí)服務(wù)熱線：028-86922220，官方網(wǎng)址：www.muchs.cn

述性統(tǒng)計(jì)信息。同樣地，在R中完成這個(gè)任務(wù)有若干種方法。我們將以獲取變速箱類型各水平的

描述性統(tǒng)計(jì)量開始。

vars<- c("mpg","hp","wt")

> aggregate(mtcars[vars],by=list(am=mtcars$am),median)

  am  mpg  hp   wt

1  0 17.3 175 3.52

2  1 22.8 109 2.32

注意list(am=mtcars$am)的使用。如果使用的是list(mtcars$am)，則am列將被標(biāo)注為

Group.1而不是am。你使用這個(gè)賦值指定了一個(gè)更有幫助的列標(biāo)簽。如果有多個(gè)分組變量，可以

使用by=list(name1=groupvar1, name2=groupvar2, ... , groupvarN)這樣的語句。

> aggregate(mtcars[vars],by=list(am=mtcars$am,gear=mtcars$gear),median)

  am gear   mpg    hp    wt

1  0    3 15.50 180.0 3.730

2  0    4 21.00 109.0 3.315

3  1    4 25.05  79.5 2.260

4  1    5 19.70 175.0 2.770

------------------------------------------------

doBy包和psych包也提供了分組計(jì)算描述性統(tǒng)計(jì)量的函數(shù)。同樣地，它們未隨基本安裝發(fā)布，

必須在首次使用前進(jìn)行安裝。doBy包中summaryBy()函數(shù)的使用格式為

> summaryBy(mpg+hp+wt~am,data=mtcars,FUN=mtstat)

  am mpg.n mpg.mean mpg.stdev mpg.skew mpg.kurtosis hp.n  hp.mean hp.stdev   hp.skew hp.kurtosis wt.n wt.mean  wt.stdev   wt.skew wt.kurtosis

1  0    32 20.09062  6.026948 0.610655    -0.372766   32 146.6875 68.56287 0.7260237  -0.1355511   32 3.21725 0.9784574 0.4231465 -0.02271075

>

> describe.by(mtcars[vars],mtcars$am)

group: 0

    vars  n   mean    sd median trimmed   mad   min    max  range skew kurtosis    se

mpg    1 32  20.09  6.03  19.20   19.70  5.41 10.40  33.90  23.50 0.61    -0.37  1.07

hp     2 32 146.69 68.56 123.00  141.19 77.10 52.00 335.00 283.00 0.73    -0.14 12.12

wt     3 32   3.22  0.98   3.33    3.15  0.77  1.51   5.42   3.91 0.42    -0.02  0.17

------------------------------------------------------------------------------

的reshape包靈活地按組導(dǎo)出描述性統(tǒng)計(jì)量

> dfs <- melt(mtcars,measure.vars=c("mpg","hp","wt"),id.vars=c("am","cyl"))

> cast(dfs,am+cyl+variable~.,dstats)

  am cyl variable  n       mean         sd

1  0   4      mpg 11  26.663636  4.5098277

2  0   4       hp 11  82.636364 20.9345300

3  0   4       wt 11   2.285727  0.5695637

4  0   6      mpg  7  19.742857  1.4535670

5  0   6       hp  7 122.285714 24.2604911

6  0   6       wt  7   3.117143  0.3563455

7  0   8      mpg 14  15.100000  2.5600481

8  0   8       hp 14 209.214286 50.9768855

9  0   8       wt 14   3.×××14  0.7594047

網(wǎng)頁標(biāo)題：R語言筆記分組計(jì)算描述性統(tǒng)計(jì)量aggregate,summaryBy,describe.by
URL鏈接：http://www.muchs.cn/article2/ishgoc.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián)，為您提供網(wǎng)站制作、企業(yè)建站、網(wǎng)站策劃、關(guān)鍵詞優(yōu)化、外貿(mào)網(wǎng)站建設(shè)、搜索引擎優(yōu)化

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主，如果涉及侵權(quán)請盡快告知，我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場，如需處理請聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載，或轉(zhuǎn)載時(shí)需注明來源： 創(chuàng)新互聯(lián)