similar to: Summing up matrices in a list

Following the recommendation of Prof. Ripley, I have created the Makevars file with the line: PKG_LIBS="-L/usr/lib/libm -lm -L/usr/local/lib/libgsl -lgsl -L/usr/local/lib /libgslcblas -lgslcblas" in the working directory. Now I have the code file Example3.c which computes the Bessel function value (the example is taken from the GSL reference book). I am running: R CMD SHLIB Example3.c

Dear experts! I have a minimization problem with non-linear constraint and Objective function(theta)=lambda*(Constr)^2-f(x,theta). Theta is a vector of parameters. I'd like to stop the optimization after the value of the constraint is less or equal some constant value, say d, and save the last computed value of the function. For this purpose, I thought to define the Objective function like

Dear all, I'd like to find the mle esttimates using the mle function mle(negloglik, start = list(), fixed=list(), method="..."). I am using the L-BGFS-B method and I don't supply the gradient function. Is there a way to print the gradients found at the solution value? I am using R-1.9.1 on Windows and on Unix. Thank you in advance, Victoria Landsman. [[alternative

Is there a simple way to calculate summary statistics for all the matrices or dataframes in a list? For example: > z <- list(matrix(c(2,2,2,2), ncol = 2), matrix(c(4,4,4,4), ncol = 2)) > z [[1]] [,1] [,2] [1,] 2 2 [2,] 2 2 [[2]] [,1] [,2] [1,] 4 4 [2,] 4 4 > I would like to calculate, for example, the mean value for each cell. I can do that the hard

Dear experts! Is there a package that enables to create the bivariate log-normal variables? Thanks a lot, Vicky Landsman. [[alternative HTML version deleted]]

Dear all, I need an advice in the following problem. I have to maximize two functions of the form f1(x)=f(y1,x,alpha1,beta1) and f2(x)=(y2,x,alpha2,beta2), the maximization is with respect to alpha1, alpha2, beta1, beta2. I can maximize each function separately using nlm. The problem is that I have to add the constraint of the form g(alpha1)=g(alpha2). The total number of parameters is

Hi R, I have a list of matrices. I need to get the sum of all the matrices in the list. Example: a=b=c=d=matrix(1:4,2,2) l=list(a,b,c,d) I need: > a+b+c+d [,1] [,2] [1,] 4 12 [2,] 8 16 Something like do.call("+",l) is not working...why is this? I may not be knowing the number of matrices in the list... Thanks, Shubha This e-mail

Dear all, I'd like to find a minimum of (-loglik) function which is a function of k parameters. I'd like to run the minimization algorithm for the different subsets of the parameters and assign the fixed values to the complementary subset. How should I define my (-loglik) function such that it can be passed to the optim or other optimization function? Much thanks for any suggestions.

>From the `R Data Import/Export' manual, shipped with R, as well as available from the official R web site (last three paragraphs of Section 3.1, describing functions in the `foreign' package): Function read.S which can read binary objects produced by S-PLUS 3.x, 4.x or 2000 on (32-bit) Unix or Windows (and can read them on a di erent OS). This is able to read many but not all S

Dear all, I need to compute (numerically) the two-dimensional integral: int(int f(x,y)dy)dx. What is the more efficient(fast) way to do it? Is adapt function appropriate for this problem? I will much appreciate your help and attention. Vicky. [[alternative HTML version deleted]]

Dear R-list people, I asked a similar question a few hours before. I will try to be more specific. We like to add the GSL library to the file SHLIB in order to make it possible to run the C code using GSL functions from R. We read that the path to the libgsl.a should be added to the line shlib_libadd=' ' in the file SHLIB but it does not work on our system. Dyn.load fails with error

Can someone help me calculate correlations for grouped values? Heres what my first few line of data look like: > head(cmexpr) LLID GMID CEXPR MEXPR 1 1005 10831 2.057462 -0.08486 2 1005 10831 2.057515 -0.08486 3 1005 10831 2.057462 0.01209 4 1005 10831 2.057515 0.01209 5 1005 10836 2.050980 0.17237 6 1005 10836 2.018576 0.17237 LLID is gene id, GMID is cell line id, the EXPR

Deall all, We need to perform a statistical analysis of a large database (40,000 entries with approximately 500 fields in each entry) currently handled in Oracle. The data contains categorical variables only. At the current stage we suggest classification and clustering analysis. We are planning to perform the analysis in R and would be very grateful for any

Hi How do I make do.call() take "+" as a function for a list of more than two elements? Toy problem follows: f <- function(i){matrix((1:6)^i,2,3)} # Thus f() returns a matrix of size 2x3; I want to add a whole bunch of such matrices, # as in f(1) + f(2) + f(3) + f(4) # But: > do.call("+",sapply(1:4,f,simplify=FALSE)) Error in do.call("+",

I have found the use of save( ) and attach( ) when supported by a pair of functions written by my colleague John Miyamoto, move( ) and rm.sv( ) quite useful in managing (1) collections of useful homebrew functions, (2) project workspaces, and (3) "packages" under development. An .Rdata file containing these and other handy functions together with a brief supporting document can be

Couldn't you do this by subtracting 0.5 + x from your y values and checking for normality with mean 0 and sd = 1 (using ks.test or another test of normality). If you fail, you'll have to do additional work to find out whether pairs with some particular x value (or range of x values) is causing the problem, but I think this fits the question as stated. Of course, if you have discrete x

This is not a direct answer to your question, but if all you want to do is integrate a function numerically, you don't need to be doing it through R. You should be able to find codes that you can use directly, such as quadpack or even GSL. You may want to search on GAMS (Guide to Available Mathematical Software, http://gams.nist.gov/). Andy > From: Victoria Landsman > > Dear all,

Adjust drivers that support the 'FLOW_DISSECTOR_KEY_META' key to reject filters that try to match on the newly added layer 2 miss option. Add an extack message to clearly communicate the failure reason to user space. Example: # tc filter add dev swp1 egress pref 1 proto all flower skip_sw l2_miss true action drop Error: mlxsw_spectrum: Can't match on "l2_miss". We have

Dear all, I have tried a lot of things to make broadvoice work with asterisk , but I failed each time. Please suggest a good service providers that I can use with asterisk for outbound and inbound calls. -- With regards, Vicky Shrestha System Director WorldLink Communications Jawalakhel , Kathmandu, Nepal

Just did some testing.. root at vicky:~# cat /var/log/auth.log | grep "Set" May 30 21:41:05 vicky sshd[1568]: Set /proc/self/oom_adj from -17 to -17 May 30 21:41:07 vicky sshd[1574]: Set /proc/self/oom_adj to -17 root at vicky:~# ps faux | grep 1574 root 1574 0.0 0.0 70488 3404 ? Ss 21:41 0:00 \_ sshd: root at pts/1 root at vicky:~# ps faux | grep