создавая одну и ту же переменную в разных фреймах данных - PullRequest
Новая
       14

создавая одну и ту же переменную в разных фреймах данных

2 голосов
/ 03 октября 2019

У меня 14 фреймов данных с одинаковыми именами столбцов и одинаковым количеством столбцов. Я хотел бы создать одну и ту же переменную (hi130_eur = hi130 / fixrate) в каждом из этих фреймов данных.

Сначала я попробовал с циклом: 

countries<-list(country2, country3, country4, country6, country8, country10, country11, country12,
country13, country14, country15, country51, country55, country57)

for(i in 1:length(countries)) { 
     countries[i] <- countries[i] %>% group_by(wave) %>% mutate(hi130_eur=hi130/fixrate) 
  }

, где wave в group_by - это имя столбца в каждом data.frame, определяющем количество волн обзора. Тем не менее, он сообщает о следующей ошибке:

Ошибка в UseMethod ("group_by_"): нет применимого метода для group_by_, примененного к объекту класса "list"

Итак, я попробовал также с lapply следующим образом: 

lapply(countries, function(x) { 
x<-x %>% group_by(wave) %>% mutate(hi100_eur=hi130/fixrate)
   cbind(x, hi130_eur)
   return(x)
   })

, но сообщает, что hi130_eur не найден.
Есть какие-либо предложения о том, как правильно создать одну и ту же переменную в каждом фрейме данных?

