Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- id = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12)
- strata = c(6, 6, 6, 7, 7, 7, 8, 8, 8, 8, 8, 8)
- weight = c(60, 75, 85, 140, 170, 175, 270, 310, 325, 785, 1450, 3920)
- fpc = c(8, 8, 8, 7, 7, 7, 6, 6, 6, 6, 6, 6)
- answer = c("2", "2", "3", "1", "2", NA, NA, 2, "3", NA, "1", NA)
- df = data.frame(id, strata, weight, fpc, answer)
- df <- df[complete.cases(df), ]
- dstrat<-svydesign(id=~1,strata=~strata, weights=~weight, data=df, fpc=~fpc)
- svymean(~answer, dstrat)
- mean SE
- answer1 0.60803 0.2573
- answer2 0.23518 0.1755
- answer3 0.15679 0.1479
- answer1 = good
- answer2 = neutral
- answer3 = bad
- library(survey)
- id <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12)
- strata <- c(6, 6, 6, 7, 7, 7, 8, 8, 8, 8, 8, 8)
- weight <- c(60, 75, 85, 140, 170, 175, 270, 310, 325, 785, 1450, 3920)
- fpc <- c(8, 8, 8, 7, 7, 7, 6, 6, 6, 6, 6, 6)
- answer <- c("2", "2", "3", "1", "2", NA, NA, 2, "3", NA, "1", NA)
- df <- data.frame(id=id, strata=strata, weight=weight, fpc=fpc, answer=answer)
- # this is probably a mistake
- df <- df[complete.cases(df), ]
- # in most data sets, you should be using na.rm=TRUE later
- # and not tossing out statements before the `svydesign` gets run
- df$ones <- as.numeric( df$answer %in% 1 )
- df$threes <- as.numeric( df$answer %in% 3 )
- dstrat<-svydesign(id=~1,strata=~strata, weights=~weight, data=df, fpc=~fpc)
- a <- svymean( ~ ones + threes , dstrat , na.rm = TRUE )
- svycontrast(a, list(avg=c(0,0), diff=c(1,-1)))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement