similar to: Convergent series

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

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

This patch series silences some unknown cap warnings and fixes up coding style (patch 1+4). The most important part of this series are the two patches in the middle. They rework the state_fb code, so that we are able to render to not 64 byte aligned targets, as this is the only real use-case for render temporaries this allows us to drop temp code completely and simplifies a lot

Dear Asterisk-Users: Question: ======== How do I get asterisk to pass DNID/RDNIS information between asterisk machines using iax2, in a Dial(IAX2...) command ? Setup: ===== I have two asterisk boxes, MASTER and SLAVE. MASTER is running 1.2.0 and SLAVE is running 1.2.1. The main box handles incoming calls on a multiple lines (both via hardware connection to our internal PBX and calls

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

Hi, I have a problem with R, after an update: this piece of code: cat("creating resolver data frame\n"); dfG<-cbind(dfG,2^(RGN$G)) dfR<-cbind(dfR,2^(RGN$R)) suddenly, creates values to inf. Is the syntax changed for the ^ symbol in the latest R downloadable version? Kind regards, wilfred -- Dr. Ir. Wilfred Van IJcken Labmanager Genomics Erasmus Center for Biomics,

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

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

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:

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

Hi there, I'm one of the developers of Hero Lab ( http://www.wolflair.com/hero_lab ), and I'm trying to investigate a problem that our users are having running Hero Lab under WINE. The user says he receives the following messages before Hero Lab hangs or crashes: > ~/.wine/drive_c/HeroLab$ wine HeroLab.exe > fixme:wininet:InternetSetOptionW Option INTERNET_OPTION_CONNECT_TIMEOUT

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

Thanks. Yes. Your approach can identify: Glaxy ace S 5830 and S 5830 Glaxy ace But you can not identify using same program: Iphone 4S 16 G Iphone 4S 16G How should I solve both in same time. Kind regards,Tammy [[alternative HTML version deleted]]

Hi R users: I'm trying to make a function where two of the parameters are functions, but I don't know how to put each set of parameters for each function. What am I missing? I try this code: f2<-function(n=2,nsim=100,fun1=rnorm,par1=list(),fun2=rnorm,par2=list()){ force(fun1) force(fun2) force(n) p1<-unlist(par1) p2<-unlist(par2) force(p1) force(p2)