其中包括Hmisc、pastecs和psych。由于這些包并未包括在基礎(chǔ)安裝中。
創(chuàng)新互聯(lián)-專業(yè)網(wǎng)站定制、快速模板網(wǎng)站建設(shè)、高性價(jià)比明溪網(wǎng)站開發(fā)、企業(yè)建站全套包干低至880元,成熟完善的模板庫,直接使用。一站式明溪網(wǎng)站制作公司更省心,省錢,快速模板網(wǎng)站建設(shè)找我們,業(yè)務(wù)覆蓋明溪地區(qū)。費(fèi)用合理售后完善,十載實(shí)體公司更值得信賴。
> library(Hmisc)
> describe(mtcars[vars])
mtcars[vars]
3 Variables 32 Observations
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------
mpg
n missing unique Info Mean .05 .10 .25 .50 .75 .90 .95
32 0 25 1 20.09 12.00 14.34 15.43 19.20 22.80 30.09 31.30
lowest : 10.4 13.3 14.3 14.7 15.0, highest: 26.0 27.3 30.4 32.4 33.9
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------
hp
n missing unique Info Mean .05 .10 .25 .50 .75 .90 .95
32 0 22 1 146.7 63.65 66.00 96.50 123.00 180.00 243.50 253.55
lowest : 52 62 65 66 91, highest: 215 230 245 264 335
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------
wt
n missing unique Info Mean .05 .10 .25 .50 .75 .90 .95
32 0 29 1 3.217 1.736 1.956 2.581 3.325 3.610 4.048 5.293
lowest : 1.513 1.615 1.835 1.935 2.140, highest: 3.845 4.070 5.250 5.345 5.424
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------
>
--------------------------------------------------------------------
包中有一個(gè)名為stat.desc()的函數(shù),它可以計(jì)算種類繁多的描述性統(tǒng)計(jì)量。使用
格式為:
> stat.desc(mtcars[vars],basic=TRUE,desc=TRUE,norm=FALSE,p=0.95)
其中的x是一個(gè)數(shù)據(jù)框或時(shí)間序列。若basic=TRUE(默認(rèn)值),則計(jì)算其中所有值、空值、缺失
值的數(shù)量,以及最小值、最大值、值域,還有總和。若desc=TRUE(同樣也是默認(rèn)值),則計(jì)算
中位數(shù)、平均數(shù)、平均數(shù)的標(biāo)準(zhǔn)誤、平均數(shù)置信度為95%的置信區(qū)間、方差、標(biāo)準(zhǔn)差以及變異系
數(shù)。最后,若norm=TRUE(不是默認(rèn)的),則返回正態(tài)分布統(tǒng)計(jì)量,包括偏度和峰度(以及它們
的統(tǒng)計(jì)顯著程度)和Shapiro–Wilk正態(tài)檢驗(yàn)結(jié)果。這里使用了p值來計(jì)算平均數(shù)的置信區(qū)間(默
認(rèn)置信度為0.95)
> library(pastecs)
> stat.desc(mtcars[vars])
mpg hp wt
nbr.val 32.0000000 32.0000000 32.0000000
nbr.null 0.0000000 0.0000000 0.0000000
nbr.na 0.0000000 0.0000000 0.0000000
min 10.4000000 52.0000000 1.5130000
max 33.9000000 335.0000000 5.4240000
range 23.5000000 283.0000000 3.9110000
sum 642.9000000 4694.0000000 102.9520000
median 19.2000000 123.0000000 3.3250000
mean 20.0906250 146.6875000 3.2172500
SE.mean 1.0654240 12.1203173 0.1729685
CI.mean.0.95 2.1729465 24.7195501 0.3527715
var 36.3241028 4700.8669355 0.9573790
std.dev 6.0269481 68.5628685 0.9784574
coef.var 0.2999881 0.4674077 0.3041285
------------------------------------------------
library(psych)
describe(mtcars[vars])
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
在前面的示例中,psych包和Hmisc包均提供了名為describe()的函數(shù)。R如何知道該
使用哪個(gè)呢?簡言之,如代碼清單7-5所示,最后載入的程序包優(yōu)先。在這里,psych在
Hmisc之后被載入,然后顯示了一條信息,提示Hmisc包中的describe()函數(shù)被psych
包中的同名函數(shù)所屏蔽(masked)。鍵入describe()后,R在搜索這個(gè)函數(shù)時(shí)將首先找
到psych包中的函數(shù)并執(zhí)行它。如果你想改而使用Hmisc包中的版本,可以鍵入
Hmisc::describe(mt)。這個(gè)函數(shù)仍然在那里。你只是需要給予R更多信息以找到它
標(biāo)題名稱:R語言筆記計(jì)算描述性統(tǒng)計(jì)量的函數(shù)
瀏覽地址:http://jinyejixie.com/article14/possde.html
成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供外貿(mào)建站、App設(shè)計(jì)、關(guān)鍵詞優(yōu)化、網(wǎng)站導(dǎo)航、、品牌網(wǎng)站設(shè)計(jì)
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請盡快告知,我們將會在第一時(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)