similar to: reading binary file

(RFC := Request For Comments) >>>>> "Tim" == Tim Keighley <Tim.Keighley@csiro.au> >>>>> on Thu, 23 Oct 2003 11:45:22 +1000 writes: Tim> Hi Martin, Tim> In October 2000 you wrote to r-help: >>> which reminds me that I've had a desire for something like >>> the old S function [from the blue book, and

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

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

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

Is there a way to determine the SVD of a complex matrix using R? (I'm using v1.0.1 and svd() won't do the trick). I know LAPACK has a function to do this. Thanks -- Ben Stapley Biomolecular Modelling Lab Imperial Cancer Research -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

I'm trying to get a grasp of maximum-likelihood estimation and would like to find a package that performs mle (hopefully a simple example). It seems as if there are plenty of packages that make use of different types of likelihood estimators, but none are of a simple, "newbie" type. Does anyone have a suggestion for which package would be the best for a mle example? Thanks, Jeff.

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

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

Appology in advance if this has been asked and answered. I am getting different answers using Fisher's exact test in the package ctest. For example: x_cbind(c(8, 22), c(31, 29)) fisher.test(x) Fisher's Exact Test for Count Data data: x p-value = 0.04024 alternative hypothesis: two.sided However, in SAS or other packages on the net the answer I get is p-value= 0.02664. Is this

>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 Sir, I am a new user of R and I am doing a tast, which is: find the maximum likelihood estimate of the parameter of Gaussian distribution for generated 100 numbers by using >x=rnorm(100, mean=3, sd=1). I tried to use following Maximum Likelihood function >fn<-function(x) (-50*log((sd(x))^2))-50*log(sqrt(2*pi))-(1/2*((mean(x))^2))*(sum((x-(mean(x))^2)), but it did not work. I am

(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 *

I'd like to put a degree symbol in the title of a plot. Something like: plot(1,main=expression("Data at " * 45 * degree * " North")) That works fine, but I like the 45 to be the value of a variable, instead of a constant, an in: lat <- 45 plot(1,main=expression("Data at " * lat * degree * " North")) But this doesn't work of course, the graph

All, I have some data on parasites on apple leaves and want to do a goodness of fit test to a Poisson distribution. This seems to do it: mites <- c(rep(0,70), rep(1,38), rep(2,17), rep(3,10), rep(4,9), rep(5,3), rep(6,2), rep(7,1)) tab <- table(mites) NSU <- length(mites) N <-

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

At 13:51 26/08/00 +0100, Yudi Pawitan wrote: > >If you have an actual large table to analyse, rather >than trying to solve the space problem, you may want >to consider a MonteCarlo implementation >of the exact test. Very easy to implement in R. See, >for example, Lange's Numerical Analysis for Statisticians, >Section 21.7. > >-Yudi- Yes, but this does not seem

My attempt to get R running in Windows 3.11 didn't go very far. I have read the FAQ and install.tex carefully. What am I missing? Here is the step: 1. From statlib I got rexe.zip (and rbeta.zip, but this one seems for windows 95 only?). 2. On unzipping I got 8 files, including rb1.exe, and a sub- dir win32s, which I copied to windows/system. (There is no setup.exe,

Hello R-Help, did somebody ever implement some hazard smoothing routine in R (that allows, say, some group comparison as for the Kaplan-Meier curves.) As far as I can tell it is not in survival5 package. Thanks, -Yudi- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",