similar to: SVD of complex matrices

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

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

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

Does anyone know of a R package that will determine the convex hull of a high-dimensional dataset (say 4-10 dimensions). I know chull works for 2D data. I'm neophyte to R and convex hulls so please keep it simple. Many thanks Ben -- Ben Stapley. Biomolecular Sciences, UMIST, PO Box 88, Manchester M60 1QD. Tel 0161 200 5818 Fax 0161 236 0409

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

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

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

Hello All, I'd like to create diagrams for dynamic web-pages via a cgi-script that itself uses R. Sadly, I found out that the png and jpeg devices need an X-server :-( , so I cannot use these nice devices. Does any one have any experience using R for on-the-flycreation of graphics for web-pages pages? Any suggestions for workarounds? thanks very much for help, Arne -- Arne Mueller

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

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

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

Dear All, I'm getting confused with the concept R uses to do regression using lm. I'm afmiliar with gnuplot and the build-in fit command, but couldn't get R's lm to work on my data. I know that my data follows a powerlaw or maybe an exponential function, and I'd like to determine the best fitting factors for the formula: a*x^b where b < 0. I've tried thge follwoing:

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

Hello, I've just started R, and I'm getting a bit mad using it. I've managed to produce a barplot with the labels for the ytick marks placed horizontally (perpendicular to the y-axis) usiing par(las=1). The problem is that most of my labels are in part beyond the plotting area because they are rather long (e.g. "H. sapiens", "D. melanogaster" ...). What is the

>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

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

Hello, I've read in a data file as a frame and now I'd like the columns to be rows and the rows to be columns. What's the easiest way to do this in R? > class(d) [1] "data.frame" > rownames(d) [1] "98x101" "40x98" "30x40" "0x30" > colnames(d) [1] "H..sapiens" "C..elegans"