Cara Praktis Membuat Diagram XY dan Menambah Kesalahan Standar

Halo semuanya,

dalam tulisan singkat ini, saya akan membagikan cara membuat plot XY dengan menambahkan kesalahan standard dari Y.

Misalkan kita mempunyai dataset sebagai berikut:

> # input data https://catalog.data.gov/dataset/in-vivo-plasma-concentration-for-lindane-after-6-hour-exposure-in-human-skin
> time <- c(1.029411765,
+           2.205882353,
+           3.235294118,
+           4.117647059,
+           5.073529412,
+           6.102941176,
+           7.205882353,
+           8.161764706,
+           12.05882353,
+           24.19117647,
+           32.42647059,
+           48.38235294,
+           79.04411765
+ )
> cp <- c(0.118055556,
+           0.309027778,
+           0.295138889,
+           0.309027778,
+           0.302083333,
+           0.336805556,
+           0.392361111,
+           0.395833333,
+           0.340277778,
+           0.173611111,
+           0.180555556,
+           0.131944444,
+           0.076388889
+ )
> sem <- c(0.045861286,
+          0.201592285,
+          0.208510595,
+          0.2209973,
+          0.21122832,
+          0.2431814,
+          0.248706332,
+          0.272317046,
+          0.24019267,
+          0.138240919,
+          0.136563753,
+          0.076203673,
+          0.048711957
+   )
> dataset01 <- data.frame(cbind(time, cp, sem))
> dataset01
        time         cp        sem
1   1.029412 0.11805556 0.04586129
2   2.205882 0.30902778 0.20159229
3   3.235294 0.29513889 0.20851059
4   4.117647 0.30902778 0.22099730
5   5.073529 0.30208333 0.21122832
6   6.102941 0.33680556 0.24318140
7   7.205882 0.39236111 0.24870633
8   8.161765 0.39583333 0.27231705
9  12.058824 0.34027778 0.24019267
10 24.191176 0.17361111 0.13824092
11 32.426471 0.18055556 0.13656375
12 48.382353 0.13194444 0.07620367
13 79.044118 0.07638889 0.04871196

# time adalah waktu (jam), cp adalah kadar obat dalam plasma (mg/L), 
# dan sem adalah standard error measurement dari cp

maka kita dapat membuat plot dengan standard error menggunakan ggplot sebagai berikut:

> library(ggplot2)
> ggplot(dataset01, aes(x = time, y = cp)) +
+   geom_point(size = 3) +
+   geom_line() +
+   geom_errorbar(aes(ymin = cp - sem, ymax = cp + sem), width = 0.1) +
+   labs(x = "Time", y = "CP", title = "") +
+   theme_bw()

Mudah sekali bukan?

Demikian yang dapat saya bagikan dalam tulisan singkat ini. Semoga bermanfaat.