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
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
|
#' Calculate Estimated Population
#'
#' @description \code{aw_calculate} multiplies the given \code{value} by the area weight. This
#' is the fourth step in the interpolation process after \link{aw_weight}.
#'
#' @usage aw_calculate(.data, value, areaWeight, newVar)
#'
#' @param .data A given intersected dataset
#' @param value A column within \code{source} to be interpolated
#' @param areaWeight The name of the variable containing area weight per feature
#' @param newVar Optional; a new field name to store the interpolated value in. If not specified,
#' the \code{value} argument will be used as the new field name.
#'
#' @return An intersected file of class sf with a new field of interest recalculated with area weight
#'
#' @examples
#' library(dplyr)
#'
#' race <- select(ar_stl_race, GEOID, TOTAL_E)
#' wards <- select(ar_stl_wards, WARD)
#'
#' wards %>%
#' aw_intersect(source = race, areaVar = "area") %>%
#' aw_total(source = race, id = GEOID, areaVar = "area", totalVar = "totalArea",
#' weight = "sum", type = "extensive") %>%
#' aw_weight(areaVar = "area", totalVar = "totalArea", areaWeight = "areaWeight") -> intersect
#'
#' aw_calculate(intersect, value = "TOTAL_E", areaWeight = "areaWeight")
#'
#' @importFrom dplyr mutate
#' @importFrom glue glue
#' @importFrom rlang :=
#' @importFrom rlang enquo
#' @importFrom rlang quo
#' @importFrom rlang quo_name
#' @importFrom rlang sym
#'
#' @export
aw_calculate <- function(.data, value, areaWeight, newVar){
# save parameters to list
paramList <- as.list(match.call())
# check for missing parameters
if (missing(.data)) {
stop("A sf object containing intersected data must be specified for the '.data' argument.")
}
if (missing(value)) {
stop("A variable name must be specified for the 'value' argument.")
}
if (missing(areaWeight)) {
stop("A variable name must be specified for the 'areaWeight' argument.")
}
# nse
if (!is.character(paramList$areaWeight)) {
areaWeightQ <- rlang::enquo(areaWeight)
} else if (is.character(paramList$areaWeight)) {
areaWeightQ <- rlang::quo(!! rlang::sym(areaWeight))
}
areaWeightQN <- rlang::quo_name(rlang::enquo(areaWeight))
if (!is.character(paramList$value)) {
valsQ <- rlang::enquo(value)
} else if (is.character(paramList$value)) {
valsQ <- rlang::quo(!! rlang::sym(value))
}
valsQN <- rlang::quo_name(rlang::enquo(value))
if (missing(newVar)) {
newFieldQN <- valsQN
} else if (!missing(newVar)){
newFieldQN <- rlang::quo_name(rlang::enquo(newVar))
}
# check variables
if(!!valsQN %in% colnames(.data) == FALSE) {
stop(glue::glue("Variable '{var}', given for the value, cannot be found in the given intersected object.",
var = valsQN))
}
if (!!areaWeightQN != "...areaWeight"){
if(!!areaWeightQN %in% colnames(.data) == FALSE) {
stop(glue::glue("Variable '{var}', given for the area weight, cannot be found in the given intersected object.",
var = areaWeightQN))
}
}
# recalculate source values of interest using area weight and assign as new field
out <- dplyr::mutate(.data, !!newFieldQN := !!valsQ * !!areaWeightQ)
# return output
return(out)
}
|