Каждый фрейм данных имеет следующую структуру: 

 structure(list(country = c(2L, 2L, 2L, 2L, 2L, 2L), wave = c(1L, 
1L, 1L, 1L, 1L, 1L), hid = c(7601L, 8401L, 10701L, 15701L, 15701L, 
16501L), hg004 = c(1.07908, 1.47759, 1.24982, 0.94703001, 0.94703001, 
0.92088997), hg005 = c(-9L, -9L, -9L, -9L, -9L, -9L), hg006 = c(-9L, 
-9L, -9L, -9L, -9L, -9L), hg007 = c(-9L, -9L, -9L, -9L, -9L, 
-9L), hg014 = c(1994L, 1994L, 1994L, 1994L, 1994L, 1994L), hd001 = c(1L, 
1L, 1L, 3L, 3L, 1L), hd002 = c(1L, 1L, 1L, 2L, 2L, 1L), hd003 = c(1L, 
1L, 1L, 2L, 2L, 1L), hd004 = c(1, 1, 1, 2.2, 2.2, 1), hd005 = c(1, 
1, 1, 1.8, 1.8, 1), hd006 = c(2L, 1L, 3L, 8L, 8L, 2L), hd006a = c(2L, 
3L, 4L, 11L, 11L, 5L), hd006b = c(1L, 2L, 3L, 11L, 11L, 3L), 
hd007 = c(-8L, -8L, -8L, -8L, -8L, -8L), hd008 = c(-8L, -8L, 
-8L, -8L, -8L, -8L), hd009 = c(-8L, -8L, -8L, -8L, -8L, -8L
), hd010 = c(-8L, -8L, -8L, -8L, -8L, -8L), hi001 = c(4L, 
1L, 1L, 1L, 1L, 1L), hi020 = c(0.58311999, 0.57414001, 0.57635999, 
0.6074, 0.6074, 0.56733), hi100 = c(97000L, 126550L, 90000L, 
249500L, 249500L, 156166L), hi110 = c(13000L, 106000L, 90000L, 
185500L, 185500L, 156166L), hi111 = c(13000L, 106000L, 90000L, 
185500L, 185500L, 156166L), hi112 = c(0L, 0L, 0L, 0L, 0L, 
0L), hi120 = c(0L, 0L, 0L, 0L, 0L, 0L), hi121 = c(0L, 0L, 
0L, 0L, 0L, 0L), hi122 = c(0L, 0L, 0L, 0L, 0L, 0L), hi122g = c(0L, 
0L, 0L, 0L, 0L, 0L), hi123 = c(0L, 0L, 0L, 0L, 0L, 0L), hi130 = c(84000L, 
20550L, 0L, 64000L, 64000L, 0L), hi131 = c(84000L, 0L, 0L, 
0L, 0L, 0L), hi132 = c(0L, 6150L, 0L, 0L, 0L, 0L), hi133 = c(0L, 
0L, 0L, 64000L, 64000L, 0L), hi134 = c(0L, 0L, 0L, 0L, 0L, 
0L), hi135 = c(0L, 0L, 0L, 0L, 0L, 0L), hi136 = c(0L, 0L, 
0L, 0L, 0L, 0L), hi137 = c(0L, 0L, 0L, 0L, 0L, 0L), hi138 = c(0L, 
14400L, 0L, 0L, 0L, 0L), ha005 = c(4L, 4L, 4L, 4L, 4L, 2L
), ha006 = c(2L, 2L, 2L, 3L, 3L, 3L), ha007 = c(3L, 3L, 3L, 
4L, 4L, 4L), ha008 = c(1L, 1L, 1L, 1L, 1L, 1L), ha009 = c(1L, 
1L, 1L, 1L, 1L, 1L), ha010 = c(1L, 1L, 1L, 1L, 1L, 1L), ha011 = c(1L, 
1L, 1L, 1L, 1L, 1L), ha012 = c(2L, 1L, 1L, 1L, 1L, 1L), ha013 = c(2L, 
1L, 2L, 2L, 2L, 1L), ha014 = c(1L, 2L, 2L, 1L, 1L, 2L), ha015b = c(-8L, 
-8L, -8L, -8L, -8L, -8L), ha022 = c(2L, 2L, 2L, 2L, 2L, 2L
), ha023 = c(1L, 2L, 2L, 1L, 1L, 1L), ha024a = c(1L, -8L, 
-8L, 1L, 1L, 1L), ha024b = c(1180L, -8L, -8L, 1000L, 1000L, 
800L), ha025 = c(-8L, -8L, -8L, -8L, -8L, -8L), ha026 = c(-8L, 
1L, 1L, -8L, -8L, -8L), ha027 = c(-8L, 3400L, 2150L, -8L, 
-8L, -8L), ha029 = c(-8L, 1L, 2L, -8L, -8L, -8L), ha030 = c(-8L, 
2L, 2L, -8L, -8L, -8L), ha031 = c(-8L, 1L, 1L, -8L, -8L, 
-8L), ha032 = c(-8L, 2L, 2L, -8L, -8L, -8L), ha033 = c(-8L, 
2L, 2L, -8L, -8L, -8L), ha034 = c(-8L, 2L, 2L, -8L, -8L, 
-8L), ha035 = c(-8L, 2L, 2L, -8L, -8L, -8L), hl001 = c(-8L, 
-8L, -8L, -8L, -8L, -8L), hl002 = c(-8L, -8L, -8L, -8L, -8L, 
-8L), hl003 = c(-8L, -8L, -8L, -8L, -8L, -8L), pid = c(76101, 
84101, 107101, 157101, 157102, 165101), pg002 = c(1.07218, 
1.4762599, 1.26633, 0.94999999, 0.94999999, 0.87747997), 
pg003 = c(1.07218, 1.4762599, 1.26633, 0.94999999, 0.94999999, 
0.87747997), pd003 = c(42L, 70L, 22L, 33L, 33L, 46L), pd004 = c(1L, 
1L, 2L, 2L, 1L, 2L), pe001 = c(7L, 1L, 7L, 1L, 1L, 1L), pe003 = c(5L, 
1L, 3L, 1L, 1L, 1L), pe004 = c(-8L, 1L, -8L, 1L, 1L, 1L), 
pe005 = c(-8L, 30L, -8L, 37L, 37L, 37L), pe005b = c(-8L, 
-8L, -8L, -8L, -8L, -8L), pe005c = c(-8L, 1L, -8L, 1L, 1L, 
1L), pe014 = c(-8L, -8L, -8L, 2L, 2L, -8L), pe015 = c(-8L, 
-8L, -8L, 0L, 0L, -8L), pe024 = c(-8L, -8L, -8L, -8L, -8L, 
-8L), pe025 = c(-8L, -8L, -8L, -8L, -8L, -8L), pe026 = c(-8L, 
2L, -8L, 2L, 2L, 2L), pe030 = c(-8L, 2L, -8L, 2L, 2L, 2L), 
pu001 = c(1L, 2L, 2L, 2L, 2L, 2L), pu004 = c(1L, -8L, 2L, 
-8L, -8L, -8L), pu004a = c(1L, -8L, 2L, -8L, -8L, -8L), ps002 = c(20L, 
-8L, 20L, -8L, -8L, -8L), ps006 = c(-8L, -8L, -8L, -8L, -8L, 
-8L), pi001 = c(4L, 1L, 1L, 5L, 1L, 1L), pi100 = c(97000L, 
126550L, 90000L, 118000L, 131500L, 156166L), pi110 = c(13000L, 
106000L, 90000L, 54000L, 131500L, 156166L), pi112 = c(0L, 
0L, 0L, 0L, 0L, 0L), pi121 = c(0L, 0L, 0L, 0L, 0L, 0L), pi122a = c(0L, 
0L, 0L, 0L, 0L, 0L), pi123 = c(0L, 0L, 0L, 0L, 0L, 0L), pi130 = c(84000L, 
20550L, 0L, 64000L, 0L, 0L), pi131 = c(84000L, 0L, 0L, 0L, 
0L, 0L), pi132 = c(0L, 6150L, 0L, 0L, 0L, 0L), pi133 = c(0L, 
0L, 0L, 64000L, 0L, 0L), pi134 = c(0L, 0L, 0L, 0L, 0L, 0L
), pi135 = c(0L, 0L, 0L, 0L, 0L, 0L), pi136 = c(0L, 0L, 0L, 
0L, 0L, 0L), pi137a = c(0L, 0L, 0L, 0L, 0L, 0L), pi138a = c(0L, 
14400L, 0L, 0L, 0L, 0L), pt022 = c(2L, 2L, 2L, 2L, 2L, 1L
), pt023 = c(-9L, 20L, 19L, 28L, 20L, 23L), pt024 = c(-8L, 
-8L, -8L, -8L, -8L, -8L), ph006 = c(2L, 2L, 1L, 2L, 2L, 2L
), ph007 = c(0L, 0L, 2L, 0L, 0L, 0L), ph008 = c(-8L, -8L, 
-8L, -8L, -8L, -8L), ph009 = c(-8L, -8L, -8L, -8L, -8L, -8L
), ph012 = c(-8L, -8L, -8L, -8L, -8L, -8L), ph013 = c(3L, 
3L, 3L, 3L, 3L, 1L), ph022 = c(-8, -8, -8, -8, -8, -8), pr006 = c(4L, 
4L, 4L, 1L, 1L, 2L), pr007 = c(-8L, -8L, -8L, 3L, 3L, -8L
), pr008 = c(-8L, -8L, -8L, -8L, -8L, 1L), weight_prsn = c(1.07908, 
1.47759, 1.24982, 2.84109, 2.84109, 0.92088997), fixrate = c(1, 
1, 1, 1, 1, 1), hi100_eur = c(97000, 126550, 90000, 249500, 
249500, 156166), pareto = c(378337.646856238, 645715.836409999, 
373027.903072689, 392088.525948343, 459127.982320312, 363013.212050222
), gpd = c(NA, NA, NA, NA, NA, NA)), row.names = c(NA, 6L
), class = "data.frame")

