Skip to menu

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

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


안녕하세요. 공단자료를 가지고  표본조사를 하고 있습니다. 
그중 제조업을 추출하려고 하는데  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
No. Subject Author Date
1153 webrPSM 설치 문제 [1] hawk1227 2023.01.03
1152 정회원용서버 ExploPlots 오류 [1] 버블 2023.01.03
1151 생존분석에서 세군이상 시 사후 검정이 가능할지 여쭙습니다 [1] 정망고 2023.01.03
1150 gaze에서 show.stats Arguments가 작동하지 않습니다. [1] 김찬식 2022.12.31
1149 서버 연결이 자꾸 끊깁니다. [1] secret wave 2022.12.30
1148 PSM, full matching, cox proportional hazard model [4] file 주니 2022.12.28
1147 비밀번호 찾기 오류. [2] ezra006 2022.12.24
1146 그래프 만들던 중 오류 질문드립니다.. [1] file 헝힝헝 2022.12.14
1145 multipleROC 패키지 질문 [2] 강우성 2022.12.14
1144 웹에서 하는 R통계 (무료버전) 업로드 가능한 파일 용량이 궁금합니다. [2] 성열 2022.12.10
1143 데이터프레임 자료를 ggplot2를 이용해 꺾은선그래프 그리는 방법 [1] 다다다다빈 2022.12.06
1142 서버에서 table 만들기 에러 [1] 버블 2022.12.04
1141 mytable 문의드립니다. [1] file seyi0918 2022.11.30
1140 교수님 안녕하세요 moonbook 패키지 mycsv 오류가 나서 문의드립니다. [1] holyhyun 2022.11.15
1139 autoReg(fit, uni=TRUE)시, multivariable에 들어갈 변수들을 manual로 선택할 수 있을까요? [6] 김찬식 2022.11.14
1138 processR 패키지에서 Measurment model 정의가 가능한지요. chrik 2022.11.13
1137 Incidence rate per person-years 및 rate ratio에 대한 95% CI 산출 방법 문의 [1] file 꾸이뭉 2022.11.11
1136 서열로짓 시 평균한계효과 요다 2022.10.30
1135 서버문제로 ppt 다운로드가 안된다고 뜹니다. [5] liverim 2022.10.26
1134 메타분석 문의드립니다. file neurojang 2022.10.23
Write Tag