similar to: convolution of the double exponential distribution

Given a numeric vector of observations, does R have any generic way to estimate the parameters of commonly used distributions (exponential, gamma, etc.) without numerically optimizing the likelihood function? Thanks, David _______________________________________ David R. Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International Bioinformatics & Exploratory Research 7250

Does anyone know how to use the equivalent of debug() on an S4 method? I would like R to enter the browser not for the generic function, but for the method of the class that I specify. Thanks, David _____________________________ David Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International Bioinformatics & Exploratory Research 7250 NW 62nd Ave., PO Box 552 Johnston,

Is there a function that provides the Fisher information matrix for a generalized linear model? I do not see how to access the off-diagonal matrix elements of the value returned by glm. (I'm particularly interested in logistic regression.) If not, what is a good way to use R to compute Hessians or other partial derivatives of log likelihoods? I would appreciate any guidance. David

Is there a way I can get R to automatically enter the browser inside a user-defined function on the generation of an error? Specifically, I'm trying to debug this: Error in as.double.default(sapply(lis, FUN)) : (list) object cannot be coerced to double In addition: There were 38 warnings (use warnings() to see them) > traceback() 8: as.double.default(sapply(lis, FUN)) 7:

While I think total compatibility is not possible or even very important I have become aware recently of some fairly low level incompatibilities that have rather serious consequences for interchange between systems. Some of these are possibly well known, but I put them on on record here just in case. (I do not expect any action on these, of course, but people may want to note them, and correct

While I think total compatibility is not possible or even very important I have become aware recently of some fairly low level incompatibilities that have rather serious consequences for interchange between systems. Some of these are possibly well known, but I put them on on record here just in case. (I do not expect any action on these, of course, but people may want to note them, and correct

After reading the help page on set.seed, I am unsure about how to solve the following problem. I need to call function f a thousand times. The list of values returned by f, should be as random as possible. f calls g twice: f <- function(){g1 <- g(1); g2 <- g(2); c(g1; g2)} The function g in turn calls sample and returns a number, but also depends on its argument, so, starting from the

Is there a way to call an R function to work in the background within the same R session? What I have in mind is the equivalent of adding '&' at the end of a UNIX command, so that R can process an intensive command while I execute other R commands. I want to do this in the same R session, rather than with multiple calls to R, so that all my modified objects will stay in the same

Is there any R function that can display a Postscript file that is already in the working directory? For example, if 'graph.ps' is such a file, I'd like to type something like this: > plot.postscript.file(file = 'graph.ps') If no such function exists, I'd be interested in a way to use existing R functions to do this under UNIX or Windows, preferably without a system

Are there any R functions or packages that can compute distributions, expectations, or quantiles of order statistics (or sample quantiles or extreme values) for a given distribution such as a normal distribution? Both exact and asymptotic calculations are of interest. I am already aware of the 'quantile' function of 'stats'. David _______________________________________ David R.

Is there a way I can have R automatically execute the commands in a source file without ever having to use R interactively? If so, what arguments should I pass to the UNIX call to R? I need to do this to run several R jobs in parallel. An alternative may be to have R and an R package behave as a stand-alone application that can be called from the UNIX command line. Is there any documentation on

R never gave me any serious problems until I tried to use it on Windows XP. I had 1.8.1 installed and working successfully for a few weeks. This morning when I launched it, I received this error message: "R for Windows GUI front-end has encountered a problem and needs to close. We are sorry for the inconvenience." I downloaded 1.9.0 and installed it, but when I launch it, I get this

What are some effective ways to leverage the strengths of R and Mathematica for the analysis of a single data set? More specifically, are there any functions that can assist with any of the following? 1. Calling an R function from Mathematica. 2. Calling a Mathematica function from R. 3. Using XML or another reliable data format to pass vectors, matrices, and/or lists from one environment to

Does anyone know how to interrupt R in RedHat? Neither control-c nor Esc is working. What I don't want to do is close the window or kill the entire R process. Thanks, David This communication is for use by the intended recipient and ...{{dropped}}

I am having trouble with the library function of R Version 1.5.0 for Mac OS X Carbon when I try to load libraries that did not come with the R installation. For example, I put the Biobase library from bioconductor.org in the 'library' folder. R cannot find it, even though it can find the pre-installed libraries: > library() # this lists standard libraries such as 'boot',

Is there an R function that automatically writes all input and output to a file? I would at least like it to log all the commands I enter, and to preferably also write the standard output to the file as well as to the screen. (The ideal would be to write the input to one file and both the input and output to another file.) I tried R2HTML for this, but I could not get it to work consistently. I am

Could anyone recommend a package for visualizing a likelihood function of two scalar parameters? I would prefer a three-dimensional plot similar to the kind Mathematica is known for, perhaps generated by a package not specific to likelihood. David ______________________________ David R. Bickel Ottawa Institute of Systems Biology BMI Dept., University of Ottawa http://www.oisb.ca/members.htm

Is there any automatic mechanism for extracting a likelihood or test statistic distribution (PDF or CDF) from an object of class "htest" or from another object of a general class encoding a hypothesis test result? I would like to have a function that takes "x", an object of class "htest", as its only argument and that returns the likelihood or test statistic

What are some reliable R functions that can compute the value of a convergent series? David -- David R. Bickel, PhD Associate Professor Ottawa Institute of Systems Biology Biochem., Micro. and I. Department Mathematics and Statistics Department University of Ottawa 451 Smyth Road Ottawa, Ontario K1H 8M5 http://www.statomics.com Office Tel: (613) 562-5800 ext. 8670 Office Fax: (613) 562-5185

Neither biocLite nor the GUI menus can install packages on my system. Here is relevant output: > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 11.1 year 2010 month 05 day 31 svn rev 52157 language R version.string R version 2.11.1 (2010-05-31) > source("http://bioconductor.org/biocLite.R") BioC_mirror =