Duplicate alternatives
Posted: Thu Nov 08, 2012 9:49 am
I have run the following design and the output looks good, except I am getting some duplicates. Specifically, all attribute levels for alternatives 2 and 3 (or 4 and 5) are identical in some rows. I do not think this should be happening. Can someone please advise? (I only show the code for the first segment because the others are almost identical and my problem does not have to do with pivoting). Many thanks!
Design
;alts(Luxury) = car1, car2, car3, car4, car5
;alts(Trucks) = car1, car2, car3, car4, car5
;alts(CarsMid) = car1, car2, car3, car4, car5
;alts(CarsEcon) = car1, car2, car3, car4, car5
;rows = 120
;eff = fish(mnl,d,mean)
;bdraws = halton(500)
;fisher(fish) = Module7(Luxury[.08], Trucks[.445], CarsMid[.152], CarsEcon[.323])
;model(Luxury):
U(car1)= a1*range_gas[300] + b1[(n,-.000142,.0000137)]*price.ref[78219]+ b3[(n,-.174,.0154)]*refuel_gas.ref[.234] + b33[(u,0,0)]*refuel_electric.ref[.234] /
U(car2)=b1*price.piv[5%,15%,25%,50%]+ b2ev[(n,.0263,.002287)]*range_ev[50,75,100,200]+ b3*refuel_electric.piv[-30%,-40%,-50%,-60%]+ b5*hov[0,1]+b6*range_ev*refuel_electric + b7[.001]*range_ev*refuel_gas + b8*range_ev*price/
U(car3)=b1*price.piv[5%,15%,25%,50%]+ b2ev[(n,.0263,.002287)]*range_ev[50,75,100,200]+b3*refuel_electric.piv[-30%,-40%,-50%,-60%]+ b5*hov[0,1]+b6*range_ev*refuel_electric + b7[.001]*range_ev*refuel_gas + b8*range_ev*price/
U(car4)= a1*range_gas[300] + b1*price.piv[5%,15%,25%,50%] + b2phev[(n,.0263,.002287)]*range_phev[10,20,40,60] + b3a*refuel_electric.piv[-30%,-40%,-50%,-60%] + b3b*refuel_gas.piv[0%,-10%,-20%,-30%] + b5*hov[0,1] + b6*range_phev*refuel_electric + b7[.001]*range_phev*refuel_gas + b8*range_phev*price/
U(car5)= a1*range_gas[300] + b1*price.piv[5%,15%,25%,50%]+ b2phev[(n,.0263,.002287)]*range_phev[10,20,40,60]+b3a*refuel_electric.piv[-30%,-40%,-50%,-60%] + b3b*refuel_gas.piv[0%,-10%,-20%,-30%] + b5*hov[0,1] + b6*range_phev*refuel_electric + b7[.001]*range_phev*refuel_gas + b8*range_phev*price
Design
;alts(Luxury) = car1, car2, car3, car4, car5
;alts(Trucks) = car1, car2, car3, car4, car5
;alts(CarsMid) = car1, car2, car3, car4, car5
;alts(CarsEcon) = car1, car2, car3, car4, car5
;rows = 120
;eff = fish(mnl,d,mean)
;bdraws = halton(500)
;fisher(fish) = Module7(Luxury[.08], Trucks[.445], CarsMid[.152], CarsEcon[.323])
;model(Luxury):
U(car1)= a1*range_gas[300] + b1[(n,-.000142,.0000137)]*price.ref[78219]+ b3[(n,-.174,.0154)]*refuel_gas.ref[.234] + b33[(u,0,0)]*refuel_electric.ref[.234] /
U(car2)=b1*price.piv[5%,15%,25%,50%]+ b2ev[(n,.0263,.002287)]*range_ev[50,75,100,200]+ b3*refuel_electric.piv[-30%,-40%,-50%,-60%]+ b5*hov[0,1]+b6*range_ev*refuel_electric + b7[.001]*range_ev*refuel_gas + b8*range_ev*price/
U(car3)=b1*price.piv[5%,15%,25%,50%]+ b2ev[(n,.0263,.002287)]*range_ev[50,75,100,200]+b3*refuel_electric.piv[-30%,-40%,-50%,-60%]+ b5*hov[0,1]+b6*range_ev*refuel_electric + b7[.001]*range_ev*refuel_gas + b8*range_ev*price/
U(car4)= a1*range_gas[300] + b1*price.piv[5%,15%,25%,50%] + b2phev[(n,.0263,.002287)]*range_phev[10,20,40,60] + b3a*refuel_electric.piv[-30%,-40%,-50%,-60%] + b3b*refuel_gas.piv[0%,-10%,-20%,-30%] + b5*hov[0,1] + b6*range_phev*refuel_electric + b7[.001]*range_phev*refuel_gas + b8*range_phev*price/
U(car5)= a1*range_gas[300] + b1*price.piv[5%,15%,25%,50%]+ b2phev[(n,.0263,.002287)]*range_phev[10,20,40,60]+b3a*refuel_electric.piv[-30%,-40%,-50%,-60%] + b3b*refuel_gas.piv[0%,-10%,-20%,-30%] + b5*hov[0,1] + b6*range_phev*refuel_electric + b7[.001]*range_phev*refuel_gas + b8*range_phev*price