similar to: handling recursion relation

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

Dear All! If I pass an object to an assignment function I cannot get it's name by deparse(substitute(argument)), but I get *tmp* and I found no way to get the original name, in the example below it should be "va1". Is there a way? Thanks, Heinz ## example 'fu1<-' <- function(var, value) { print(c(name.of.var=deparse(substitute(var))))} fu1(va1) <- 3 name.of.var

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.

Hello All: Can anybody tell me what is the problem with my program please. I have an error message as appears below. My program is: ifn <- "Jul08_09.LM" data <- read.table(ifn) ofn <- "Jul.png" bitmap(ofn, type = "png256", width = 30, height = 30, pointsize = 30, bg = "white",res=50) par(mar=c(5, 5, 3, 2),lwd=5)

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

channel <- odbcConnectExcel("nuova tabella terapia occupazionale mod.xls") > ## list the spreadsheets > sqlTables(channel) TABLE_CAT TABLE_SCHEM TABLE_NAME 1 c:\\TABELLE DEFINITIVE\\nuova tabella terapia occupazionale mod <NA> 'emi tot 2006 OAI_60g TO FU1$' 2 c:\\TABELLE DEFINITIVE\\nuova

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 all, I have this simple barplot code: ifn <- "id.dat" dat <- read.table(ifn) ofn <- "id.png" bitmap(ofn, type = "png256", width = 30, height = 30, pointsize = 30, bg = "white",res=50) par(mar=c(5, 5, 3, 2),lwd=5) par(cex.main=1.6,cex.lab=1.6,cex.axis=1.6) names(dat)<-c("NumberOfPeople","Average")

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 there, I have the following info in MyData: 01 0.40 02 0.40 03 0.40 04 0.35 05 0.34 06 0.33 I want to read them using the following code: ifn <- "MyData" dat <- read.table(ifn) names(dat1)<-c("Code","M") names(dat2)<-c("Code","M") I want to have the first three lines in

Hello All: I am having the following data file named "data.dat" Number of people Number of pets Age of trees ifn <- "data.dat" dat <- read.table(ifn) colnames(dat)<-c("Variables") I want R to read all these in a string but when I ask R to read these, it gives me error because there more than one word in each line. Any comments to solve this problem?

Hello All: Could anybody tell me please what is the meaning of this error message: ESP Ghostscript 815.02: Unrecoverable error, exit code 1 Error in plot.xy(xy, type, ...) : ignoring SIGPIPE signal Calls: plot -> plot.default -> plot.xy Execution halted My code has to go to 18 loops for three times and in each loop it has to read a file which is 100-300 Megs. It always crashes at loop 14

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

Dear Fellows, I would like to know how to uncompress a gz file at the R console. I could not find out any help from the R-help archive. Thanks for your great help. Best Regards, Wonsang You ----- -- Wonsang You Special Lab Non-Invasive Brain Imaging Leibniz Institute for Neurobiology http://www.ifn-magdeburg.de -- View this message in context:

I am unable to get a my Draytek working with our Asterisk server. I can make/recieve calls but get no audio. I have tried the various codecs at the Vigor end but still getting nothing. I looked at sip debug (below) but am new to Asterisk and don't really know what I am looking for. Asterisk works fine with XLITE so I know my installation is ok. Sip read: INVITE

Dear fellows, I am a novice in R. I would like to combine a matrix and a vector. Assume that we have the matrix a and the vector b with same length of column. a<-matrix(seq(1:10),nrow=2,ncol=5,byrow=TRUE) a= 1 2 3 4 5 6 7 8 9 10 b<-t(c(11,12,13,14,15)) b= 11 12 13 14 15 Then, I want to combine a and b as follows. c= 1 2 3 4 5 6 7 8 9

Hello, I am trying to write math-type on a plot. Due to space limitations on the plot, I want 2 short expressions written on top of each other. It is certainly possible to write them in two separate calls, but that involves fine-tuning locations and a lot of trial and error (and I'm trying to write a general purpose function). Here's where I've gotten to:

Hi, R-Helpers, I would like to ask about multiple graphs in one figure. I tried to execute the following codes. xlim <- c(1,100) ylim <- c(1,4) plot(NA, xlim=xlim, ylim=ylim) > x <- c(1:100) for(j in seq(1,10,by=1)) { y <- j*x^2+log(j) lines(x, y) } In the above codes, I had to arbitrarily set up the coordinate range of the figure in advance before

How might you fit a generalized linear model (glm) with variance = mu+theta*mu^2 (where mu = mean of the exponential family random variable and theta is a parameter to be estimated)? This appears in Table 2.7 of Fahrmeir and Tutz (2001) Multivariate Statisticial Modeling Based on Generalized Linear Models, 2nd ed. (Springer, p. 60), where they compare "log-linear model fits to

Hello, In the past I have used "expression" to include greek letters in axis labels, but this time I need to include the greek letter as part of a legend. Basically, I need to create the following vector to rename the levels of a factor: c("Interferon-gamma", "IL-10", "IL-5"), where "gamma" obviously needs to be printed as the greek letter