similar to: Extending the accuracy of exp(1) in R

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

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

Hi R-users, I was plotting the differences of the variances of the three estimators- T^(1), T^(2), T^(3), ofcourse taking two at a time. I was using the expression() in the legend function in order to show which line correspond to which of the difference, but the following that I had used didn't gave desired result. I would be grateful, if you help me out. plot(n, 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]]

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

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

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)

Do you know if either of the patches referenced above can be applied to the 3.6p1 code shipped by Red Hat (RHEL3) and Suse (SLES8)? It appears they were integrated in the 3.7p1 and 3.9p1 releases, respectively. Steven Vandor Computer Associates 425 825 2799 vandor at ca.com This message may contain confidential, privileged or other legally protected information.? It is intended for the addressee

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

Hi Preston, >>> As noted in the GEP FAQ, GEPs don't support variable-length arrays; >> >> that's not quite right. The problem is only with arrays of variable length >> arrays, and more generally with arrays where the element type has variable >> size (this occurs with Ada, which has all kinds of funky variable sized types, >> for example). > >

Duncan Sands wrote: >> As noted in the GEP FAQ, GEPs don't support variable-length arrays; > > that's not quite right. The problem is only with arrays of variable length > arrays, and more generally with arrays where the element type has variable > size (this occurs with Ada, which has all kinds of funky variable sized types, > for example). You're right, though

I am trying to produce rather mundane output of the form e.g. pi, e = 3.14 2.718 The closest result I achieved so far with print() is: > print (c(pi, exp(1)), digits = 3) [1] 3.14 2.72 > print(c("pi, e =", pi, exp(1)), digits = 3) [1] "pi, e =" "3.14159265358979" "2.71828182845905" I understand that c() promotes floats to strings and