similar to: Equivalent of TensorRank

Hi. Can anyone suggest a simple way to obtain in R a list of vector coordinates of the following form? The code below is Mathematica. In[5]:= Flatten[Table[{i,j,k},{i,3},{j,4},{k,5}], 2] Out[5]= {{1,1,1},{1,1,2},{1,1,3},{1,1,4},{1,1,5},{1,2,1},{1,2,2},{1,2,3},{1 ,2,4},{1,2, 5},{1,3,1},{1,3,2},{1,3,3},{1,3,4},{1,3,5},{1,4,1},{1,4,2},{1,4,3}, {1,4,

Hi, I've got an ugly but fairly simple function: mdevstdev <- function(a){ l <- dnorm(a)/(1-pnorm(a)) integrand <- function(z)(abs(z-l)*dnorm(z)) inted <- integrate(integrand, a, Inf) inted[[1]]/((1- pnorm(a))*sqrt((1 + a*l - l^2))) } I wanted to quickly produce a graph of this over the range [-3,3] so I used: plotit <-function(x=seq(-3,3,0.01),...){

Hi People, I have just started building my own RPMS on both Fedora and CentOS and generally things have gone well. Currently I am trying to create RPMS for some commercial software that we have purchased. Step 1 was install the software using its JAVA Based installer ensuring that all files were installed into a particular directory in /usr/local. Step 2 was create a tar file of the

Hi, this is not an R-help post, but I found this extract below that was written by a leading mathematician back in 1999 when he was talking about statistics and computing. I found it interesting to share and I ask your opinion do you think this still holds today or things have changed? Thanks. ?...we would also like to mention that in our opinion Mathematica provides an excellent and indeed

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

Dear R-list members, I have a problem with translating a mathematica script into R. The whole script is at the end of the email (with initial values for easy reproduction) and can be pasted directly into R. The problematic part (which is included below of course) is <--- Original Mathematica ---> (* p_svbar *) UiA = Ni (Dsi - 2Di A + A^2)/2; UiiA = Nii (Dsii - 2Dii A + A^2)/2; psvbar =

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

Full_Name: Skylab Gupta Version: R version 2.5.1 (2007-06-27) OS: Windows XP Submission from: (NULL) (216.82.144.137) Hello, I have been playing around with the statistical distributions in R. I think the computations for the cumulative distribution function of the students t distribution in R are not very accurate. For instance, the cdf of a students t distribution with 13 degrees of freedom

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

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

I just tried rendering a plot of 166 data point (time series) on my PPC mac, and compared it with the same data rendered by Mathematica. The R plotting is very slow to paint and repaint, compared to Mathematica. This also makes it extremely awkward to use in other documentation such as pdfLatex. On closer inspection, it looks like the Mathematica plot is just a polyline, while in R it is a

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

Hi, I need to solve an optimization problem in R having linear objective function and quadratic constraints(number of variables is around 80). What are the possible choices to do this in R. optim() function only allows box constrained problems. Is it possible in nlm()? Or please tell me if there is any other routine. Thanks Amit

I know that at the end of an R session I'm given the option to save the current *state* of the session. But I would like to save the entire sequence of inputs that took place during the session, so that I can play them back later, and not only be left in the same state I was at the time of saving the session, but be able to see the entire history of the session (inputs and outputs). (This is

I think I should have posted this question here as well. I am posting my question here since it is R related. Please see below. I originally posted this to sci.stat.math "Nasser Abbasi" <nma at 12000.org> wrote in message news:Mm4gj.28843$R92.4987 at newsfe16.phx... > > I think R documentation is a bit hard for me to sort out at this time. > > I was wondering if

I'm making some illustrations and it would be convenient to automatically shade the overlapping portions of circles. These illustrations are for Social Choice theory, a field in political science and economics. I've wrestled together some examples so you can see what I mean, but have not mastered the "color overlapping sections" problem (as you will see):

Using Mathematica, I've performed a hierarchical clustering of 3,107 U. S. counties based on 1995-2000 intercounty migration, using an algorithm of my own devising (see http://arxiv.org/abs/0807.1550). I can also generate the associated 3,107 x 3,107 matrix of ultrametric distances. Unfortunately, as far as I have been able to determine, the Mathematica hierarchical clustering package only

Hi, I've used the Mathematica to produce 3D graphics, contour plots of a bivariate normal distribution Now I want make these graphics in R, but i do not know how. I would like to: - Plot a 3D graph for some different variance matrix - Plot the contour plots - Find and try to plot (in the 3d graph ou contour plot) the (1-a)% confidence region based in a chi-square(a) with the degrees of

Hello I am trying to evaluate an Incomplete gamma function in R. Library Zipfr gives the Igamma function. From Mathematica, I have: "Gamma[a, z] is the incomplete gamma function." In[16]: Gamma[9,11.1] Out[16]: 9000.5 Trying the same in R, I get > Igamma(9,11.1) [1] 31319.5 OR > Igamma(11.1,9) [1] 1300998 I know I have to understand the theory and the math behind it rather