multiH<- glm(H ~sex+ region+ age+ incm+ edu+ occp + marri_1 + BD1_11+ BS3_1 +HE_BMI ,family=binomial(), na.action="na.exclude",data=hn1)
# sex with 1,2/ region for 1~10/ age for 30~99/ incm for 1~4/ occp for 1~3/ marri_1 for 1,2/ BD1_11 for 1,2/ BS3_1 for 1,2/ HE_BMI for 0~99
我的数据中没有NA。
还为某些变量使用了因子代码。
我不知道解决此错误问题的方法。
哪个是
#variable lengths differ for "...."
也许数据框应该相同，但不是......？
任何帮助将不胜感激:)
Error in model.frame.default(formula = H ~ sex + region + age + incm + :
variable lengths differ (found for 'sex')
> str(hn1)
'data.frame': 27887 obs. of 1 variable:
$ sex.age.incm.edu.occp.marri_1.BD1_11.BS3_1.H.Var: Factor w/ 19847 levels "1,0,0,0,0,0,0,0,0,0",..: 5222 13979 585 10962 13430 2938 12857 6137 15667 66 ...
这是我得到的错误标志......
>hn1<-read.csv("C:\\Users\\hayon14\\Desktop\\HN.csv", header=T, dec=".",sep="\t"
> H<-factor(hn1$H)
> region<-factor(hn1$region)
> sex<-factor(hn1$sex)
> age<-factor(hn1$age)
> incm<-factor(hn1$incm)
> edu<-factor(hn1$edu)
> occp<-factor(hn1$occp)
> marri_1<-factor(hn1$marri_1)
> BD1_11<-factor(hn1$BD1_11)
> BS3_1<-factor(hn1$BS3_1)
> HE_BMI<-factor(hn1$HE_BMI)
> multiH<- glm(H ~sex+ region+ age+ incm+ edu+ occp + marri_1 + BD1_11+ BS3_1 +HE_BMI ,family=binomial(), na.action="na.exclude",data=hn1)
这是我的完整代码。