use data,clear
/*
 qhapipf a1 a2, ipf(l1) reg(1) qt(y)
 qhapipf a1 a2, ipf(l1) reg([l1+l1]) qt(y)
 qhapipf a1 a2, ipf(l1) reg([l1a+l1b]) qt(y)
 qhapipf a1 a2, ipf(l1) reg([l1*l1]) qt(y)
 qhapipf a1 a2, ipf(l1) reg([l1a*l1b]) qt(y)
*/
 qhapipf a1 a2, ipf(l1) reg([l1+l1]) model(0) qt(y)
 qhapipf a1 a2, ipf(l1) reg(1) model(1) lrtest(0,1) qt(y)
 /*
qhapipf a1 a2 b1 b2, ipf(l1*l2) reg([l1+l1]+[l2+l2]) qt(y)
 qhapipf a1 a2 b1 b2, ipf(l1*l2) reg([l1+l1]*[l2+l2]) qt(y)
 qhapipf a1 a2 b1 b2, ipf(l1*l2) reg([l1*l1]*[l2*l2]) qt(y)
 qhapipf a1 a2 b1 b2, ipf(l1*l2) reg([l1a*l1b]*[l2a*l2b]) qt(y)
*/