nwos_mean: NWOS Mean
In bbutler01/nwos: NWOS Weights and Estimates

Description Usage Arguments Details Value References Examples

Calculate means for the NWOS.

1	nwos_mean(weight, area = 1, domain = 1, variable)

`weight`	vector of weights.
`area`	vector of area (e.g., forest acres in the stratum). Default = 1 (i.e., estimates are in terms of ownerships).
`domain`	binary vector (i.e., dummay variable) with 1 indicating inclusion in the domain of interest and 0 otherwise. Default = 1 (i.e., all ownerships are included).
`variable`	vector of values for the variable of interest.

This function needs to be run by stratum (e.g., family forest ownerships in a state).

Mean of variable of interest.

Butler, B.J. In review. Weighting for the US Forest Service, National Woodland Owner Survey. U.S. Department of Agriculture, Forest Service, Northern Research Station. Newotwn Square, PA.

wi <- tbl_df(read.csv("data/wi.csv")) %>% mutate(ROW_NAME = row.names(wi), AC_WOOD = ACRES_FOREST, FFO = if_else(LAND_USE == 1 & OWN_CD == 45 & AC_WOOD >= 1, 1, 0), RESPONSE = if_else(RESPONSE_PROPENSITY >= 0.5, 1, 0), RESPONSE = if_else(is.na(RESPONSE_PROPENSITY), 0, RESPONSE))
WI_FFO_AREA <- nwos_stratum_area(stratum = wi$FFO, point.count = wi$POINT_COUNT, state.area = 33898733)
WI_FFO_RR <- nwos_response_rate(stratum = wi$FFO, point.count = wi$POINT_COUNT, response = wi$RESPONSE)
wi$WEIGHT <- nwos_weights(stratum = wi$FFO, point.count = wi$POINT_COUNT, response = wi$RESPONSE, area = wi$AC_WOOD, stratum.area = WI_FFO_AREA, response.rate = WI_FFO_RR)
nwos_mean(weight = wi$WEIGHT, variable = wi$Y_3)
nwos_mean(weight = wi$WEIGHT, area = wi$AC_WOOD, variable = wi$Y_3)