สถิติเชิงพรรณนาด้วย R

R สำหรับผู้ใช้ SAS

Melinda Higgins, PhD

Research Professor/Senior Biostatistician Emory University

การโหลดชุดข้อมูล CSV จากภายนอก

  • ชุดข้อมูล Abalone มีการวัด 9 ค่า:
    • ความยาว
    • เส้นผ่านศูนย์กลาง
    • ความสูง
    • น้ำหนักทั้งตัว
    • น้ำหนักเนื้อ
    • น้ำหนักเปลือก
    • น้ำหนักอวัยวะภายใน
    • เพศ (วัยอ่อน, เพศเมีย, เพศผู้)
    • จำนวนวงแหวน

ภาพหอยเป๋าฮื้อ

  • สำหรับหอยเป๋าฮื้อ 4,177 ตัว
R สำหรับผู้ใช้ SAS

การโหลดชุดข้อมูล CSV จากภายนอก

  • ชุดข้อมูล abalone มีให้ในรูปแบบ CSV (comma separated value)
  • ใช้ฟังก์ชัน read_csv() จากแพ็กเกจ readr เพื่อโหลดข้อมูล CSV

โลโก้ readr hexsticker

R สำหรับผู้ใช้ SAS

โค้ด SAS proc import สำหรับไฟล์ CSV เทียบกับโค้ด R ที่ใช้ฟังก์ชัน read csv จากแพ็กเกจ readr

R สำหรับผู้ใช้ SAS

SAS proc import เทียบเท่ากับฟังก์ชัน read csv จากแพ็กเกจ readr ของ R

R สำหรับผู้ใช้ SAS

SAS proc import เทียบเท่ากับฟังก์ชัน read csv จากแพ็กเกจ readr ของ R

R สำหรับผู้ใช้ SAS

SAS proc import เทียบเท่ากับฟังก์ชัน read csv จากแพ็กเกจ readr ของ R

  • ตัวดำเนินการกำหนดค่า <- นำผลลัพธ์จาก readr::read_csv ไปเก็บไว้ในอ็อบเจ็กต์ abalone
  • ตอนนี้ abalone ถูกบันทึกไว้ใน global environment แล้ว
R สำหรับผู้ใช้ SAS

SAS proc contents เทียบเท่ากับฟังก์ชัน str structure ของ R

str(abalone)

Classes 'tbl_df', 'tbl' and 'data.frame':    4177 obs. of  9 variables:
 $ sex          : chr  "M" "M" "F" "M" ...
 $ length       : num  0.455 0.35 0.53 0.44 0.33 0.425 0.53 0.545 ...
 $ diameter     : num  0.365 0.265 0.42 0.365 0.255 0.3 0.415 0.425  ...
 $ height       : num  0.095 0.09 0.135 0.125 0.08 0.095 0.15 0.125 ...
 $ wholeWeight  : num  0.514 0.226 0.677 0.516 0.205 ...
 $ shuckedWeight: num  0.2245 0.0995 0.2565 0.2155 0.0895 ...
 $ visceraWeight: num  0.101 0.0485 0.1415 0.114 0.0395 ...
 $ shellWeight  : num  0.15 0.07 0.21 0.155 0.055 0.12 0.33 0.26 ...
 $ rings        : int  15 7 9 10 7 8 20 16 9 19 ...
R สำหรับผู้ใช้ SAS

SAS proc contents เทียบเท่ากับฟังก์ชัน dim และ names ของ R

# แสดงมิติของชุดข้อมูล abalone
dim(abalone)

4177    9

# องค์ประกอบหรือตัวแปรในชุดข้อมูล abalone
names(abalone)

"sex"           "length"        "diameter"     "height"        "wholeWeight"   
"shuckedWeight" "visceraWeight" "shellWeight"  "rings"
R สำหรับผู้ใช้ SAS

เนื้อหาชุดข้อมูลและชนิดตัวแปร

  • head() และ tail() แสดงแถวบนสุดและล่างสุด 6 แถวตามค่าเริ่มต้น
  • เปลี่ยนจำนวนแถวที่แสดงได้โดยเพิ่มอาร์กิวเมนต์ที่สองในฟังก์ชัน
