Skip to menu

표본설계시 subset 명령어 문의드립니다.

네잎클로바 2020.05.09 08:32 Views : 9371


안녕하세요. 공단자료를 가지고  표본조사를 하고 있습니다. 
그중 제조업을 추출하려고 하는데  svy함수가 계속 에러가 납니다. 
제가 subset 함수부분이 잘못되어서 그럴까요? 도움을 부탁드립니다.
인자가 많은데요. 설문문항이 봐야할 부분이 7항( Q62_13,Q62_14,Q62_6, Q62_16, Q62_110, Q25_1_1, Q25_2_1) 정도 됩니다.
변수가 길어지는데 subset에 다 넣을수 있을까요?

Q62_13(요통), Q25_2_1(소음)
Q62_13<-recode(Q62_1_3, "1=1;2:9=2;NA=NA")
Q25_2_1<-recode(Q25_2, "1:5=1;6:9=2;NA=NA")

i <- subset(kwcs5th_190924, TINDUST==3, cut= c(Q62_13,Q25_1_1))
i <- svydesign(data=i, ids=kwcs5th_190924$ID, strata=NULL, variables = NULL, weights=kwcs5th_190924$WT2)


svytable(~Q62_13+Q25_1_1, design=i)
svyby(~Q62_13, ~Q25_1_1, design=i, svymean, na.rm=TRUE, data=kwcs5th_190924, level=0.95) 

콘솔창에서 이런 에러가 나옵니다.
subset 
 class(TINDUST)
[1] "numeric"
> i <- subset(kwcs5th_190924, TINDUST==3, cut= c(Q62_13,Q25_1_1))
> i <- svydesign(data=i, ids=kwcs5th_190924$ID, strata=NULL, variables = NULL, weights=kwcs5th_190924$WT2)
>  
> svytable(~Q62_13+Q25_1_1, design=i)
Error in model.frame.default(formula = weights ~ Q62_13 + Q25_1_1, data = model.frame(design)) : 
  변수의 길이들이 다릅니다 ('Q25_1_1'에서 찾았습니다)
> svyby(~Q62_13, ~Q25_1_1, design=i, svymean, na.rm=TRUE, data=kwcs5th_190924, level=0.95) 
Error in x$fpc$sampsize[i, , drop = FALSE] : 
  (subscript) logical subscript too long


I want to
List
Write Tag