>>>>> "MW" == Morten Welinder <[EMAIL PROTECTED]> >>>>> on Tue, 24 Jul 2007 17:18:09 +0200 (CEST) writes:
MW> Full_Name: Morten Welinder MW> Version: 2.2.0 MW> OS: Linux MW> Submission from: (NULL) (216.223.241.229) MW> The qt function for df=1 is implemented as... MW> q = - tan((P+1) * M_PI_2); MW> Adding 1 kills accuracy for P near 0. A better way is to use MW> q = cot(P * M_PI_2); MW> or MW> q = 1/tan(P * M_PI_2); Yes, indeed, thanks a lot, Morten!! ((BTW: you should really upgrade your version of R)) Best regards, Martin ______________________________________________ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel