Hello, The function qt returns NaN for degrees of freedom <1. For example:> qt(0.5,0.5)[1] NaN Warning message: In qt(p, df, lower.tail, log.p) : NaNs produced But qt(0.5,0.5) should be 0, since the distribution is symmetric.> pt(0,0.5)[1] 0.5 It actually fails with any value, as long as df<1. Is this a bug, or is there some fundamental reason why this cannot be computed? Thank you, Enrico [[alternative HTML version deleted]]