How to include missing values in svytable in R to match Stata’s tab command?
Sure, here's a Stack Overflow question you could use:
I'm trying to replicate a Stata command in R using the survey package but I'm not getting the same results. In Stata, I use the following command to get weighted counts including missing values:
tab fp_message_radio [iweight=wt], missing
In R, I have created the survey design object and used svytable to get weighted counts for fp_message_radio:
# Create the survey design object
DHSdesign <- svydesign(ids = ~1, data = IRdata, weights = ~wt)
# Use svytable to get weighted counts for fp_message_radio
svytable(~fp_message_radio, DHSdesign)
However, this does not include the missing values. How can I modify my R code to include missing values and make it identical to the Stata command?
Solution
from the ?svytable help page, it looks like ... parameters get passed to the xtabs function, so would the addNA parameter from that function produce the results you expect? i've also included two lines to show the overall total.. thanks!!
library(survey)
data(api)
x <- apiclus1
x[ 1 , 'stype' ] <- NA
x[ 2 , 'sch.wide' ] <- NA
dclus1<-svydesign(id=~dnum, weights=~pw, data=x, fpc=~fpc)
svytable(~sch.wide+stype, dclus1)
svytable(~sch.wide+stype, dclus1,addNA=TRUE)
dclus1 <- update( dclus1 , one = 1 )
svytable( ~ one , dclus1 )
- summing across multiple columns named with a pattern with a condition on another column
- Reference current data.frame in dplyr flow
- How can I use dplyr summarize to create a count of all items and fraction of items in a certain group?
- Split one table to multiple tables
- Can I run an SQL update statement using only dplyr syntax in R
- iteratively assign values with case_when and for loops
- Why does a given attribute disappears after using dplyr::mutate?
- convenient way to select columns within a function based on regex within dplyr
- Adding column if it does not exist
- using apply and filter to R dataframe
- Use a for loop for a mutate function with an operation
- Alternatives to if_all and if_any for Data Manipulation in R
- How to select last N observation from each group in dplyr dataframe?
- data.table's equivalent of dplyr's bang bang (!!)
- Finding the most recent update within a table
- How to include missing values in svytable in R to match Stata’s tab command?
- How do I drop duplicates within groups that were collected within 60min of one another?
- How to avoid transposition of duplicates into lists with pivot_wider?
- R dplyr function ceases to work when using map to iterate over vector of column names
- in deply , how to select variable which contains string [cat] or [b]?
- How to fill NA's in a column if the column contains a number but leave NA if a column is all NAs
- Error: Data source must be a dictionary (dplyr)
- Is there a way to rearrange columns uniformly between multiple tables and plot?
- R: How to average a column between specific values in another?
- sf: Find bounding boxes for groups and comparison with ArcGIS Pro results
- dplyr::case_when throws "no non-missing arguments" with na.rm=TRUE even when compensating for all NA cases
- R: How to calculate cumulative mean over first N rows and fill result down column?
- Elongate data but also keep one of the original variable in R
- How to use if statements inside of groups in a dataframe in R?
- Collapsing rows by user with dplyr