similar to: How to run "mathematica" or "c" programs in R?

I think there is a question in R. I donot know the reason. This is my data about comulative percentage figure. The result is not right. The first point (no=1,leiji=26.94350) in the plot figure was showen in a lower location. Why? Thanks ! > fre no leiji 1 1 26.94350 2 2 46.86401 3 3 60.59032 4 4 72.17355 5 5 77.85521 6 6 82.05853 7 7 85.56495 8 8 87.64378 9 9

I want to get the data in a circle from a large square plot. If the data includes only two columns or more columns which is number style, I can do it. But if the data includes other columns which is not number style (e.g. species name), I donot know how to do it. And I donot want to delete the columns because it is useful to me. e.g. species x y a 12

Hey, folk. I am trying to get many figures by using the function "xyplot" in the library "lattice". I tried to using the loop to finish it quickly. But I cannot open the saved file after I run the program. But if I use the function "plot" to get other simply figures, it can work. So I want to ask how I can get the saved figures by using "xyplot". The

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

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

Dear All, I just read an announcement saying that Mathematica is launching a version working with Nvidia GPUs. It is claimed that it'd make it ~10-100x faster! http://www.physorg.com/news146247669.html I was wondering if you are aware of any development going into this direction with R? Thanks for sharing your thoughts, Best wishes, Emmanuel

I can load the package in other R file, but I can not do it in the file. How to solve the question? > local({pkg <- select.list(sort(.packages(all.available = TRUE))) + if(nchar(pkg)) library(pkg, character.only=TRUE)}) Loading required package: tcltk Loading required package: car Error in parse(file, n, text, prompt) : syntax error in "*" Error: .onAttach failed in

Hi! I want to choose part of the data. But it has a question. > attach(PIKO) > PIKO[1:5,] code species dbh x y tag status branch 92 10142 PIKO 38.9 6.05 12.81 10165 alive 0 109 10213 PIKO 41.0 6.71 26.21 10202 alive 0 135 10222 PIKO 48.3 18.98 21.28 10214 alive 0 146 10223 PIKO 47.0 19.45 22.83 10224 alive 0 147 10223 PIKO 20.0 19.74

Hi,everyone. I use the following program calculates Fisher's alpha from counts of individuals and species. The program is wrote by Prof. Kyle Harm. However, when I run the program, it can work very quickly sometimes, but it can not work very well sometimes. It depends on the counts of individuals and species. For example, > calc.alpha(1000,70) [1] 17.14375 > calc.alpha(10000,70) [1]

I want to replace some objects in one row or column.For example, One colume: a,b,a,c,b,b,a,a,c. I want to replace "a" with "1", "b" with "2", and "c" with "3". Like this: 1,2,1,3,2,2,1,1,3. How to do it? I donot know how to do it. Thanks. [[alternative HTML version deleted]]

Hi all, Having seen a few messages regarding the ASUS Eee PC (p701) and R, and then buying one and trying to get R installed on it, I thought there might be a few people on the list who would appreciate a cookbook on how to do it. It's a nice little piece of machinery, running Xandros 4 Linux (seems to be compatible with Debian Etch) and set up for the person who wants to do all those

I want to load many files in the R. The names of the files are "Sim1.txt", " Sim2.txt", "Sim3.txt", "Sim4.txt", "Sim5.txt" and so on. Can I read them at one time? What should I do? I can give the same names in R. Thanks. For example: > tst=paste("Sim",1:20,".txt",sep="") # the file names > tst [1]

Hi! I want to use the package "SPLANCE" to do Ripley's K analysis I have a basic question about the package. I try to find how to do it, but I cannot. I hope you can help me. I donot know how to change my data form to "point data set". my data form is: > qumo[1:5,] code species dbh x y tag status branch 223 10312 QUMO 54.5 7.83 44.80 10306 alive

I want to gain thousands of random sampling data by randomizing the presence-absence data. Meantime, one important limition is that the row and column sums must be fixed. For example, the data "tst" is following: site1 site2 site3 site4 site5 site6 site7 site8 1 0 0 0 1 1 0 0 0 1 1 1 0 1 0 1 1 0 0 0 1 0 1 0 0 0 0 1 0 1 0 1 1 0 1 0 0 0 0 0 0 1 0 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 1 0 1 0

I want to transform the data by "tapply" to one dataframe. But I can not get it. For example: > tst=tapply(point,pp,length) > tst[1:10] p1 p10 p100 p1000 p1001 p1002 p1003 p1004 p1005 p1006 1 5 1 8 6 5 8 7 4 4 > res=as.data.frame(tst) # I try to transform it > res[1:10,] p1 p10 p100 p1000 p1001 p1002 p1003 p1004 p1005 p1006 1 5 1 8 6 5

Hi. I want to get a series data just like this: sp1 sp2 1 2 1 3 1 4 1 5 2 3 2 4 2 5 3 4 3 5 4 5 I can get one part of the data every time, but I can get the data directly. I try to use the loop, but it can not work. How to get it ? Thanks. pair.fn=function(i,sp=5){res=data.frame(sp1=rep(i,(sp-i)),sp2=(i+1):sp) return(res)} pair.fn(1) pair.fn(2) ..... [[alternative HTML

If I do not add "ress" into the data frame "res", there is no quote in the data frame. However, I add "ress", all column were found the quote. How to remove it? If you can delete the quote in the file "ress", that is better. Thanks. > ress[1:10] [1] "ABHO.ABNE" "ABHO.ACBA" "ABHO.ACGI" "ABHO.ACKO"