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
1213 교수님 지난 번 문의 드렸던 PSM 재문의 드립니다. [3] secret 긍정적인마음가짐 2023.06.27
1212 안녕하세요 ^^ odds ratio 와 관련된 문의드리고자 합니다 ^^ [2] 성열 2023.06.26
1211 mytable() 함수에서 na처리와 pvalue 계산 문의드립니다. [1] hjk34 2023.06.23
1210 교수님 안녕하세요 ^^ web-R table 작성관련 문의드립니다. [2] file mkw7654 2023.06.21
1209 Methods for continuous variables과 P value corrected for multiple comparisons 방법 궁금합니다. [1] 호호셀 2023.06.19
1208 ordinal logistic regression analysis 문의 [1] 버블 2023.06.18
1207 e-mail 변경 문의 [1] secret kimdo 2023.06.16
1206 안녕하세요 number at risk 관련 문의 드립니다 [1] secret cebukim 2023.06.13
1205 Propensity matching score 에 대해 문의드립니다. [1] saga 2023.06.08
1204 콕스비례위험모델의 적용에 관련해서 질문 드립니다. [2] secret kng1032 2023.06.07
1203 기초통계 계산 [2] secret leesav 2023.06.07
1202 정회원 후원하였습니다. [2] secret toast 2023.06.04
1201 logistic regression analysis 가 에러가 뜹니다. [1] secret Ikarus86 2023.06.02
1200 PMS caliper 에 대해 여쭤봅니다. [2] secret pjhaaa 2023.05.30
1199 ggGam에 대하여 [2] leesav 2023.05.22
1198 R 에서 PFS 의 confidence interval 을 구하고 싶을 때 어떻게 하면 될까요? [1] ddeeny 2023.05.21
1197 메타분석 관련 문의드립니다 secret 화이팅 2023.05.12
1196 교수님, 새로운 앱 개발 문의드립니다! [4] swpapa 2023.05.08
1195 교수님, 안녕하세요? 전용 앱 관련 문의드립니다! [1] swpapa 2023.05.07
1194 PSM 시 범주형 변수가 총 10가지 일 때 어떻게 돌려야 하는지요 [1] 긍정적인마음가짐 2023.05.06
Write Tag