Cài các gói lệnh cần thiết để vẽ hình
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr 1.1.4 ✔ readr 2.1.5
## ✔ forcats 1.0.0 ✔ stringr 1.5.1
## ✔ ggplot2 3.5.2 ✔ tibble 3.2.1
## ✔ lubridate 1.9.4 ✔ tidyr 1.3.1
## ✔ purrr 1.0.4
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library(scales)
##
## Attaching package: 'scales'
##
## The following object is masked from 'package:purrr':
##
## discard
##
## The following object is masked from 'package:readr':
##
## col_factor
library(plyr)
## ------------------------------------------------------------------------------
## You have loaded plyr after dplyr - this is likely to cause problems.
## If you need functions from both plyr and dplyr, please load plyr first, then dplyr:
## library(plyr); library(dplyr)
## ------------------------------------------------------------------------------
##
## Attaching package: 'plyr'
##
## The following objects are masked from 'package:dplyr':
##
## arrange, count, desc, failwith, id, mutate, rename, summarise,
## summarize
##
## The following object is masked from 'package:purrr':
##
## compact
Sử dụng nền trắng đơn giản để vẽ hình
theme_set(theme_bw())
Dẫn đến thư mục chứa phai *.txt
setwd("D:\\VIETNAMESE_CHINESE_TONE_ANALYSIS\\DATA\\TRUNG\\CTONES\\")
Đọc dữ liệu thanh điệu sau khi chay script praatSauc, có tên
Chinese_ma định dạng *.txt
PS <- read.csv ("D:\\VIETNAMESE_CHINESE_TONE_ANALYSIS\\DATA\\TRUNG\\CTONES\\Chinese_ma.txt")
Xem 6 dòng đầu
head(PS)
## Filename var1 var2 var3 Label seg_Start seg_End t t_ms f0
## 1 TQN_ma_ambinh TQN ma ambinh cvo 0.159 0.485 1 0.165 152.479
## 2 TQN_ma_ambinh TQN ma ambinh cvo 0.159 0.485 2 0.170 152.513
## 3 TQN_ma_ambinh TQN ma ambinh cvo 0.159 0.485 3 0.175 152.482
## 4 TQN_ma_ambinh TQN ma ambinh cvo 0.159 0.485 4 0.180 152.450
## 5 TQN_ma_ambinh TQN ma ambinh cvo 0.159 0.485 5 0.185 152.372
## 6 TQN_ma_ambinh TQN ma ambinh cvo 0.159 0.485 6 0.190 152.257
## F1 F2 F3 B1 B2 B3 H1u H2u H4u H2Ku
## 1 706.486 949.108 2656.574 192.924 198.974 663.830 38.666 35.425 44.817 12.582
## 2 735.513 982.621 2747.883 165.743 171.434 591.193 37.997 36.038 46.284 16.570
## 3 743.585 1008.544 2830.994 139.080 121.074 528.836 38.066 36.910 45.952 17.050
## 4 747.667 1032.319 2911.086 122.732 86.962 484.677 38.467 37.267 45.155 15.505
## 5 752.407 1053.926 2970.541 119.537 73.871 433.556 38.560 37.355 45.088 14.562
## 6 757.344 1070.883 2986.533 127.975 71.005 342.683 38.513 37.408 44.776 14.512
## H5Ku A1u A2u A3u H1H2u H2H4u H1A1u H1A2u H1A3u H2KH5Ku
## 1 -0.534 46.253 44.858 15.552 3.241 -9.392 -7.587 -6.192 23.114 13.116
## 2 4.171 48.836 48.227 20.364 1.959 -10.246 -10.839 -10.230 17.633 12.398
## 3 3.645 49.899 49.066 23.395 1.156 -9.042 -11.833 -11.000 14.671 13.406
## 4 3.723 50.024 49.210 24.139 1.200 -7.888 -11.558 -10.743 14.328 11.782
## 5 4.787 49.962 49.518 24.913 1.205 -7.732 -11.402 -10.958 13.647 9.775
## 6 7.033 49.935 50.717 26.335 1.105 -7.368 -11.422 -12.204 12.178 7.478
## H1c H2c H4c A1c A2c A3c H1H2c H2H4c H1A1c H1A2c H1A3c
## 1 38.050 32.824 31.377 28.307 29.607 39.792 5.226 1.447 9.744 8.443 -1.741
## 2 37.418 33.596 33.506 28.978 31.046 43.087 3.822 0.090 8.440 6.373 -5.669
## 3 37.499 34.518 33.260 28.529 29.099 45.232 2.981 1.258 8.970 8.400 -7.733
## 4 37.910 34.917 32.608 27.782 26.758 45.361 2.993 2.309 10.128 11.152 -7.451
## 5 38.015 35.059 32.870 27.673 25.953 45.132 2.957 2.189 10.342 12.062 -7.117
## 6 37.981 35.166 32.969 28.333 26.997 44.207 2.815 2.197 9.648 10.984 -6.226
## CPP HNR05 HNR15 HNR25 HNR35
## 1 16.966 3.914 3.589 3.611 3.569
## 2 21.408 5.029 7.027 7.012 6.964
## 3 25.730 10.067 10.013 9.997 9.923
## 4 25.214 18.742 12.917 12.929 12.813
## 5 28.757 24.671 15.493 15.500 15.359
## 6 29.312 28.323 17.760 17.736 17.587
Đổi tên biến ‘var1’, ‘var2’, ‘var3’ thành ‘Repetition’, ‘Target’,
‘Tone’
colnames(PS)[2:4] <- c('Repetition', 'Target', 'Tone')
Xem 6 dòng đầu
head(PS)
## Filename Repetition Target Tone Label seg_Start seg_End t t_ms
## 1 TQN_ma_ambinh TQN ma ambinh cvo 0.159 0.485 1 0.165
## 2 TQN_ma_ambinh TQN ma ambinh cvo 0.159 0.485 2 0.170
## 3 TQN_ma_ambinh TQN ma ambinh cvo 0.159 0.485 3 0.175
## 4 TQN_ma_ambinh TQN ma ambinh cvo 0.159 0.485 4 0.180
## 5 TQN_ma_ambinh TQN ma ambinh cvo 0.159 0.485 5 0.185
## 6 TQN_ma_ambinh TQN ma ambinh cvo 0.159 0.485 6 0.190
## f0 F1 F2 F3 B1 B2 B3 H1u H2u
## 1 152.479 706.486 949.108 2656.574 192.924 198.974 663.830 38.666 35.425
## 2 152.513 735.513 982.621 2747.883 165.743 171.434 591.193 37.997 36.038
## 3 152.482 743.585 1008.544 2830.994 139.080 121.074 528.836 38.066 36.910
## 4 152.450 747.667 1032.319 2911.086 122.732 86.962 484.677 38.467 37.267
## 5 152.372 752.407 1053.926 2970.541 119.537 73.871 433.556 38.560 37.355
## 6 152.257 757.344 1070.883 2986.533 127.975 71.005 342.683 38.513 37.408
## H4u H2Ku H5Ku A1u A2u A3u H1H2u H2H4u H1A1u H1A2u
## 1 44.817 12.582 -0.534 46.253 44.858 15.552 3.241 -9.392 -7.587 -6.192
## 2 46.284 16.570 4.171 48.836 48.227 20.364 1.959 -10.246 -10.839 -10.230
## 3 45.952 17.050 3.645 49.899 49.066 23.395 1.156 -9.042 -11.833 -11.000
## 4 45.155 15.505 3.723 50.024 49.210 24.139 1.200 -7.888 -11.558 -10.743
## 5 45.088 14.562 4.787 49.962 49.518 24.913 1.205 -7.732 -11.402 -10.958
## 6 44.776 14.512 7.033 49.935 50.717 26.335 1.105 -7.368 -11.422 -12.204
## H1A3u H2KH5Ku H1c H2c H4c A1c A2c A3c H1H2c H2H4c H1A1c
## 1 23.114 13.116 38.050 32.824 31.377 28.307 29.607 39.792 5.226 1.447 9.744
## 2 17.633 12.398 37.418 33.596 33.506 28.978 31.046 43.087 3.822 0.090 8.440
## 3 14.671 13.406 37.499 34.518 33.260 28.529 29.099 45.232 2.981 1.258 8.970
## 4 14.328 11.782 37.910 34.917 32.608 27.782 26.758 45.361 2.993 2.309 10.128
## 5 13.647 9.775 38.015 35.059 32.870 27.673 25.953 45.132 2.957 2.189 10.342
## 6 12.178 7.478 37.981 35.166 32.969 28.333 26.997 44.207 2.815 2.197 9.648
## H1A2c H1A3c CPP HNR05 HNR15 HNR25 HNR35
## 1 8.443 -1.741 16.966 3.914 3.589 3.611 3.569
## 2 6.373 -5.669 21.408 5.029 7.027 7.012 6.964
## 3 8.400 -7.733 25.730 10.067 10.013 9.997 9.923
## 4 11.152 -7.451 25.214 18.742 12.917 12.929 12.813
## 5 12.062 -7.117 28.757 24.671 15.493 15.500 15.359
## 6 10.984 -6.226 29.312 28.323 17.760 17.736 17.587
Chuyển đơn vị thời gian giây sang mili giây cho dễ hiểu
PS$seg_Start <- PS$seg_Start * 1000
PS$seg_End <- PS$seg_End * 1000
PS$t_ms <- PS$t_ms * 1000
Xem 6 dòng đầu
head(PS)
## Filename Repetition Target Tone Label seg_Start seg_End t t_ms f0
## 1 TQN_ma_ambinh TQN ma ambinh cvo 159 485 1 165 152.479
## 2 TQN_ma_ambinh TQN ma ambinh cvo 159 485 2 170 152.513
## 3 TQN_ma_ambinh TQN ma ambinh cvo 159 485 3 175 152.482
## 4 TQN_ma_ambinh TQN ma ambinh cvo 159 485 4 180 152.450
## 5 TQN_ma_ambinh TQN ma ambinh cvo 159 485 5 185 152.372
## 6 TQN_ma_ambinh TQN ma ambinh cvo 159 485 6 190 152.257
## F1 F2 F3 B1 B2 B3 H1u H2u H4u H2Ku
## 1 706.486 949.108 2656.574 192.924 198.974 663.830 38.666 35.425 44.817 12.582
## 2 735.513 982.621 2747.883 165.743 171.434 591.193 37.997 36.038 46.284 16.570
## 3 743.585 1008.544 2830.994 139.080 121.074 528.836 38.066 36.910 45.952 17.050
## 4 747.667 1032.319 2911.086 122.732 86.962 484.677 38.467 37.267 45.155 15.505
## 5 752.407 1053.926 2970.541 119.537 73.871 433.556 38.560 37.355 45.088 14.562
## 6 757.344 1070.883 2986.533 127.975 71.005 342.683 38.513 37.408 44.776 14.512
## H5Ku A1u A2u A3u H1H2u H2H4u H1A1u H1A2u H1A3u H2KH5Ku
## 1 -0.534 46.253 44.858 15.552 3.241 -9.392 -7.587 -6.192 23.114 13.116
## 2 4.171 48.836 48.227 20.364 1.959 -10.246 -10.839 -10.230 17.633 12.398
## 3 3.645 49.899 49.066 23.395 1.156 -9.042 -11.833 -11.000 14.671 13.406
## 4 3.723 50.024 49.210 24.139 1.200 -7.888 -11.558 -10.743 14.328 11.782
## 5 4.787 49.962 49.518 24.913 1.205 -7.732 -11.402 -10.958 13.647 9.775
## 6 7.033 49.935 50.717 26.335 1.105 -7.368 -11.422 -12.204 12.178 7.478
## H1c H2c H4c A1c A2c A3c H1H2c H2H4c H1A1c H1A2c H1A3c
## 1 38.050 32.824 31.377 28.307 29.607 39.792 5.226 1.447 9.744 8.443 -1.741
## 2 37.418 33.596 33.506 28.978 31.046 43.087 3.822 0.090 8.440 6.373 -5.669
## 3 37.499 34.518 33.260 28.529 29.099 45.232 2.981 1.258 8.970 8.400 -7.733
## 4 37.910 34.917 32.608 27.782 26.758 45.361 2.993 2.309 10.128 11.152 -7.451
## 5 38.015 35.059 32.870 27.673 25.953 45.132 2.957 2.189 10.342 12.062 -7.117
## 6 37.981 35.166 32.969 28.333 26.997 44.207 2.815 2.197 9.648 10.984 -6.226
## CPP HNR05 HNR15 HNR25 HNR35
## 1 16.966 3.914 3.589 3.611 3.569
## 2 21.408 5.029 7.027 7.012 6.964
## 3 25.730 10.067 10.013 9.997 9.923
## 4 25.214 18.742 12.917 12.929 12.813
## 5 28.757 24.671 15.493 15.500 15.359
## 6 29.312 28.323 17.760 17.736 17.587
Tách phụ âm “m” ở cột cuối
PS$Onset <- substring(PS$Target, 1, 1)
PS$Onset <- as.factor(PS$Onset)
Xem 6 dòng đầu
head(PS)
## Filename Repetition Target Tone Label seg_Start seg_End t t_ms f0
## 1 TQN_ma_ambinh TQN ma ambinh cvo 159 485 1 165 152.479
## 2 TQN_ma_ambinh TQN ma ambinh cvo 159 485 2 170 152.513
## 3 TQN_ma_ambinh TQN ma ambinh cvo 159 485 3 175 152.482
## 4 TQN_ma_ambinh TQN ma ambinh cvo 159 485 4 180 152.450
## 5 TQN_ma_ambinh TQN ma ambinh cvo 159 485 5 185 152.372
## 6 TQN_ma_ambinh TQN ma ambinh cvo 159 485 6 190 152.257
## F1 F2 F3 B1 B2 B3 H1u H2u H4u H2Ku
## 1 706.486 949.108 2656.574 192.924 198.974 663.830 38.666 35.425 44.817 12.582
## 2 735.513 982.621 2747.883 165.743 171.434 591.193 37.997 36.038 46.284 16.570
## 3 743.585 1008.544 2830.994 139.080 121.074 528.836 38.066 36.910 45.952 17.050
## 4 747.667 1032.319 2911.086 122.732 86.962 484.677 38.467 37.267 45.155 15.505
## 5 752.407 1053.926 2970.541 119.537 73.871 433.556 38.560 37.355 45.088 14.562
## 6 757.344 1070.883 2986.533 127.975 71.005 342.683 38.513 37.408 44.776 14.512
## H5Ku A1u A2u A3u H1H2u H2H4u H1A1u H1A2u H1A3u H2KH5Ku
## 1 -0.534 46.253 44.858 15.552 3.241 -9.392 -7.587 -6.192 23.114 13.116
## 2 4.171 48.836 48.227 20.364 1.959 -10.246 -10.839 -10.230 17.633 12.398
## 3 3.645 49.899 49.066 23.395 1.156 -9.042 -11.833 -11.000 14.671 13.406
## 4 3.723 50.024 49.210 24.139 1.200 -7.888 -11.558 -10.743 14.328 11.782
## 5 4.787 49.962 49.518 24.913 1.205 -7.732 -11.402 -10.958 13.647 9.775
## 6 7.033 49.935 50.717 26.335 1.105 -7.368 -11.422 -12.204 12.178 7.478
## H1c H2c H4c A1c A2c A3c H1H2c H2H4c H1A1c H1A2c H1A3c
## 1 38.050 32.824 31.377 28.307 29.607 39.792 5.226 1.447 9.744 8.443 -1.741
## 2 37.418 33.596 33.506 28.978 31.046 43.087 3.822 0.090 8.440 6.373 -5.669
## 3 37.499 34.518 33.260 28.529 29.099 45.232 2.981 1.258 8.970 8.400 -7.733
## 4 37.910 34.917 32.608 27.782 26.758 45.361 2.993 2.309 10.128 11.152 -7.451
## 5 38.015 35.059 32.870 27.673 25.953 45.132 2.957 2.189 10.342 12.062 -7.117
## 6 37.981 35.166 32.969 28.333 26.997 44.207 2.815 2.197 9.648 10.984 -6.226
## CPP HNR05 HNR15 HNR25 HNR35 Onset
## 1 16.966 3.914 3.589 3.611 3.569 m
## 2 21.408 5.029 7.027 7.012 6.964 m
## 3 25.730 10.067 10.013 9.997 9.923 m
## 4 25.214 18.742 12.917 12.929 12.813 m
## 5 28.757 24.671 15.493 15.500 15.359 m
## 6 29.312 28.323 17.760 17.736 17.587 m
Tách nguyên tâm chính ở cột cuối cùng
PS$Vowel <- as.factor(substring(PS$Target, 2, 2))
PS$Vowel <- as.factor(PS$Vowel)
Xem 6 dòng đầu
head(PS)
## Filename Repetition Target Tone Label seg_Start seg_End t t_ms f0
## 1 TQN_ma_ambinh TQN ma ambinh cvo 159 485 1 165 152.479
## 2 TQN_ma_ambinh TQN ma ambinh cvo 159 485 2 170 152.513
## 3 TQN_ma_ambinh TQN ma ambinh cvo 159 485 3 175 152.482
## 4 TQN_ma_ambinh TQN ma ambinh cvo 159 485 4 180 152.450
## 5 TQN_ma_ambinh TQN ma ambinh cvo 159 485 5 185 152.372
## 6 TQN_ma_ambinh TQN ma ambinh cvo 159 485 6 190 152.257
## F1 F2 F3 B1 B2 B3 H1u H2u H4u H2Ku
## 1 706.486 949.108 2656.574 192.924 198.974 663.830 38.666 35.425 44.817 12.582
## 2 735.513 982.621 2747.883 165.743 171.434 591.193 37.997 36.038 46.284 16.570
## 3 743.585 1008.544 2830.994 139.080 121.074 528.836 38.066 36.910 45.952 17.050
## 4 747.667 1032.319 2911.086 122.732 86.962 484.677 38.467 37.267 45.155 15.505
## 5 752.407 1053.926 2970.541 119.537 73.871 433.556 38.560 37.355 45.088 14.562
## 6 757.344 1070.883 2986.533 127.975 71.005 342.683 38.513 37.408 44.776 14.512
## H5Ku A1u A2u A3u H1H2u H2H4u H1A1u H1A2u H1A3u H2KH5Ku
## 1 -0.534 46.253 44.858 15.552 3.241 -9.392 -7.587 -6.192 23.114 13.116
## 2 4.171 48.836 48.227 20.364 1.959 -10.246 -10.839 -10.230 17.633 12.398
## 3 3.645 49.899 49.066 23.395 1.156 -9.042 -11.833 -11.000 14.671 13.406
## 4 3.723 50.024 49.210 24.139 1.200 -7.888 -11.558 -10.743 14.328 11.782
## 5 4.787 49.962 49.518 24.913 1.205 -7.732 -11.402 -10.958 13.647 9.775
## 6 7.033 49.935 50.717 26.335 1.105 -7.368 -11.422 -12.204 12.178 7.478
## H1c H2c H4c A1c A2c A3c H1H2c H2H4c H1A1c H1A2c H1A3c
## 1 38.050 32.824 31.377 28.307 29.607 39.792 5.226 1.447 9.744 8.443 -1.741
## 2 37.418 33.596 33.506 28.978 31.046 43.087 3.822 0.090 8.440 6.373 -5.669
## 3 37.499 34.518 33.260 28.529 29.099 45.232 2.981 1.258 8.970 8.400 -7.733
## 4 37.910 34.917 32.608 27.782 26.758 45.361 2.993 2.309 10.128 11.152 -7.451
## 5 38.015 35.059 32.870 27.673 25.953 45.132 2.957 2.189 10.342 12.062 -7.117
## 6 37.981 35.166 32.969 28.333 26.997 44.207 2.815 2.197 9.648 10.984 -6.226
## CPP HNR05 HNR15 HNR25 HNR35 Onset Vowel
## 1 16.966 3.914 3.589 3.611 3.569 m a
## 2 21.408 5.029 7.027 7.012 6.964 m a
## 3 25.730 10.067 10.013 9.997 9.923 m a
## 4 25.214 18.742 12.917 12.929 12.813 m a
## 5 28.757 24.671 15.493 15.500 15.359 m a
## 6 29.312 28.323 17.760 17.736 17.587 m a
Kiểm tra cấu trúc dữ liệu
xtabs(~Tone+Repetition, PS)
## Repetition
## Tone TQN
## ambinh 65
## duongbinh 73
## khuthanh 35
## thuongthanh 85
xtabs(~Tone+Target, PS)
## Target
## Tone ma
## ambinh 65
## duongbinh 73
## khuthanh 35
## thuongthanh 85
Căn chỉnh dữ liệu sao cho tất cả các đường thanh điệu bắt đầu tại t
= 1
PS$min_t <- ave(PS$t, PS$Filename, FUN=min)
PS$corr_t <- PS$t - PS$min_t + 1
Xem 6 dòng đầu
head(PS)
## Filename Repetition Target Tone Label seg_Start seg_End t t_ms f0
## 1 TQN_ma_ambinh TQN ma ambinh cvo 159 485 1 165 152.479
## 2 TQN_ma_ambinh TQN ma ambinh cvo 159 485 2 170 152.513
## 3 TQN_ma_ambinh TQN ma ambinh cvo 159 485 3 175 152.482
## 4 TQN_ma_ambinh TQN ma ambinh cvo 159 485 4 180 152.450
## 5 TQN_ma_ambinh TQN ma ambinh cvo 159 485 5 185 152.372
## 6 TQN_ma_ambinh TQN ma ambinh cvo 159 485 6 190 152.257
## F1 F2 F3 B1 B2 B3 H1u H2u H4u H2Ku
## 1 706.486 949.108 2656.574 192.924 198.974 663.830 38.666 35.425 44.817 12.582
## 2 735.513 982.621 2747.883 165.743 171.434 591.193 37.997 36.038 46.284 16.570
## 3 743.585 1008.544 2830.994 139.080 121.074 528.836 38.066 36.910 45.952 17.050
## 4 747.667 1032.319 2911.086 122.732 86.962 484.677 38.467 37.267 45.155 15.505
## 5 752.407 1053.926 2970.541 119.537 73.871 433.556 38.560 37.355 45.088 14.562
## 6 757.344 1070.883 2986.533 127.975 71.005 342.683 38.513 37.408 44.776 14.512
## H5Ku A1u A2u A3u H1H2u H2H4u H1A1u H1A2u H1A3u H2KH5Ku
## 1 -0.534 46.253 44.858 15.552 3.241 -9.392 -7.587 -6.192 23.114 13.116
## 2 4.171 48.836 48.227 20.364 1.959 -10.246 -10.839 -10.230 17.633 12.398
## 3 3.645 49.899 49.066 23.395 1.156 -9.042 -11.833 -11.000 14.671 13.406
## 4 3.723 50.024 49.210 24.139 1.200 -7.888 -11.558 -10.743 14.328 11.782
## 5 4.787 49.962 49.518 24.913 1.205 -7.732 -11.402 -10.958 13.647 9.775
## 6 7.033 49.935 50.717 26.335 1.105 -7.368 -11.422 -12.204 12.178 7.478
## H1c H2c H4c A1c A2c A3c H1H2c H2H4c H1A1c H1A2c H1A3c
## 1 38.050 32.824 31.377 28.307 29.607 39.792 5.226 1.447 9.744 8.443 -1.741
## 2 37.418 33.596 33.506 28.978 31.046 43.087 3.822 0.090 8.440 6.373 -5.669
## 3 37.499 34.518 33.260 28.529 29.099 45.232 2.981 1.258 8.970 8.400 -7.733
## 4 37.910 34.917 32.608 27.782 26.758 45.361 2.993 2.309 10.128 11.152 -7.451
## 5 38.015 35.059 32.870 27.673 25.953 45.132 2.957 2.189 10.342 12.062 -7.117
## 6 37.981 35.166 32.969 28.333 26.997 44.207 2.815 2.197 9.648 10.984 -6.226
## CPP HNR05 HNR15 HNR25 HNR35 Onset Vowel min_t corr_t
## 1 16.966 3.914 3.589 3.611 3.569 m a 1 1
## 2 21.408 5.029 7.027 7.012 6.964 m a 1 2
## 3 25.730 10.067 10.013 9.997 9.923 m a 1 3
## 4 25.214 18.742 12.917 12.929 12.813 m a 1 4
## 5 28.757 24.671 15.493 15.500 15.359 m a 1 5
## 6 29.312 28.323 17.760 17.736 17.587 m a 1 6
Biểu đồ mỗi âm tiết
f0chart <- PS %>%
filter(Tone %in% c("ngang", "huyền", "sắc", "hỏi", "ngã", "nặng")) %>%
ggplot(aes(x = corr_t, y = f0, color = Tone, linetype = Tone, group = Filename)) +
geom_line(size = 0.5) +
labs(x = "Thời gian từ đầu nguyên âm (ms)", y = "f0 (Hz)", title = "Đường thanh điệu 6 thanh tiếng Trung") +
scale_color_manual(values = c("#e41a1c", "#377eb8", "#4daf4a", "#984ea3", "#ff7f00", "black")) +
scale_linetype_manual(values = c("solid", "dashed", "solid", "dashed", "solid", "dashed")) +
theme_minimal() +
theme(plot.title = element_text(hjust = 0.5))
## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
Xoá các dòng F0 = 0 để tránh ảnh hưởng biểu đồ
PS <- subset(PS, f0 != 0)
Vẽ lại biểu đồ mà không có các giá trị bằng 0
f0chart <- PS %>%
filter(Tone %in% c("ngang", "huyền", "sắc", "hỏi", "ngã", "nặng")) %>%
ggplot(aes(x = corr_t, y = f0, color = Tone, linetype = Tone, group = Filename)) +
geom_line(size = 0.5) +
labs(x = "Thời gian từ đầu nguyên âm (ms)", y = "f0 (Hz)", title = "Đường thanh điệu 6 thanh tiếng Trung") +
scale_color_manual(values = c("#e41a1c", "#377eb8", "#4daf4a", "#984ea3", "#ff7f00", "black")) +
scale_linetype_manual(values = c("solid", "dashed", "solid", "dashed", "solid", "dashed")) +
theme_minimal() +
theme(plot.title = element_text(hjust = 0.5))
Đảm bảo onset và vowel
PS$Onset <- as.factor(substring(PS$Target, 1, 1))
PS$Vowel <- as.factor(substring(PS$Target, 2, 2))
Loại bỏ các dòng có giá trị NA ở Onset hoặc Vowel
PS$Onset <- as.factor(substring(PS$Target, 1, 1))
PS$Vowel <- as.factor(substring(PS$Target, 2, 2))
Phân tích theo âm đầu vKiểm tra dữ liệu có đa dạng về Onset và
Vowelà nguyên âm
table(PS$Onset)
##
## m
## 217
table(PS$Vowel)
##
## a
## 217
Tạo biểu đồ với facet_grid()
table(PS$Onset)
##
## m
## 217
table(PS$Vowel)
##
## a
## 217
Vẽ biểu đồ
f0chart <- PS %>%
filter(Tone %in% c("ambinh", "duongbinh", "khuthanh", "thuongthanh")) %>%
ggplot(aes(x = corr_t, y = f0, color = Tone, linetype = Tone, group = Filename)) +
geom_line(size = 0.5) +
facet_grid(Onset ~ Vowel) +
labs(x = "Thời gian từ đầu nguyên âm (ms)", y = "f0 (Hz)", title = "Đường thanh điệu 4 thanh tiếng Trung") +
scale_color_manual(values = c("#e41a1c", "#377eb8", "#4daf4a", "#984ea3")) +
scale_linetype_manual(values = c("solid", "dashed", "solid", "dashed")) +
theme_minimal() +
theme(plot.title = element_text(hjust = 0.5))
print(f0chart)
