> with(qdata,tapply(response,IND=list(gram),mean))
     gram    ungram 
0.6862745 0.7058824 
> with(qdata,tapply(response,IND=list(int),mean))
       hi        lo 
0.6699346 0.7222222 

NP3:
g            0.02832    0.03771    0.75
i            0.02887    0.03790    0.76
gxi          0.08473    0.03755    2.26
V3:
g            0.02832    0.03771    0.75
i            0.02887    0.03790    0.76
gxi          0.08473    0.03755    2.26
V1:
g           -0.09896    0.03684   -2.69
i           -0.10995    0.03678   -2.99
gxi          0.08941    0.03680    2.43
Post-V1:
g           -0.08354    0.05377   -1.55
i           -0.05829    0.05346   -1.09
gxi          0.01961    0.05413    0.36


i            0.02887    0.03790    0.76  -0.04768  0.10581
i            0.02887    0.03790    0.76	 -0.04790  0.10637
i           -0.10995    0.03678   -2.99	 -0.18514 -0.03476
i           -0.05829    0.05346   -1.09	 -0.16730  0.05092

g           -0.09896    0.03684   -2.69  -0.17386 -0.02384
g           -0.08354    0.05377   -1.55	 -0.19153  0.02655


\begin{tabular}{rrrr}
  \hline
 & V1 & lower & upper \\
  \hline
  V1      & 0.09896 &  0.03684  & 2.69 & 0.01985 & 0.16949 \\
  Post-V1 & 0.08354 &  0.05377  & 1.55 & $-$0.02945 & 0.18896 \\
   \hline
\end{tabular}