메뉴 건너뛰기

웹에서 하는 R 통계


안녕하세요. 공단자료를 가지고  표본조사를 하고 있습니다. 
그중 제조업을 추출하려고 하는데  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


번호 제목 글쓴이 날짜 조회 수
공지 묻고답하기 게시판을 운영합니다 [1] cardiomoon 2015.04.20 911
595 질문드립니다 new qubic 2020.05.26 1
594 heatmap [2] new zkbi 2020.05.26 1
593 아주 기본적인 질문드립니다. [1] update jjunyon 2020.05.25 7
592 웹R expoPlot 관련 질문드립니다 [1] updatefile sungmas 2020.05.23 14
591 생존분석에 관해 질문드립니다. [3] update 투명한나비 2020.05.22 22
590 ROC multiple predictor model에서 cutoff 질문드립니다 [1] 박소현 2020.05.22 18
589 선생님, 안녕하십니까. ARCL 관련 질문 드립니다. chrik 2020.05.22 10
588 diagnostic test [1] 김날 2020.05.21 9
587 P for trend [1] update Ryan 2020.05.19 23
586 PSM (propensity score matching) 결과 해석 관련 질문 드립니다. nsstp 2020.05.18 147
585 ggpubr에 대해 질문드립니다. [2] 송이싸리버섯 2020.05.18 21
584 ggplot 그래프 질문입니다. [7] R입문기초 2020.05.17 25
583 PSM 에 대해서 문의 드립니다 [1] file jaljary 2020.05.14 40
582 정회원 web-R ver.4.0 survival graph [1] file 딸기공주 2020.05.13 32
581 propensity score matching 에서 no units matched 라고 나옵니다. [2] file nsstp 2020.05.13 28
580 정회원 web-R ver.4.0 survival graph [2] 딸기공주 2020.05.13 24
579 cutoff 질문 루아 2020.05.12 17
578 propensity score matching 질문드립니다. [5] nsstp 2020.05.12 40
» 표본설계시 subset 명령어 문의드립니다. [1] 네잎클로바 2020.05.09 22
576 Disconnected from the server [3] docmk81 2020.05.04 43