您好,登錄后才能下訂單哦!
在比較多組個(gè)體或觀測(cè)時(shí),關(guān)注的焦點(diǎn)經(jīng)常是各組的描述性統(tǒng)計(jì)信息,而不是樣本整體的描
述性統(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
免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng),如果涉及侵權(quán)請(qǐng)聯(lián)系站長郵箱:is@yisu.com進(jìn)行舉報(bào),并提供相關(guān)證據(jù),一經(jīng)查實(shí),將立刻刪除涉嫌侵權(quán)內(nèi)容。