Add ncp argument qt R package

[jeremyirilli]

library(intarith)
qt(...(0.2, 0.3), ...(10, 11))
-0.8790578...-0.5399379

qt(...(0.2, 0.3), ...(10, 11), 0)
-0.8790578...-0.5399379

qt(...(0.2, 0.3), ...(10, 11), 1)
0.1615167...0.4845898

qt(...(0.2, 0.3), ...(10, 11), 1, lower.tail=F)
1.589613...1.966734

qt(...(log(0.2), log(0.3)), ...(10, 11), 1, lower.tail=F, log.p=T)
1.589613...1.966734

[mgondan]

dt oder qt?

[jeremyirilli]

qt, dt gibt es schon

[mgondan]

bei dt bitte beachten, dass min max dort nicht funktioniert, hast Du das beachtet?

[jeremyirilli]

Ja, dafür hatte ich optimize_ definiert, um zuerst jeweils das Maximum zu finden und dann bei eval_min_max zu berücksichtigen.

dt1(A1...A2, Df1...Df2, Ncp1...Ncp2, Log, Res) :-
    optimize_(dt, [A1...A2, Df1...Df2, Ncp1...Ncp2, Log], true, Maxima),
    eval_min_max(dt, [A1...A2, Df1...Df2, Ncp1...Ncp2, log=Log], [mode(Maxima)], Res).