dat<-read.csv("",header=TRUE)

V_func<-function(commute_time,commute_time_var,age,income,work_time){
 commute_time_ <- -0.000167
 commute_time_var_ <- -9.38*(10**-8)
 age_ <- 0.0809
 work_time_ <- 7.42*(10**-6)
 income_ <-0.00289
 V <- commute_time_*commute_time + commute_time_var_*commute_time_var + age_*age + income_*income + work_time_*work_time
 tau1 <- 4.05
 tau2 <- 5.32
 tau3 <- 6.28
 tau4 <- 7.49
 if(V<tau1){
  outp <- 0
 }elif(V<tau2){
  outp <- 1
 }elif(V<tau3){
  outp <- 2
 }elif(V<tau4){
  outp <- 3
 }else{
  outp <- 4
 }
 return(outp)
}

cnt0<-c(0,0,0,0,0)
for(i in range(1:length(dat))){
 subs<-V_func(dat$commute_time[i],dat$commute_time_var[i],dat$age[i],dat$income[i],dat$work_time[i])
 cnt0[subs] <- cnt0[subs] + 1
}


cnt1<-c(0,0,0,0,0)
for(i in range(1:length(dat))){
 subs<-V_func(dat$commute_time[i],dat$commute_time_var[i],dat$age[i],dat$income[i],dat$work_time[i]*0.8)
 cnt1[subs] <- cnt1[subs] + 1
}


cnt2<-c(0,0,0,0,0)
for(i in range(1:length(dat))){
 subs<-V_func(dat$commute_time[i],dat$commute_time_var[i],dat$age[i],dat$income[i],dat$work_time[i]*1.2)
 cnt2[subs] <- cnt2[subs] + 1
}