# แสดงแถวล่างสุด 7 แถวของ abalone
tail(abalone, 7)

# A tibble: 7 x 9
  sex   length diameter height wholeWeight shuckedWeight visceraWeight shellWeight rings
  <chr>  <dbl>    <dbl>  <dbl>       <dbl>         <dbl>         <dbl>       <dbl> <dbl>
1 M      0.55     0.43   0.13        0.840         0.316         0.196       0.240    10
2 M      0.56     0.43   0.155       0.868         0.4           0.172       0.229     8
3 F      0.565    0.45   0.165       0.887         0.37          0.239       0.249    11
4 M      0.59     0.44   0.135       0.966         0.439         0.214       0.260    10
5 M      0.6      0.475  0.205       1.18          0.526         0.288       0.308     9
6 F      0.625    0.485  0.15        1.09          0.531         0.261       0.296    10
7 M      0.71     0.555  0.195       1.95          0.946         0.376       0.495    12
R สำหรับผู้ใช้ SAS

การจัดการข้อมูลด้วยแนวทาง dplyr

โลโก้ dplyr hexsticker

ในคอร์สนี้ จะใช้ฟังก์ชัน dplyr ต่อไปนี้:

  • %>% คือตัวดำเนินการ pipe จากแพ็กเกจ magrittr ที่มาพร้อมกับ dplyr
  • arrange() จะเรียงลำดับข้อมูลตามตัวแปรหนึ่งตัวหรือมากกว่า
  • pull(x) จะดึงตัวแปรคอลัมน์ x ออกจากชุดข้อมูล
  • select(x,y,z) จะเลือกตัวแปรมากกว่าหนึ่งตัวจากชุดข้อมูล
R สำหรับผู้ใช้ SAS

ฟังก์ชัน arrange ของ dplyr และแนวทางการใช้ pipe %>%

ฟังก์ชัน arrange ของ R แบบใช้และไม่ใช้ตัวดำเนินการ pipe

R สำหรับผู้ใช้ SAS

ฟังก์ชัน arrange ของ dplyr และแนวทางการใช้ pipe %>%

ฟังก์ชัน arrange ของ R แบบใช้และไม่ใช้ตัวดำเนินการ pipe

R สำหรับผู้ใช้ SAS

ฟังก์ชัน arrange ของ dplyr และแนวทางการใช้ pipe %>%

ฟังก์ชัน arrange ของ R แบบใช้และไม่ใช้ตัวดำเนินการ pipe

R สำหรับผู้ใช้ SAS

ฟังก์ชัน arrange ของ dplyr และแนวทางการใช้ pipe %>%

ฟังก์ชัน arrange ของ R แบบใช้และไม่ใช้ตัวดำเนินการ pipe

R สำหรับผู้ใช้ SAS

เรียงลำดับหอยเป๋าฮื้อตามเส้นผ่านศูนย์กลาง

# เรียงลำดับชุดข้อมูล abalone ตามมิติเส้นผ่านศูนย์กลาง
abalone %>%
  arrange(diameter)

# A tibble: 4,177 x 9
   sex   length diameter height wholeWeight shuckedWeight visceraWeight shellWeight rings
   <chr>  <dbl>    <dbl>  <dbl>       <dbl>         <dbl>         <dbl>       <dbl> <dbl>
 1 I      0.075    0.055  0.01       0.002         0.001         0.0005      0.0015     1
 2 I      0.11     0.09   0.03       0.008         0.0025        0.002       0.003      3
 3 I      0.13     0.095  0.035      0.0105        0.005         0.0065      0.0035     4
 4 I      0.13     0.1    0.03       0.013         0.0045        0.003       0.004      3
 5 I      0.15     0.1    0.025      0.015         0.0045        0.004       0.005      2
 6 I      0.155    0.105  0.05       0.0175        0.005         0.0035      0.005      4
 7 I      0.14     0.105  0.035      0.014         0.0055        0.0025      0.004      3
 8 I      0.17     0.105  0.035      0.034         0.012         0.0085      0.005      4
 9 I      0.14     0.105  0.035      0.0145        0.005         0.0035      0.005      4
