f@0
|
1 library(tidyverse)
|
f@0
|
2
|
f@0
|
3 feats_mode <- 'original'
|
f@0
|
4 feats_lev <- 'f'
|
f@0
|
5
|
f@0
|
6 feats_path <- paste0('~/Repositories/research/datasets/features/gtzan/',
|
f@0
|
7 feats_mode, '/scattering/')
|
f@0
|
8
|
f@0
|
9 feats <- read.csv(file = paste0(feats_path, feats_lev, '.csv'))
|
f@0
|
10
|
f@0
|
11 feats_desc <-
|
f@0
|
12 feats %>%
|
f@0
|
13 group_by(ex_id, gt) %>%
|
f@0
|
14 summarise_all(funs(mean, median, var, min, max,
|
f@0
|
15 skew = psych::skew,
|
f@0
|
16 kurt = psych::kurtosi)) %>%
|
f@0
|
17 ungroup()
|
f@0
|
18
|
f@0
|
19 dif <- function(x) x - lag(x)
|
f@0
|
20
|
f@0
|
21 feats_dif <-
|
f@0
|
22 feats %>%
|
f@0
|
23 group_by(ex_id, gt) %>%
|
f@0
|
24 mutate_all(funs(dif)) %>%
|
f@0
|
25 na.omit() %>%
|
f@0
|
26 ungroup()
|
f@0
|
27 names(feats_dif) <- c(names(feats_dif)[1:2],
|
f@0
|
28 paste0(names(feats_dif)[3:length(names(feats_dif))],
|
f@0
|
29 "_dif"))
|
f@0
|
30
|
f@0
|
31 feats_dif2 <-
|
f@0
|
32 feats_dif %>%
|
f@0
|
33 group_by(ex_id, gt) %>%
|
f@0
|
34 mutate_all(funs(dif)) %>%
|
f@0
|
35 na.omit() %>%
|
f@0
|
36 ungroup()
|
f@0
|
37 names(feats_dif2) <- c(names(feats_dif2)[1:2],
|
f@0
|
38 paste0(names(feats_dif2)[3:length(names(feats_dif2))],
|
f@0
|
39 "2"))
|
f@0
|
40
|
f@0
|
41 feats_desc <-
|
f@0
|
42 cbind(feats_desc,
|
f@0
|
43 feats_dif %>%
|
f@0
|
44 group_by(ex_id, gt) %>%
|
f@0
|
45 summarise_all(funs(mean, var)) %>%
|
f@0
|
46 ungroup() %>%
|
f@0
|
47 select(-ex_id, -gt),
|
f@0
|
48 feats_dif2 %>%
|
f@0
|
49 group_by(ex_id, gt) %>%
|
f@0
|
50 summarise_all(funs(mean, var)) %>%
|
f@0
|
51 ungroup() %>%
|
f@0
|
52 select(-ex_id, -gt)
|
f@0
|
53 )
|
f@0
|
54
|
f@0
|
55
|
f@0
|
56 file_name <- paste0(feats_path, feats_lev, '_desc.csv')
|
f@0
|
57 write.csv(x = feats_desc, file = file_name, row.names = F)
|
f@0
|
58
|