search for: pnchisq

...*************************************************** *** src/library/base/man/NChisquare.orig Sat Jul 26 10:50:54 1997 --- src/library/base/man/NChisquare Sat Jul 26 11:26:23 1997 *************** *** 1,26 **** TITLE(NChisquare @@ The Non-Central Chi-Square Distribution) USAGE( ! dnchisq(x, df) ! pnchisq(x, df) ! qnchisq(x, df) ! rnchisq(x, df) ) ALIAS(dnchisq) ALIAS(pnchisq) ALIAS(qnchisq) ALIAS(rnchisq) DESCRIPTION( ! These functions provide information about the non-central ! chi-square distribution with LANG(df) degrees of freedom. ! LANG(dnchisq) gives the density, LANG(pnchisq) gi...

...<krcabrer@epm.net.co> writes: > Hello developers and users: > > My system fails (the computer freezes) when I use the ncp parameter, > with the lower.tail=FALSE option in the qchisq function. > > qchisq(0.025,31,ncp=1,lower.tail=FALSE) Yup, that's a bug. We have in pnchisq.c 48 for (ux = 1.0; pnchisq(ux, n, lambda, lower_tail, log_p) < p; ux *= 2); 49 for (lx = ux; pnchisq(lx, n, lambda, lower_tail, log_p) > p; lx *= 0.5); but if we look at the opposite tail, we also need to reverse that logic since pnchisq(x,...) is now decreasing in x. Ot...

...velopers and users: >> >> My system fails (the computer freezes) when I use the ncp parameter, >> with the lower.tail=FALSE option in the qchisq function. >> >> qchisq(0.025,31,ncp=1,lower.tail=FALSE) PD> Yup, that's a bug. We have in pnchisq.c PD> 48 for (ux = 1.0; pnchisq(ux, n, lambda, lower_tail, log_p) < p; PD> ux *= 2); PD> 49 for (lx = ux; pnchisq(lx, n, lambda, lower_tail, log_p) > p; PD> lx *= 0.5); PD> but if we look at the opposite tail, we also need to reverse that PD&gt...

Full_Name: Alan Bain Version: 2.4.0 OS: XP Submission from: (NULL) (155.140.122.227) Code for pnchisq contains following if (tSml) { if (x> f+ theta+ 5*sqrt(2*(f+2*theta) ))){ return 1.; /* better than 0 --- but definately FIXME */ } } This needs to check which tail has been requested; it is only correct if the default lower_tail=1 has been requested; for upper tail should return 0 if m...

Hi, Thanks, everyone for all the help! So, here is my calling function in C (called test.c): #include<stdio.h> #include<stdlib.h> #include<Rmath.h> int main(void) { printf("%f \n",pchisq(2.,7., 1, 0)); printf("%f \n",pnchisq(2.,7.,0., 1, 0)); return EXIT_SUCCESS; } I compile using: gcc test.c -I/usr/lib/R/include -L/usr/lib/R/lib -lm -lR However, running ./a.out gives me: 1.000000 0.040160 The first is wrong, but the second non-central is correct, and matches the answer from R. Inciden...

...nbinom.o dnchisq.o dnf.o dnorm.o dnt.o dpois.o dt.o dunif.o dweibull.o fmax2.o fmin2.o ftrunc.o gamma.o gamma_cody.o i1mach.o imax2.o imin2.o lbeta.o lgamma.o lgammacor.o mlutils.o pbeta.o pbinom.o pcauchy.o pchisq.o pexp.o pf.o pgamma.o pgeom.o phyper.o plnorm.o plogis.o pnbeta.o pnbinom.o pnchisq.o pnf.o pnmath.o pnorm.o pnt.o polygamma.o ppois.o pt.o ptukey.o punif.o pweibull.o qbeta.o qbinom.o qcauchy.o qchisq.o qexp.o qf.o qgamma.o qgeom.o qhyper.o qlnorm.o qlogis.o qnbeta.o qnbinom.o qnchisq.o qnf.o qnorm.o qnt.o qpois.o qt.o qtukey.o qunif.o qweibull.o sign.o stirlerr.o toms7...

...alled > > | test.c): > > | > > | #include<stdio.h> > > | #include<stdlib.h> > > | #include<Rmath.h> > > | > > | int main(void) { > > | printf("%f \n",pchisq(2.,7., 1, 0)); > > | printf("%f \n",pnchisq(2.,7.,0., 1, 0)); > > | return EXIT_SUCCESS; > > | } > > | > > | I compile using: > > | > > | gcc test.c -I/usr/lib/R/include -L/usr/lib/R/lib -lm -lR > > | > > | However, running > > | ./a.out > > | > > | gives me...

Full_Name: Drouilhet R?my Version: 1.8.1 OS: Linux Submission from: (NULL) (195.221.43.136) qchisq(1,10) works well but qchisq(1,10,ncp=0) does not work whereas ncp=0 is the default value of the function qchisq(1,10). (of course, 10 will be replaced by any integer value). Let us notice that this bug occurs only when applying probability one. (qchisq(seq(0,.9,.1),10,ncp=0) works very well).

Full_Name: Richard Mott Version: 1.9.0 OS: Windows XP Submission from: (NULL) (81.178.233.208) Shouldn't these give the same answer? > pchisq(67.60644,df=1,lower.tail=F,ncp=0) [1] 3.219647e-15 > pchisq(67.60644,df=1,lower.tail=F) [1] 1.996145e-16 >

