similar to: R stat functions do not work as stated on the mannual (PR#7419)

Dear R Developers: I have been playing with R, release 2.0.1 for a week now and have detected = that all stat functions related to distribution probabilities have the same= problem: 1.- The log.p parameter of all distribution functions, when set to TRUE, re= turns a extrange value. Accoding to the manual, when set to true, it should return log(p) probabili= ty. So to my understanding, setting

> From: casadoj@ecc.es > > > Dear R Developers: > > I have been playing with R, release 2.0.1 for a week now and > have detected = > that all stat functions related to distribution probabilities > have the same= > problem: > > 1.- According to the manual the log.p parameter is always the > last one. > 2.- When you use the software, the last

Dear Sir I have successfuly downloaded R package. Kindly let me know where i will get _/*user mannual*/_ for the same. awaiting a prompt reply with regards Pankaj @#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@ Pankaj Kumar Research Scholar Climatology & Hydrometeorology Divison Indian Institute Of Tropical Meteorology Homi Bhabha Road, Pune-411008 India Phone No. +91-20-5893600 Ext. 361

Hola a todos Estoy cargando una base de datos a R llamada "Base", desde excel. A la hora de utilizar funciones como: >pairs(Base) Error in pairs.default(Base) : non-numeric argument to 'pairs' Una forma de arreglar este problema es utilizando la funcion as.numeric(), pero me toca hacerlo variable por variable: >Base$Variable1<-as.numeric( Base$Variable1)

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 <<insert bug report here>> Reproduced on Debian and Windows ... On 2.4.x if you execute set.seed(12345) t.1 <- rt(n = 1000, df = 20) set.seed(12345) t.2 <- rt(n = 1000, df = 20, ncp = 0) all.equal(t.1, t.2) ## Not close to true This appears to be due to the fact that in 2.4.x rt is now rt function (n, df, ncp = 0) { if

help(qt) states that: "ncp non-centrality parameter delta; currently except for rt(), only for abs(ncp) <= 37.62" so I would expect that calling qt with non-centrality parameter exceeding 37.62 should fail, instead e.g. calling > mapply(function(x) qt(p = 0.9, df = 55, ncp = x),35:45) gives: [1] 40.21448 41.35293 42.49164 43.68862 44.82945 45.97048 47.11170 48.25310 [9]

The function 'pchisq' from the 'stats' library gives a wrong result if the argument equals exactly zero: # Upper tail of central 1-df chi^2 distribution > pchisq(1 , 1, ncp=0, lower.tail = F, log.p = FALSE) [1] 0.3173105 > pchisq(0.5 , 1, ncp=0, lower.tail = F, log.p = FALSE) [1] 0.4795001 > pchisq(0.01 , 1, ncp=0, lower.tail = F, log.p = FALSE) [1]

Bug Description: Problem with the function pf() when the non-centrality parameter is large. Here is a sample command. You should see a smooth line from 0 to about 55, and then the values of pf() go crazy from 55 to 100. ############################ ncp <- seq(0,100,length=200) plot(ncp,pf(5,7,2,ncp=ncp)) ############################ Version: platform = i686-pc-linux-gnu arch = i686 os =

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

I cannot find the solutions in the mannual about "import and export of data". For example, I run >chisq.test(x,p=probs) chi-squared..... ......... How can I write the output into a file? Is there some file operations in R? Actually, I only want to save the P-value. But I don''t know how I can extract the P-value field only from the output of chisq.test. So, I think I

Dear all, I want to see the source codes for "dchisq(x, df, ncp=0, log = FALSE)", but cannot find it. I input "dchisq" in the R interface, and then enter, the following message return: > dchisq /*****************************************************/ function (x, df, ncp = 0, log = FALSE) { if (missing(ncp)) .Internal(dchisq(x, df, log)) else

Dear List: The square of the noncentral t-statistic with noncentrality parameter \delta is a noncentral F with noncentrality parameter \lambda=\delta^2. So, t^2_{\nu,\delta} = F_{1,\nu,\lambda=\delta^2}. Consequently, it should follow that t^2_{1-\alpha/2,\nu,\delta} = f_{1-alpha,1,\vu,\lambda=\delta^2}. However, this is not what is happening with the following code. The central

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)

Has any advancement been made to get 3102 operational in either a SIP or H323 asterisk environment. A post back in time mentioned a downloader service. >From the posts and articles I have read, the NCP is acting like a bootp and tftp server which uploads the configuration to the phone?? Am I close? if so, where does one get the SIP image for he 3102 and 2102 phones? I had 8 donated, but they

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 >

Dear people, I've a problem in doing a power calculation. In Fryer and Nicholson (1993), ICES J. mar. Sci. 50: 161-168 page 164 an example is given with the following characteristics T=5, points in time R=5, replicates Var.within=0.1 q=10, a 10% increase per year The degrees of freedom for the test are calculated as Vl=T*R-2=23 and the non-centrality parameter Dl=4.54. Using this they get a

Full_Name: David Clayton Version: 1.8.1 OS: Linux Submission from: (NULL) (131.111.126.242) qchisq behaves very strangely when ncp is passed as zero (forcing internal qnchisq to be called) when first argument is small. Eg > qchisq(1-1e-6, 1, ncp=0, lower.tail=TRUE) qchisq(1-1e-6, 1, ncp=0, lower.tail=TRUE) [1] 1024 while, if ncp is unspecified, > qchisq(1-1e-6, 1) qchisq(1-1e-6, 1)

Hello Everybody, Is there a function available for calculating the density of a non-central t distribution? (dt does not accept the ncp option). Thanks in advance! -- Wolfgang Viechtbauer "Are you not thinking what I am not thinking?" -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

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).