create_vtl_corpus/analysis.R at main · quantling/create_vtl_corpus · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50


# zipf-plots

wav_orig <- dir('vtl_corpus1.0/wav_original/')
tokens_orig <- tools::file_path_sans_ext(wav_orig)  # remove .wav
tokens_orig <- substring(tokens_orig, 8)  # remove numbers in the beginning

types_orig <- table(tokens_orig)
types_orig <- sort(types_orig, decreasing=T)

wav_synth <- dir('vtl_corpus1.0/wav_synthesized/')
tokens_synth <- tools::file_path_sans_ext(wav_synth)  # remove .wav
tokens_synth <- substring(tokens_synth, 8)  # remove numbers in the beginning

types_synth <- table(tokens_synth)
types_synth <- sort(types_synth, decreasing=T)

pdf('figs/zipf.pdf', height=4, width=8)
par(mfcol=c(1,2))

plot(log10(types_orig), log='x', axes=F, xlab="words", ylab="occurrences", main='Recorded')
axis(1, at=1:length(types_orig), labels=names(types_orig))
yticks <- c(1:9, 1:9*10, 1:9*100, 1:2*1000)
axis(2, at=log10(yticks), labels=yticks)

plot(log10(types_synth), log='x', axes=F, xlab="words", ylab="occurrences", main='Synthesized')
axis(1, at=1:length(types_synth), labels=names(types_synth))
yticks <- c(1:9, 1:9*10, 1:9*100, 1:2*1000)
axis(2, at=log10(yticks), labels=yticks)
dev.off()


# durations

dur_orig <- rep(NA, length(wav_orig))

for (ii in 1:length(wav_orig)) {
	name <- wav_orig[ii]
	wav <- readWave(paste0('vtl_corpus1.0/wav_original/', name))
	dur_orig[ii] <- length(wav@left) / wav@samp.rate
}

dur_synth <- rep(NA, length(wav_synth))

for (ii in 1:length(wav_synth)) {
	name <- wav_synth[ii]
	wav <- readWave(paste0('vtl_corpus1.0/wav_synthesized/', name))
	dur_synth[ii] <- length(wav@left) / wav@samp.rate
}