search for: pawitan

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 Sorry, Mark, the program worked for Rw1011, but this is what happened in Rw1021: > dnbinom(0:5,.9,.4) [1] 0.4383833 0.0...

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 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help&...

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

...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 Department of Statistics, UCD Dublin 4, Ireland Ph : 353-1-706 7641 Fax: 353-1-706 1186 ------------------------------ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send &qu...

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

...from Explorer), move it to somewhere sensible. Then, click the right mouse button on the shortcut icon, and go to Properties, then Shortcut section, that's where the command line is (look for Target), where you can add the vsize and nsize option. -Yudi- ------------------------------ Yudi Pawitan: yudi at stat.ucc.ie Department of Statistics, UCC Cork, Ireland Ph : 353-21-902 906 Fax: 353-21-271 040 ------------------------------ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send &...

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

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 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 Department of Statistics, UCD Dublin 4, Ireland Ph : 353-1-706 7641 Fax: 353-1-706 1186 ------------------------------ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send &qu...

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.

...: ld returned 1 exit status MAKE.EXE[2]: *** [acepack.dll] Error 1 MAKE.EXE[1]: *** [src/acepack.dll] Error 2 MAKE.EXE: *** [pkg-acepack] Error 2 C:\R\rw1000\src\gnuwin32>path PATH=C:\R\TOOLS;C:\H-I586~1\BIN;C:\GCC-29~1.2\BIN;C:\PERL\BIN;C:\WINDOWS;C:\ WIND OWS\COMMAND;C:\DMI95\WIN32\BIN Yudi Pawitan yudi at stat.ucc.ie Department of Statistics UCC Cork, Ireland Ph 353-21-490 2906 Fax 353-21-427 1040 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help&...

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

...2) } x1_ 1:m ; x2_ 1:n x12 _ outer(x1,x2,FUN='fn') filt _ x12/sum(x12) image(filt) } > source('fnfn.r') > m_10; n_10; spar_2 > fnfn(m,n,spar) Error in get(x, envir, mode, inherits) : variable "fn" was not found ------------------------------ Yudi Pawitan: yudi at ucd.ie Department of Statistics, UCD Dublin 4, Ireland Ph : 353-1-706 7641 Fax: 353-1-706 1186 ------------------------------ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send &qu...

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

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

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

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

...en <maj at stats.waikato.ac.nz> | Subject: [R] Nested source()s | To: R-help <r-help at stat.math.ethz.ch> | Message-ID: <40F2533C.9080709 at stats.waikato.ac.nz> | Content-Type: text/plain; charset=us-ascii; format=flowed | | 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 abo...