Displaying 1 result from an estimated 1 matches for "ptd3".
Did you mean:
pt3
1999 Mar 10
0
[R] bug? and New bug. --- patch for pt() only (PR#138)
...it's even harder to decide
WHEN to switch from the (Taylor kind) series to the asymptotic formula
]
This has to wait for a while, unless someone else....
BTW, S-plus also fouls up completely here:
df3 <- c(10^(3:20))
## completely wrong on S; taking forever on R
plot(df3, ptd3 <- pt(1.96,df=df3), log='x', type='b')
ptd3
giving
[1] 9.748634e-01 9.749882e-01 9.750007e-01 9.750020e-01
[5] 9.750021e-01 9.750021e-01 9.750021e-01 9.750023e-01
[9] 9.750054e-01 9.750436e-01 9.748710e-01 9.725377e-01
[13] 8.914032e-01 -4.261217e+06...