similar to: exact test for large tables. Was: unexpected R crash - again

This came up on r-help but indicates a bug. dnbinom(x,n,p) calls dbinom_raw(n-1,...) which returns 0 for n<1. -thomas ---------- Forwarded message ---------- Date: Thu, 08 Feb 2001 17:10:23 +0000 From: Yudi Pawitan <yudi@stat.ucc.ie> To: Mark Myatt <mark@myatt.demon.co.uk> Cc: R-Help <r-help@stat.math.ethz.ch> Subject: Re: [R] Goodness of fit to Poisson / NegBinomial

I'm putting a title over different plots, where 'that' varies. How do I make the value appear? Example: that <- 1 plot(1:10) title(expression(paste(hat(theta),'= ',that))) Yudi Pawitan yudi at stat.ucc.ie Department of Statistics UCC Cork, Ireland Ph 353-21-490 2906 Fax 353-21-427 1040 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

>I use WinNT. You have to "launch" R from a DOS-shell window; you 1st >change to the directory where you have RGui.exe, then type : > Rgui --vsize 15M --nsize 1000k > >and R opens with increased memore size. It works here with my NT box, and >will probably do with Win95 (though I think that the DOS-shells are not >strictly similar in both OSs but this may not be a

Dear R-helpers: (this is part of a bigger program) the following fails as a function, but runs OK if we comment out the fnfn_ function() line. Any hint would be appreciated. -Yudi- R : Copyright 1998, The R Development Core Team Version 0.63.0 Beta (Nov 13, 1998) -- on WIndows3.11 fnfn _ function (m=10,n=10,spar=2) { fn _ function(u,v){ uc_ u-floor(m/2)-1 vc_ v-floor(n/2)-1

(Sorry, I had the previous emails on a different machine.) Thanks to Prof Ripley for diagnosing the problem with inconsistent capitals under Win98. It's 'interesting': * I had ...\acepack\Src\Ace.f Avas.f and Rlsmo.f (tar'ed in a Unix machine and Ws_FTP'ed to the PC) * when the directory name was changed to ...\acepack\src it still didn't work *

Dear All, How can I get the same stream of random numbers in R and S? I know one of the generators in R is Super-Duper also used in S-plus, but the .Random.seed's are set differently. I tried to re-exrpess that in S-plus into two seeds according to the formula in VR, but I can't get the same random stream in R. Thanks, ------------------------------ Yudi Pawitan: yudi at ucd.ie

Dear All, How can I get the same stream of random numbers in R and S? I know one of the generators in R is Super-Duper also used in S-plus, but the .Random.seed's are set differently. I tried to re-exrpess that in S-plus into two seeds according to the formula in VR, but I can't get the same random stream in R. Thanks, ------------------------------ Yudi Pawitan: yudi at ucd.ie

Sorry, but I lost this thread, so I sending this as a new message. This is really a follow-up to a post from a couple days ago saying that fisher.test from the ctest library crashed on the following data set: > T [,1] [,2] [1,] 2 1 [2,] 2 1 [3,] 4 0 [4,] 8 0 [5,] 6 0 [6,] 0 0 [7,] 1 0 [8,] 1 1 [9,] 7 1 [10,] 8 2 [11,]

At 11:57 25/10/00 +0100, Brian Ripley wrote: >> Date: mer., 25 oct. 2000 12:38:55 +0200 >> From: Emmanuel Paradis <paradis@isem.univ-montp2.fr> > >> I think it would be nice to have par(ask=T) set by default in termplot(), >> like it is in plot.lm(). > >Well, it isn't really the default in plot.lm, the default for `ask' being > >interactive()

Is there facility in R to read binary file? In Splus scan() has 'width=' argument and AsciiToInt() command to translate the binary data. Is there something similar in R? Thanks, -Yudi- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

Hi, I think it would be nice to have par(ask=T) set by default in termplot(), like it is in plot.lm(). Best, Emmanuel Paradis -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !)

I had an error message while running a macro from Yudi Pawitan's web site: > source("ex2-13.r") Error in parse(file, n, text, prompt) : syntax error on line 2 Inspecting ex2-13.r I found that the error was generated by another source() command. Clearly R does not like nested source()s, which is fair enough when you think about it. Still it's something that you might want

Ape is an R package for "analyses of phylogenetics and evolution". The first version (0.1) has been released on 27 August 2002 and is available on CRAN. >From the 'Description' file of version 0.1: Ape provides functions for reading, and plotting phylogenetic trees in parenthetic format (standard Newick format), analyses of comparative data in a

There seems to be a bug with var() when the argument is a vector with exactly three values of 1e30 (or close to this value). This does not happen with twice, four (or more) times this value, or another value. > var(rep(1e30, 3)) [1] 2.971056e+28 > var(rep(1.2e30, 3)) [1] 2.971056e+28 > var(rep(0.9e30, 3)) [1] 2.971056e+28 > var(rep(0.8e30, 3)) [1] 0 > var(rep(1e29, 3)) [1] 0 >

Dear All, I have two distributions which I don't their nature. I want to check whether they come from the same distribution. I know that I can use KS test however the standart function ks.test applies only the ks test for testing the difference of two samples (non-parametric). By the way the distribution are of Euclidean distances. One of observed and the other of shuffled data. Thanks, Ron

Hi, Can anybody tell me how to open new a new window for plot()? Thanks. Minghua

Hi, If I have a CSV file which has several comments at the top, and the data start immediately after the line: @DATA Is it possible to use the scan() command to get the CSV data into R, by only reading the lines after @DATA? If so, how can I do it? Cheers, Kevin ------------------------------------------------------------------------------ /* Time is the greatest teacher, unfortunately it

On Feb 2, 2009, at 3:00 PM, Jon Harrop wrote: > On Monday 02 February 2009 20:37:47 you wrote: >> On Feb 2, 2009, at 12:39 PM, Jon Harrop wrote: >>> On Monday 02 February 2009 06:10:26 Chris Lattner wrote: >>>> I'm seeing exactly identical .s files with -msse2 and -msse3 on the >>>> scimark version I have. Can you please send the output of:

C?dric Finet wrote: > > I thank you for your answer but I do not understand yet why the Fisher?s exact > test does not work. And why is a "negative key". > > C?dric Finet > Running the original TOMS643 fortran code (R uses an f2c translation of this) says: FEXACT ERROR: 30 Stack length exceeded in f3xact. This problem should not occur. The integer hash key is

Dear R-helpers, I'm learning how to compile R packages with fortran routines and use the acepack source as an example. After some hiccups (Guido told me about the need to set make_mode=UNIX) I got to the following: `make pkg-acepack' runs without error (see below), but it does not * compile any .f file in acepack/src * (hence) try to create any dll file or libs (I'm