...cauchy.o pcauchy.o qcauchy.o rcauchy.o dexp.o pexp.o qexp.o rexp.o dgeom.o pgeom.o qgeom.o rgeom.o dhyper.o phyper.o qhyper.o rhyper.o dnbinom.o pnbinom.o qnbinom.o rnbinom.o dpois.o ppois.o qpois.o rpois.o dweibull.o pweibull.o qweibull.o rweibull.o dlogis.o plogis.o qlogis.o rlogis.o dnchisq.o pnchisq.o qnchisq.o dnbeta.o pnbeta.o qnbeta.o pnf.o pnt.o qnf.o qnt.o ptukey.o qtukey.o toms708.o wilcox.o signrank.o ranlib libnmath.a config.status: creating src/unix/Makefile make: /usr/local/R-devel/srcunix: No such file or directory *** Error code 2 Stop in /usr/local/R-devel/src/unix. *** Error co...

There is a problem with `rpois'. It does seem to take care about the order of the arguments. This is an example: > rpois(n=1,lambda=2) [1] 3 > rpois(lambda=2,n=1) [1] 2 0 It obviously uses the first argument as the number of samples to be drawn, which is wrong. I used Version 0.49 Beta (April 23, 1997). Fredrik

There is a problem with `rpois'. It does seem to take care about the order of the arguments. This is an example: > rpois(n=1,lambda=2) [1] 3 > rpois(lambda=2,n=1) [1] 2 0 It obviously uses the first argument as the number of samples to be drawn, which is wrong. I used Version 0.49 Beta (April 23, 1997). Fredrik

I have dowloaded the Source Code of R,and I want to know the process of chi-sqared test,but how can I found it? [[alternative HTML version deleted]]

...../include -c qlogis.c -o qlogis.o gcc -g -I../include -c rlogis.c -o rlogis.o gcc -g -I../include -c dnbinom.c -o dnbinom.o gcc -g -I../include -c pnbinom.c -o pnbinom.o gcc -g -I../include -c qnbinom.c -o qnbinom.o gcc -g -I../include -c rnbinom.c -o rnbinom.o gcc -g -I../include -c pnchisq.c -o pnchisq.o gcc -g -I../include -c qnchisq.c -o qnchisq.o gcc -g -I../include -c dnorm.c -o dnorm.o gcc -g -I../include -c pnorm.c -o pnorm.o gcc -g -I../include -c qnorm.c -o qnorm.o gcc -g -I../include -c rnorm.c -o rnorm.o gcc -g -I../include -c dpois.c -o dpois.o gcc -g -I../incl...

Hi, I think I have identified a inaccuracy in dchisq when the non-centrality parameter is non-zero and large. Here's a little test: sys.dchisq.test <- function(N = 100000,mean = 0) { z <- rnorm(N,mean = mean, sd = 1) x <- z^2 xmin <- min(x) xmax <- max(x) br <- seq(xmin,xmax,length = 101) dbr <- br[2]-br[1] hist(x,br) p <- dchisq(br,df = 1,ncp =

Hi, I am wondering where the library file for R's nmath routines are? Doing a search on libR gave me the following: /usr/lib/libRKC16.so.1.2.0 /usr/lib/libRKC.so.1.2.0 /usr/lib/R/lib/libRlapack.so /usr/lib/R/lib/libR.so /usr/lib/libRKC.so.1 /usr/lib/libRKC16.so.1 None of these have the functions in nmath. Any help? Many thanks and best wishes! GT

...ct -g -I../include -I../../src/include -c qlogis.c -o qlogis.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c rlogis.c -o rlogis.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c dnchisq.c -o dnchisq.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c pnchisq.c -o pnchisq.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c qnchisq.c -o qnchisq.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c dnbeta.c -o dnbeta.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c pnbeta.c -o pnbeta.o gcc -ieee_with_inexact -g...

...lude -c dlogis.c cc -ieee_with_inexact -g -I../../src/include -c plogis.c cc -ieee_with_inexact -g -I../../src/include -c qlogis.c cc -ieee_with_inexact -g -I../../src/include -c rlogis.c cc -ieee_with_inexact -g -I../../src/include -c dnchisq.c cc -ieee_with_inexact -g -I../../src/include -c pnchisq.c cc -ieee_with_inexact -g -I../../src/include -c qnchisq.c cc -ieee_with_inexact -g -I../../src/include -c dnbeta.c cc -ieee_with_inexact -g -I../../src/include -c pnbeta.c cc -ieee_with_inexact -g -I../../src/include -c pnf.c cc -ieee_with_inexact -g -I../../src/include -c pnt.c cc -ieee_wit...

...nclude -I../../src/include -c qlogis.c -= o = qlogis.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c rlogis.c -= o = rlogis.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c dnchisq.c = -o = dnchisq.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c pnchisq.c = -o = pnchisq.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c qnchisq.c = -o = qnchisq.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c dnbeta.c -= o = dnbeta.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c pnbeta.c -= o = pnbeta.o gcc -ieee_...

...../../src/include -c plogis.c -o plogis.o gcc -g -O2 -I../include -I../../src/include -c qlogis.c -o qlogis.o gcc -g -O2 -I../include -I../../src/include -c rlogis.c -o rlogis.o gcc -g -O2 -I../include -I../../src/include -c dnchisq.c -o dnchisq.o gcc -g -O2 -I../include -I../../src/include -c pnchisq.c -o pnchisq.o gcc -g -O2 -I../include -I../../src/include -c qnchisq.c -o qnchisq.o gcc -g -O2 -I../include -I../../src/include -c dnbeta.c -o dnbeta.o gcc -g -O2 -I../include -I../../src/include -c pnbeta.c -o pnbeta.o gcc -g -O2 -I../include -I../../src/include -c pnf.c -o pnf.o gcc -g -O2 -...