similar to: Mathematica and R

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 =

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

To help me make choices regarding a platform for running high-memory R processes in parallel, I would appreciate any responses to these questions: 1. Does the amount of RAM available to an R session depend on the processor (Intel vs. Sun) or on the OS (various Linux distributions vs. Solaris)? 2. Does R have any built-in limitations of RAM available to a session? For example, could it make use

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

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,

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

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:

Is there any R function that computes the convolution of the double exponential distribution? If not, is there a good way to integrate ((q+x)^n)*exp(-2x) over x from 0 to Inf for any value of q and for any positive integer n? I need to perform the integration within a function with q and n as arguments. The function integrate() is giving me this message: "evaluation of function gave a

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

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

I've run into a problem with the datetime axis generated by axis.POSIXct. It appears a similar issue was discussed in October 2003 under the subject line "datetime data and plotting" (see https://stat.ethz.ch/pipermail/r-help/2003-October/039071.html), but I wasn't able to determine whether there is a straightforward solution. The code below produces a graph with apparently

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

Is there a function that can shade the region between two arbitrary curves? For example, how can I fill in the area between these two plotted curves? > x <- 1:10 > y <- x + rnorm(10) > z <- x + 10 + rnorm(10) > plot(x, y, type = "l", ylim = c(0,20)) > lines(x, z) I would appreciate any help. David

I have some programs which were writen in mathematica or c language, but I donot know how to use these software. So I want to run them in R. Can I do it ? How to run "mathematica" or "c" programs in R? Jian Zhang [[alternative HTML version deleted]]

David - information on calling R from within Mathematica can be found at the following link: http://www.mofeel.net/1164-comp-soft-sys-math-mathematica/13022.aspx HTH, Alan Kelly

Hi all, I am an R newbie trying to do some calculations I do in Mathematica in R on a GNU/Linux system. The main thing I am interested in doing is taking a 0.999 quantile on a data set in a file who's data is normally distributed, say foo.csv. e.g in Mathematica if I have something like this : a=Import["foo.csv"] b=Transpose[a][[1]] Quantile[b,0.999] In R I can load all the data

Dear R-list, I try to transform a mathematica script to R. #######relevant part of the Mathematica script (* p_sv *) dd = NN (DsD - DD^2); lownum = NN (L-DD)^2; upnum = NN (H-DD)^2; low = lownum/(2s^2); up = upnum/(2s^2); psv = NIntegrate[1/(s^NN) Exp[-dd/(2s^2)] (Gamma[1/2,0,up] + Gamma[1/2,0,low]),{s,sL,sH}, MinRecursion->3]; PSV = psv/Sqrt[2NN]; Print["------------- Results

Hi, This question seems to be old, however I am asking it from a different point of view. R introduces itself as a 'statistical' language and environment. There are many discussions about comparing R to MATLAB or mathematica (or other similar software). It seems to me that these other software have a broader range of applications. For example, in Mathematica one can solve a partial

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