やっぱり，よくわかってなかったんだね，久保君

> pvalue <- 2 * pt(-abs(tvalue), df.r) と計算しているのか． 予想どーり t 分布の計算してるんだけど， これって絶対値にマイナスをつける必要があったんだ …… 知らなかった．

pt(-abs(tvalue), df.r) は，tvalue の正負にかかわらず，というか，，，，具体的に，例えば tvalue = 2 であっても， pt(-2, dr.r) というのは，t 値が -2 以下となる確率ということ，その2倍は，結局 |tvalue|>2 である確率ということ。2*pt(abs(tvalue), df.r, lower.tail=FALSE) というのと同じ。