1 Ответ

1 голос
/ 03 октября 2019

Принимая во внимание комментарии, я решил переписать то, что, по моему мнению, может решить вашу проблему следующим образом.

Примеры данных

В ваших наборах данных есть 124 переменных, но для демонстрации потенциального решения вашей проблемы достаточно country, wave, hi130, and fixrate. Следовательно, два примера объектов фрейма данных (country1 и country2) показаны ниже. 

#dput(country1)
structure(list(country = c(2L, 2L, 2L, 2L, 2L, 2L), wave = c(1L, 
1L, 1L, 1L, 1L, 1L), hi130 = c(84000L, 20550L, 0L, 64000L, 64000L, 
0L), fixrate = c(1, 1, 1, 1, 1, 1)), class = "data.frame", row.names = c(NA, 
6L))
# -------------------------------------------------------------------------
#   country wave hi130 fixrate
# 1       2    1 84000       1
# 2       2    1 20550       1
# 3       2    1     0       1
# 4       2    1 64000       1
# 5       2    1 64000       1
# 6       2    1     0       1
# -------------------------------------------------------------------------
#dput(country2)
structure(list(country = c(1, 1, 1, 1, 1, 1), wave = c(2, 2, 
2, 2, 2, 2), hi130 = c(0, 59800, 20440, 19181, 121213, 0), fixrate = c(3, 
3, 3, 3, 3, 3)), class = "data.frame", row.names = c(NA, -6L))
# -------------------------------------------------------------------------
#   country wave  hi130 fixrate
# 1       1    2      0       3
# 2       1    2  59800       3
# 3       1    2  20440       3
# 4       1    2  19181       3
# 5       1    2 121213       3
# 6       1    2      0       3

Функция

Следующая функция используется для получения желаемого выхода. 

library(dplyr)
my_func <- function(df, grp="wave", hi130="hi130", fixrate="fixrate") {
  df %>%
    group_by_(.dots = grp) %>%
    mutate(hi130_eur=(hi130/fixrate))

}
#Example usage
my_func(country1)
# -------------------------------------------------------------------------
my_func(country1)
# # A tibble: 6 x 5
# # Groups:   wave [1]
#   country  wave hi130 fixrate hi130_eur
#     <int> <int> <int>   <dbl>     <dbl>
# 1       2     1 84000       1     84000
# 2       2     1 20550       1     20550
# 3       2     1     0       1         0
# 4       2     1 64000       1     64000
# 5       2     1 64000       1     64000
# 6       2     1     0       1         0
#

Применить к списку объектов фрейма данных

Вот как вы можете применить функцию к списку объектов вашего фрейма данных. 

# Add your dataframe objects to the list below
countries <-list(country1, country2) 

# lapply and save the output to df_list
df_list <- lapply(countries, my_func)

Вы можете перезаписать отдельные объекты фрейма данных, которые есть в списке countries, используя df_list следующим образом. 

for(i in seq_along(df_list)) {
  assign(paste0("country", i), df_list[[i]])
}
# -------------------------------------------------------------------------
# country2
# # A tibble: 6 x 5
# # Groups:   wave [1]
#   country  wave  hi130 fixrate hi130_eur
#     <dbl> <dbl>  <dbl>   <dbl>     <dbl>
# 1       1     2      0       3        0 
# 2       1     2  59800       3    19933.
# 3       1     2  20440       3     6813.
# 4       1     2  19181       3     6394.
# 5       1     2 121213       3    40404.
# 6       1     2      0       3        0 
#

Вы можете проверить другие объекты данных;у всех теперь будет hi130_eur столбец. Возможно, существует эффективный способ решения этой проблемы, но это тот, который у меня есть на данный момент. Смотрите ?assign, чтобы почувствовать, что он делает.

Дайте мне знать, если вы все еще сталкиваетесь с проблемой.

...