안녕하세요. 공단자료를 가지고 표본조사를 하고 있습니다.
그중 제조업을 추출하려고 하는데 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
Comment 1
-
cardiomoon
2020.05.10 09:27
No. | Subject | Author | Date |
---|---|---|---|
813 | 질문드립니다 | qubic | 2020.05.26 |
812 | heatmap [2] | zkbi | 2020.05.26 |
811 | 생존분석에 관해 질문드립니다. [3] | 투명한나비 | 2020.05.22 |
810 | ROC multiple predictor model에서 cutoff 질문드립니다 [1] | 박소현 | 2020.05.22 |
809 | 선생님, 안녕하십니까. ARCL 관련 질문 드립니다. | chrik | 2020.05.22 |
808 | diagnostic test [1] | 김날 | 2020.05.21 |
807 | P for trend [1] | Ryan | 2020.05.19 |
806 | PSM (propensity score matching) 결과 해석 관련 질문 드립니다. | nsstp | 2020.05.18 |
805 | ggpubr에 대해 질문드립니다. [2] | 송이싸리버섯 | 2020.05.18 |
804 | ggplot 그래프 질문입니다. [7] | R입문기초 | 2020.05.17 |
803 | 정회원 기간 확인 [3] | 김날 | 2020.05.14 |
802 | PSM 에 대해서 문의 드립니다 [1] | jaljary | 2020.05.14 |
801 | 교수님 안녕하십니까? 논문 작성에 대해 질문이 있습니다. [1] | 디프 | 2020.05.14 |
800 | 정회원 web-R ver.4.0 survival graph [1] | 딸기공주 | 2020.05.13 |
799 | propensity score matching 에서 no units matched 라고 나옵니다. [2] | nsstp | 2020.05.13 |
798 | Plots of simple regression 다운로드 [2] | 김밥친구 | 2020.05.13 |
797 | 정회원 web-R ver.4.0 survival graph [2] | 딸기공주 | 2020.05.13 |
796 | cutoff 질문 | 루아 | 2020.05.12 |
795 | propensity score matching 질문드립니다. [5] | nsstp | 2020.05.12 |
» | 표본설계시 subset 명령어 문의드립니다. [1] | 네잎클로바 | 2020.05.09 |
이런 경우 사용하신 데이터와 R코드로 minimal reproducible example을 같이 올려주셔야 해결 가능합니다.