Передача текущего значения ddply разделения на функцию

Question

Вот некоторые выборочные данные, для которых я хочу, чтобы закодировать пол имен с течением времени:

names_to_encode <- structure(list(names = structure(c(2L, 2L, 1L, 1L, 3L, 3L), .Label = c("jane", "john", "madison"), class = "factor"), year = c(1890, 1990, 1890, 1990, 1890, 2012)), .Names = c("names", "year"), row.names = c(NA, -6L), class = "data.frame") 

Вот минимальный набор данных по социальному обеспечению, ограниченный только те имена, с 1890 по 1990 год:

ssa_demo <- structure(list(name = c("jane", "jane", "john", "john", "madison", "madison"), year = c(1890L, 1990L, 1890L, 1990L, 1890L, 1990L), female = c(372, 771, 56, 81, 0, 1407), male = c(0, 8, 8502, 29066, 14, 145)), class = c("tbl_df", "tbl", "data.frame"), row.names = c(NA, -6L), .Names = c("name", "year", "female", "male")) 

Я определил функцию, которая подмножество данных по социальному обеспечению, приведенным в году или диапазон лет. Другими словами, он вычисляет, было ли имя мужчиной или женщиной в течение определенного периода времени, выясняя пропорцию мужских и женских родов с этим именем. Вот функция, наряду с вспомогательной функцией:

require(plyr) 
require(dplyr) 

select_ssa <- function(years) { 

    # If we get only one year (1890) convert it to a range of years (1890-1890) 
    if (length(years) == 1) years <- c(years, years) 

    # Calculate the male and female proportions for the given range of years 
    ssa_select <- ssa_demo %.% 
    filter(year >= years[1], year <= years[2]) %.% 
    group_by(name) %.% 
    summarise(female = sum(female), 
       male = sum(male)) %.% 
    mutate(proportion_male = round((male/(male + female)), digits = 4), 
      proportion_female = round((female/(male + female)), digits = 4)) %.% 
    mutate(gender = sapply(proportion_female, male_or_female)) 

    return(ssa_select) 
} 

# Helper function to determine whether a name is male or female in a given year 
male_or_female <- function(proportion_female) { 
    if (proportion_female > 0.5) { 
    return("female") 
    } else if(proportion_female == 0.5000) { 
    return("either") 
    } else { 
    return("male") 
    } 
} 

Теперь то, что я хочу сделать, это использовать plyr, в частности ddply, чтобы подмножество данных, подлежащих закодированные в год, и объединить каждый из этих частей со значением возвращенного по функции select_ssa. Это код, который у меня есть.

ddply(names_to_encode, .(year), merge, y = select_ssa(year), by.x = "names", by.y = "name", all.x = TRUE) 

При вызове select_ssa(year), эта команда работает просто отлично, если я жесткий код значение как 1890 в качестве аргумента функции. Но когда я пытаюсь передать его текущее значение year, что ddply работает с, я получаю сообщение об ошибке:

Error in filter_impl(.data, dots(...), environment()) : 
    (list) object cannot be coerced to type 'integer' 

Как я могу передать текущее значение year на ddply?

Answer 1

Я думаю, что вы делаете что-то слишком сложное, пытаясь сделать соединение внутри ddply. Если бы я должен был использовать dplyr я бы, вероятно, сделать что-то больше, как это:

names_to_encode <- structure(list(name = structure(c(2L, 2L, 1L, 1L, 3L, 3L), .Label = c("jane", "john", "madison"), class = "factor"), year = c(1890, 1990, 1890, 1990, 1890, 2012)), .Names = c("name", "year"), row.names = c(NA, -6L), class = "data.frame") 

ssa_demo <- structure(list(name = c("jane", "jane", "john", "john", "madison", "madison"), year = c(1890L, 1990L, 1890L, 1990L, 1890L, 1990L), female = c(372, 771, 56, 81, 0, 1407), male = c(0, 8, 8502, 29066, 14, 145)), class = c("tbl_df", "tbl", "data.frame"), row.names = c(NA, -6L), .Names = c("name", "year", "female", "male")) 

names_to_encode$name <- as.character(names_to_encode$name) 
names_to_encode$year <- as.integer(names_to_encode$year) 

tmp <- left_join(ssa_demo,names_to_encode) %.% 
     group_by(year,name) %.% 
     summarise(female = sum(female), 
       male = sum(male)) %.% 
     mutate(proportion_male = round((male/(male + female)), digits = 4), 
      proportion_female = round((female/(male + female)), digits = 4)) %.% 
     mutate(gender = ifelse(proportion_female == 0.5,"either", 
         ifelse(proportion_female > 0.5,"female","male"))) 

Обратите внимание, что 0.1.1 еще немного привередливо о типах столбцов, так что я должен был преобразовать их. Я думаю, что я видел некоторую активность в github, которая предположила, что это было либо исправлено в dev-версии, либо, по крайней мере, над чем-то, над чем они работают.