10 M      0.155    0.11   0.04       0.0155        0.0065        0.003       0.005      3
R สำหรับผู้ใช้ SAS

ดึงตัวแปรหนึ่งตัวจาก abalone

มาดึงตัวแปร shuckedWeight จาก abalone โดยใช้ pull() จาก dplyr

# ดึงตัวแปร shuckedWeight ออกจาก abalone
abalone %>%
  pull(shuckedWeight)

   [1] 0.2245 0.0995 0.2565 0.2155 0.0895 0.1410 0.2370 0.2940 0.2165 0.3145 0.1940 0.1675
  [13] 0.2175 0.2725 0.1675 0.2580 0.0950 0.1880 0.0970 0.1705 0.0955 0.0800 0.4275 0.3180
  [25] 0.5130 0.3825 0.3945 0.3560 0.3940 0.3930 0.3935 0.6055 0.5515 0.8150 0.6330 0.2270
  [37] 0.5305 0.2370 0.3810 0.1340 0.1865 0.3620 0.0315 0.0255 0.0175 0.0875 0.2930 0.1775
  [49] 0.0755 0.3545 0.2385 0.1335 0.2595 0.2105 0.1730 0.2565 0.1920 0.2765 0.0420 0.2460
  [61] 0.1800 0.3050 0.3020 0.1705 0.2340 0.2340 0.3540 0.4160 0.2135 0.0630 0.2640 0.1405
  [73] 0.4800 0.4740 0.4810 0.4425 0.3625 0.3630 0.2820 0.4695 0.3845 0.5105 0.3960 0.4080
  [85] 0.3800 0.3390 0.4825 0.3305 0.2205 0.3135 0.3410 0.3070 0.4015 0.5070 0.5880 0.5755
  [97] 0.2690 0.2140 0.2010 0.2775 0.1050 0.3280 0.3160 0.3105 0.4975 0.2910 0.2935 0.2610
...remaining output removed...
R สำหรับผู้ใช้ SAS

คำนวณค่าเฉลี่ยและค่ามัธยฐานของ shucked weight

# คำนวณค่าเฉลี่ยของ shuckedWeight
abalone %>%
  pull(shuckedWeight) %>%
  mean()

0.3593675

# คำนวณค่ามัธยฐานของ shuckedWeight
abalone %>%
  pull(shuckedWeight) %>%
  median()

0.336
R สำหรับผู้ใช้ SAS

เลือกตัวแปรสองตัวจาก abalone

# เลือกตัวแปรสองตัว length และ height
abalone %>%
  select(length, height)

# A tibble: 4,177 x 2
   length height
    <dbl>  <dbl>
 1  0.455  0.095
 2  0.35   0.09
 3  0.53   0.135
 4  0.44   0.125
 5  0.33   0.08
 6  0.425  0.095
 7  0.53   0.15
 8  0.545  0.125
# ... with 4,169 more rows
R สำหรับผู้ใช้ SAS

รับค่าสถิติสรุปของ length และ height

summary() แสดงค่าต่ำสุด ค่าสูงสุด ค่าเฉลี่ย ค่ามัธยฐาน และควอร์ไทล์ที่ 25 และ 75

# รับค่าสถิติสรุปของ length และ height
abalone %>%
  select(length, height) %>%
  summary()

     length          height      
 Min.   :0.075   Min.   :0.0000  
 1st Qu.:0.450   1st Qu.:0.1150  
 Median :0.545   Median :0.1400  
 Mean   :0.524   Mean   :0.1395  
 3rd Qu.:0.615   3rd Qu.:0.1650  
 Max.   :0.815   Max.   :1.1300
R สำหรับผู้ใช้ SAS

มาเรียนรู้เกี่ยวกับหอยเป๋าฮื้อกันเถอะ

R สำหรับผู้ใช้ SAS

Preparing Video For Download...