similar to: Bug in 'pchisq' for x=0.0 (PR#9485)

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 >

hello R-friends, I'm looking for a quantile function for the noncentral f-distribution in the area of equivalence hypotheses testing. Can somebody help me? Many thanks ----------------------------------------------------------------- Dipl. Inform. J. Hedderich Institut f?r Medizinische Informatik Phone : 0431 / 5973182 und Statistik im Klinikum an der CAU

Dear R experts, I'm trying to call 'pchisq' from within a C subroutine. The following error is returned: ** NON-convergence in pgamma()'s pd_lower_cf() f= nan. This error message is not printed the first time I call 'pchisq' from the C subroutine, but the second time or the next time I call 'pchisq' from within R. My session output is shown below:

Full_Name: Jerry W. Lewis Version: 2.6.1 OS: Windows XP Professional Submission from: (NULL) (24.147.191.250) pchisq(0,0,ncp=lambda) returns 0 instead of exp(-lambda/2) pchisq(x,0,ncp=lambda) returns NaN instead of exp(-lambda/2)*(1 + SUM_{r=0}^infty ((lambda/2)^r / r!) pchisq(x, df + 2r)) qchisq(.7,0,ncp=1) returns 1.712252 instead of 0.701297103 qchisq(exp(-1/2),0,ncp=1) returns 1.238938

Dear R users, I am searching for a general minimisaton function like 'ms' in R. Please give me a hint! -- ------------------------------------------------------------ Dipl. Inform. J. Hedderich Institut f?r Medizinische Informatik und Statistik im Klinikum an der CAU Phone : 0431/5973182 Brunswiker Str. 10

This is just a suggestion/wish that it would be nice for the F-distribution functions to recognize limiting cases for infinite degrees of freedom, as the t-distribution functions already do. The t-distribution functions recognize that df=Inf is equivalent to the standard normal distribution: > pt(1,df=Inf) [1] 0.8413447 > pnorm(1) [1] 0.8413447 On the other hand, pf() will accept Inf

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]]

Dear all Just for fun, I have just downloaded the paper mentioned below and checked it with R-1.6.1. Everything is ok with exception of Table 2b, where I get always 1 instead of 0.5: > pbinom(1e15,2e15,0.5) [1] 1 Which value should be correct? Best regards Christian Stratowa ============================================== Christian Stratowa, PhD Boehringer Ingelheim Austria Dept NCE Lead

Dear colleague, I not sure this R code is correctly ? I would to show the number of Sample Size at Sample Size Axis that line draw from Power Axis (80%) from R code. How I show this and select the most appropriate of this power (.79955687 - 80983575). Thank for your help and answer. Best Regards, Nikom Thanomsieng, Email: nikom at kku.ac.th .... #Power analysis: Sample size for

for example, an expression such as chisq(df=1,ncp=0) ? thanks -- View this message in context: http://www.nabble.com/how-to-calculate-chisq-value-in-R-tp15338943p15338943.html Sent from the R help mailing list archive at Nabble.com.

Hello all, Here's my query: Running R 1.6.2 on FreeBSD 5.0, and on WinXP, and I find that the following hangs the process: dchisq(alpha=0.01, df=1, ncp=295) but it does work for ncp < 294.92. Is this general? Best wishes to all, Andrew Andrew Robinson Ph: 208 885 7115 Department of Forest Resources Fa: 208 885 6226 University of Idaho E : andrewr at uidaho.edu PO

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

Dear List, If any of observed and/or expected data has less than 5 frequencies, then chisq.test (Pearson's Chi-squared Test for Count Data from package:stats) gives warning messages. For example, x<-c(10, 14, 10, 11, 11, 7, 8, 4, 1, 4, 4, 2, 1, 1, 2, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1) y<-c(9.13112391745095, 13.1626482033341, 12.6623267638188, 11.0130706413029, 9.16415925139016,

I compiled the 1.9.1 src.rpm with the standard gnu tools and it works. I tried compiling the 1.9.1 src.rpm with the Intel 8 C and FORTRAN compilers and it bombs out during the testing phase: comparing 'd-p-q-r-tests.Rout' to './d-p-q-r-tests.Rout.save' ...267c267 < df = 0.5[1] "Mean relative difference: 5.001647e-10" --- > df = 0.5[1] TRUE

Hi Folks, I note that, while the "chisq" functions dchisq(x, df, ncp=0, log = FALSE) pchisq(q, df, ncp=0, lower.tail = TRUE, log.p = FALSE) qchisq(p, df, ncp=0, lower.tail = TRUE, log.p = FALSE) rchisq(n, df, ncp=0) all have a slot for the non-centrality parameter "ncp", of the functions for the t and F distributions: dt(x, df, log = FALSE)

Hi all, I would like to estimate power and necessary sample size for a GLM with a response variable that has a poisson distribution. Do you have any suggestions for how I can do this in R? Thank you for your help. Sincerely, Craig -- Craig A. Faulhaber Department of Forest, Range, and Wildlife Sciences Utah State University 5230 Old Main Hill Logan, UT 84322 (435)797-3892

Dear List, I will appreciate any advice regarding how to convert the following numbers [I got in return by taxondive()] in numeric integers without the e.g. 6.4836e+01 abbreviations. Thank you very much in advance, Gian > taxa_dive Species Delta Delta* Lambda+ Delta+ S Delta+ Nat1 5.0000e+00 6.4836e+01 9.5412e+01 6.7753e+02 8.7398e+01 436.99 Nat2

Full_Name: Viktor Witkovsky Version: 2.9.2 OS: Windows XP Submission from: (NULL) (78.98.89.227) Hello, I have found strange behavior of the function qchisq (the non-central qchisq is based on inversion of pchisq, which is further based on pgamma). The function gives wrong results without any warning. For example: qchisq(1e-12,1,8.94^2,lower.tail=FALSE) gives 255.1840972465858 (notice that

Hi! I got a rather stupid question (I think): Is there ANY option that makes R display numericals not like "1e-8" but as "0.00000001" by default ? And I need the outcome to be really numerical, so formatC(...) which produces a character or something like this won't be acceptable. Any help on this would be appreciated, thanx. Marc

Hi, directory<- "/home/arunksa111/data.new" #first function filelist<-function(directory,number,list1){ setwd(directory) filelist1<-dir(directory) direct<-dir(directory,pattern = paste("MSMS_",number,"PepInfo.txt",sep=""), full.names = FALSE, recursive = TRUE) list1<-lapply(direct, function(x) read.table(x,header=TRUE, sep =