similar to: Help in the legend()

Hi I was trying to have a graph whose axes are X axis: m, Y axis: var[X ((a,b) in suffix, and (n,d) in the power)]. X ((a,b) in suffix, and (n,d) in the power)- X^(n,d) _ (a,b). Actually I require many plots involving different values of a,b,n,d, so need to keep this complicated notation. The expression() didn't work out for this case. Can anyone help me out. Thanks, in advance.

Hi, I want to use a more accurate value of exp(1).  The value given by R is 2.718282. I want a value which has more than 15 decimal places. Can anyone let me know how can I increase the accuracy level. Actually there are some large multipliers of exp(1) in my whole expression, and I want a more accurate result at the last step of my program, and for that I need to use highly accurate value

Hi, I wanted to compute the value of the function ifn at certain values of n. But I am receiving the following error when I was using the following code(given at the end). Error: evaluation nested too deeply: infinite recursion / options(expressions=)? I feel that since the function Grx is recursively related, perhaps making the code too complicated too handle. Can anyone let me know if

Hi, I wanted to compute the value of the function ifn at certain values of n. But I am receiving the following error when I was using the following code(given at the end). Error: evaluation nested too deeply: infinite recursion / options(expressions=)? I feel that since the function Grx is recursively related, perhaps making the code too complicated too handle. Can anyone let me know if

Hello, I want to draw a histogram of the mean of sample observations drawn from multivariate t distribution. I am getting the following error corresponding to the code I used. Though I am getting the graph, but I am curious to know the warning message. Warning messages: 1: In if (freq) x$counts else { : the condition has length > 1 and only the first element will be used 2: In if (!freq)

Hi, I was plotting or creating a map for Texas school districts using the shape file of Texas. I could not find any other helpful mail in the mailing list. txshp<-read.shape(system.file("S:\\Districts_10_11.shp", package="maptools")) Error- read.shape no found. But read.shape is there in maptools. If anyone can help me out it will be great. Thanks in advance. Shant

Hi, I was trying to have a graph whose axes are of the following type: X axis: n and Y axis: var[U ((a,b) in suffix, and (n,d) in the power)]. U ((a,b) in suffix, and (n,d) in the power)- U^(n,d) _ (a,b). Actually I require many plots involving different values of a,b,n,d, so need to keep this complicated notation. The code I used: plot(n, hn$h_pg,

Hello everyone! I want to use density function of R to compute the density at x(=0, say). But it is giving me the 5-number summary and mean of the data and densities at that point. I just want the densities at different values specified by me. Can anyone let me know how can I find that? For example Thanks in advance for your help. Shant [[alternative HTML version deleted]]

Hello R users, I have two data sets like the following. Form of dataset: data: X1 X2 X3 X4 X5 1902 RE 3 594 9 1903 RE 3 1340 7 1904 AA 3 760 14 1908 RE 4 1759 18 1909 EX 2 387 1 2901 AU 6 3116

Hello, I want to solve: x*(3^x)*log(4)-x*log(4/3)-(3^x)+1=0 for x. I used the following code, uniroot(function(x) x*(3^x)*log(4)-x*log(4/3)-(3^x)+1, lower = -2, upper = 2, tol = 0.001 ) While using this I am getting the following error. Can anyone please help me out. Error in uniroot(function(x) x * (3^x) * log(4) - x * log(4/3) - (3^x) + : f() values at end points not of opposite sign.

Hi! Can anyone please let me know what numerical integration procedure does the integrate function in package stats follow? Thank you in advance Shant [[alternative HTML version deleted]]

Hi, I want to find a value of n1. I used the following code but I am getting the error - Error in as.vector(x, mode) : cannot coerce type 'closure' to vector of type 'any' n=10 a_g<-(1/(n*(n-1)))*((pi/3)*(n+1)+(2*sqrt(3)*(n-2))-4*n+6) a_s<-function(n1) { t1=(n1-1)/2; (t1*(gamma(t1)/gamma(n1/2))^2)/2-1-a_g } xm<-solve(a_s) Can anyone help me out. Thanks in

Hello all, I was going to solve (n-m)! * (n-k)! = 0.5 *n! * (n-m-k)! for m when values of n and k are provided n1<-c(10,13,18,30,60,100,500)         # values of n kx<-seq(1,7,1);                               # values of k slv2<-lapply(n1,function(n){    slv1<-lapply(kx,function(k){              lhs<-function(m)              {                

Dear R users, I am new users of this software. I want to make box plot. Here, i have simulated data set with following commands: x<-matrix(rnorm(90),nrow=10,ncol=9) x a<-matrix(c(1,1,1,1,1,2,2,2,2,2),nrow=10,ncol=1) xx<-cbind(a,x)

n=10 x1<-(1:n)/n y1<-rnorm(n,x1^2,1) m=20 x2<-(1:m)/m The value of y2 will be rnorm (m, x2^2,1) if none of the elements of x2 is same as x1, but for every same elements in x1 and x2, the value of y2 will be same as y1. I know the following is correct, but for large vectors, this won't work efficiently as this is taking longer time. Can you please let me know how to write the

Hello: I have some trouble making a prediction from an AR(p) model. After I have the AR(p) model fitted , I want to use a new data set to make predictions. But I get the error: Error in newdata - object$x.mean : non-numeric argument to binary operator. A small version of my original data looks like: X1 X2 X3 X4 40813.65 1 10 41.86755 40813.65 1 8 41.86755 40813.66 1 8 41.86755

Dear list, I have some difficulty in manipulating list elements. More specifically, I am performing svm regression and have a list of lists, called pred.svm. The elements of the second list are 3D arrays. Thus I have pred.svm[[i]][[j]], with 1<=i<=5 and 1<=j<=20. I want to take the sum of the elements a specific array dimension across all j, for one i. Mathematically speaking, I want

Dear all, I have a data file with 3 variables (x1, x2, x3) where variable x1 have data that consists of several numbers separated by commas. id name x1 x2 x3 aa101 1,4,5 2 1 aa102 1,2,5 1 2 aa103 1,2,5 1 1 aa104 1,2,3 1 2 aa105 1,5 2 2 aa106 1,2,5 2 2 aa107 1,2,5 2 1 aa108 1,4,5 2 1 aa109 1,2 1 2 aa110 3,5 1 2 I want to count the number of data for each

please help me in this... I am unable to install pg gem file [root@vdimc04 ~]# gem install pg -v ''0.12.0'' /usr/local/lib/ruby/1.9.1/yaml.rb:56:in `<top (required)>'': It seems your ruby installation is missing psych (for YAML output). To eliminate this warning, please install libyaml and reinstall your ruby. Fetching: rake-compiler-0.8.1.gem (100%) Fetching:

OK, I''ve been putting off this question for a while now, but it eating at me, so I can''t hold off any more. I have a nice 8 gig memory stick I''ve formated with the ZFS file system. Works great on all my Solaris PC''s, but refuses to work on my Sparc processor. So I''ve formated it on my Sparc machine (Blade 2500), works great there now, but not on