Trip_hh <-merge(trip, hh, by = "HH_ID", all.x = FALSE, all.y=FALSE, sort = FALSE)
Trip_hh<-Trip_hh[which(Trip_hh$INCOME!=98&Trip_hh$INCOME!=99),]
Mode_Inc<-matrix(,ncol=10,nrow=7)
for(i in 1:10)
{
for (j in 1:7)
{
Mode_Inc[i,j]<- sum(Trip_hh$HH_WEIGHT.x[Trip_hh$INCOME==i&Trip_hh$MODE_AGG==j]/sum(Trip_hh$HH_WEIGHT.x[Trip_hh$INCOME==i]))
}
}
给我
[<-
中的错误(*tmp*
,i,j,值= sum(Trip_hh $ HH_WEIGHT.x [Trip_hh $ INCOME ==:
下标超出范围
期望的结果:10 * 7矩阵,当收入== 1和模式== 1(等等)时计算可交换权重。