R/gs_cattle.R

In adatalab/dairyCattle: Package

#' gs_cattle
#'
#' read_cattle for google spread sheet file
#' @param path google spread sheet file.
#' @param drop.zero remove the dataset that milk yeild == 0.
#' @param add add some columns for additional analysis.
#' @keywords dairy cattle
#' @import janitor
#' @import dplyr
#' @import lubridate
#' @import readxl
#' @export
#' @examples
#' df <- googlesheets::gs_title("example") %>% googlesheets::gs_read()
#' gs_cattle(path = df, drop.zero = TRUE)

gs_cattle <- function(path, drop.zero=FALSE, add=FALSE) {

  df <- path

  df <- janitor::clean_names(df,case = "lower_camel")
  df[,c(4,6,7,8,26,30)] <- lapply(df[,c(4,6,7,8,26,30)],FUN = ymd)
  df$분만후첫수정일까지일수 <- as.numeric(df$분만후첫수정일까지일수)
  df$산차 <- as.numeric(df$산차)

  if(drop.zero==TRUE) {
    df <- dplyr::filter(df, 유량 > 0)
  }

  if(add==TRUE) {

    # 4% FCM (Gaines, 1928; unit = kg/d)
    df$fcm40 <- as.numeric(0.4*df$유량 + 15*df$유지율/100*df$유량)

    # 분만예정일
    df$분만예정일 <- ifelse(
      df$최종수정일자 + 280 > df$검정일,
      df$최종수정일자 + 280,
      NA)
    df$분만예정일 <- as.Date(df$분만예정일, origin = "1970-01-01")

    # parity
    df$parity <- ifelse(
      df$산차 > 1,
      "Multiple",
      ifelse(df$산차 == 1, "First", NA))

    df$level <- ifelse(
      df$누적착유일수 < 65,
      "early",
      ifelse(df$누적착유일수 > 65 & df$누적착유일수 < 210, "mid", "late"))

    df$level <- factor(df$level, levels = c("early","mid","late"))

  }

  # attr(df, type) <- "kor"
  return(df)
}