日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當(dāng)前位置: 首頁(yè) > 编程资源 > 编程问答 >内容正文

编程问答

R语言实现描述性统计

發(fā)布時(shí)間:2024/10/8 编程问答 31 豆豆
生活随笔 收集整理的這篇文章主要介紹了 R语言实现描述性统计 小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
# 01分布 a<- runif(20) a
  • 0.0534173720516264
  • 0.0381318787112832
  • 0.253885793266818
  • 0.851638266118243
  • 0.356620342005044
  • 0.175922254333273
  • 0.270358079113066
  • 0.421792675741017
  • 0.675487545551732
  • 0.139561568852514
  • 0.649348761420697
  • 0.0383495420683175
  • 0.673801982775331
  • 0.131142142694443
  • 0.241756724659353
  • 0.205821343231946
  • 0.826634412631392
  • 0.827650502324104
  • 0.48426380334422
  • 0.385196640854701
  • # 算術(shù)平均數(shù) mean(a)

    0.385039081587456

    # 幾何平均數(shù) exp(mean(log(a)))

    0.269715541826603

    # 中位數(shù) median(a)

    0.313489210559055

    b <- sort(a) (b[10]+b[11])/2

    0.313489210559055

    # 產(chǎn)生0到10的20個(gè)數(shù)round取整 c <- round(runif(20,0,10)) c
  • 6
  • 1
  • 10
  • 9
  • 6
  • 9
  • 3
  • 10
  • 6
  • 3
  • 4
  • 4
  • 7
  • 2
  • 3
  • 3
  • 5
  • 9
  • 2
  • 0
  • # 眾數(shù) table 統(tǒng)計(jì)出現(xiàn)的次數(shù) x <- table(c) x c0 1 2 3 4 5 6 7 9 10 1 1 2 4 2 1 3 1 3 2 # 取出次數(shù)最多的 names(x)[x==max(x)]

    ‘3’

    a <- round(runif(100,1,10)) table(a) a1 2 3 4 5 6 7 8 9 10 4 8 11 7 16 10 11 8 14 11 stem(a) The decimal point is at the |1 | 00002 | 000000003 | 000000000004 | 00000005 | 00000000000000006 | 00000000007 | 000000000008 | 000000009 | 0000000000000010 | 00000000000 hist(a)

    hist(a,breaks=50)

    a<- rnorm(100,0,1) hist(a)

    # 協(xié)方差與相關(guān)系數(shù) x<-runif(20) y<-runif(20) cov(x,y)

    0.00322176575114672

    e<- runif(20,0,0.1) z<- x*3+e cov(x,z)

    0.191163714793349

    #與相關(guān)系數(shù) cor(x,z)

    0.999193595914862

    # 一元線性回歸| x<- runif(10) y<- 3*x +5+ runif(10,0,0.5) plot(x,y)

    # 線性模型 lm(y~x) Call: lm(formula = y ~ x)Coefficients: (Intercept) x 5.349 2.780 plot(x,y) abline(lm(y~x),col="red")

    # 多元回歸 x1<- runif(10) x2 <- runif(10) y <- 3*x1+5*x2+2+runif(10,0,0.1) lc<- lm(y~x1+x2) summary(lc) Call: lm(formula = y ~ x1 + x2)Residuals:Min 1Q Median 3Q Max -0.04436 -0.02495 0.00357 0.02298 0.04432 Coefficients:Estimate Std. Error t value Pr(>|t|) (Intercept) 2.08064 0.03958 52.56 2.36e-10 *** x1 2.97190 0.03884 76.51 1.71e-11 *** x2 4.97719 0.04625 107.61 1.58e-12 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1Residual standard error: 0.03303 on 7 degrees of freedom Multiple R-squared: 0.9995, Adjusted R-squared: 0.9993 F-statistic: 6916 on 2 and 7 DF, p-value: 2.91e-12

    總結(jié)

    以上是生活随笔為你收集整理的R语言实现描述性统计的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問題。

    如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò),歡迎將生活随笔推薦給好友。