qrea$hl <- (qrea$hl-coef(hl.rlm)[1]) / coef(hl.rlm)[2]
qrea$hr <- (qrea$hr-coef(hr.rlm)[1]) / coef(hr.rlm)[2]
qrea$vl <- 2* (qrea$vl-mean(qrea$vl,trim=0.2,na.rm=T)) / coef(hl.rlm)[2]
qrea$vr <- 2* (qrea$vr-mean(qrea$vr,trim=0.2,na.rm=T)) / coef(hr.rlm)[2]