similar to: [OT] "normal" (as in "Guassian")

Hi, Was wondering if it is possible to pass function name as a parameter, smth along this line param.to.pass<-c(1,'max','h') dd<-function(dfd, param=param.to.pass,...){ ttime.int <- format(ttime,fmt) data.frame( param[3] = tapply(dfd[,param[1]],ttime.int,param[3]), ...) } I know there is a as.formula expression but not quite sure if there is some

In another thread on this list, various wild allegations have been made, relating to the New York Times article on R. I object both to the subject line and to the content of several of the messages, and will not repeat or quote any of that content. It smacks to me of mischief making. Discussion has centered around the following quote from the NY Times article: ?According to them, the

Hello: would you please help me with the following glm question? for the R function glm, what I understand is: once you specify the "family", then the link function is fixed. My question is: is it possible I use, for example, "log" link function, but the estimation approach for the guassian family? Thanks, Shuangge Ma, Ph.D. ******************************************** *

Hi We are about to publish a book, which contains figures made with R plots. An important detail that we did not take into account is that the book will not be printed in 4 colors (cmyk mode), but only 2 (black +"spotcolor"). The spotcolor we use is part of the big Pantone family. The problem is that both pdf() and postscript() offer either rgb or cmyk, but no spotcolors such as

Hello I looked on the archives but did not find answer for that... Running R CMD check for a package, i get an error: Error in parse(n = -1, file = file) : unexpected symbol at 3341: } But how can I find which file is guilty? What is this 3342 referring to? Finally the solution I found is to source() every file until I find the file where the problem is... Is there some better way to do it?

Hi My goal is to do a (multiple) regression, just knowing that my Y variables will be the say k first variables of a matrix/data frame. I thought I should do it with eval(parse)) but encounter a strange problem. See: lm(y~.-y, data=freeny) #that's what I want to do in the one equation case #Problem is I don't know name of the variable... only that it is the first one... #so idea is to

hi I am seeing a strange behavior I can't understand... doing: > source("/tmp/RFile.r",echo=TRUE) Error in rep.int(c(prompt.echo, continue.echo), c(leading, length(dep) - : invalid 'times' value > traceback() 3: rep.int(c(prompt.echo, continue.echo), c(leading, length(dep) - leading)) 2: paste(rep.int(c(prompt.echo, continue.echo), c(leading, length(dep)

Hi I'm trying to install the rattle GUI on winwos 2000, with last version R 2.10. I had a first problem using the rattle package, as it was asking pkg XMl, which is no more in the CRAN repo. I instead simply install in R/libraries the XML pkg from: http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.10/ Maybe I mised something here... So now I was able to run the library(rattle) then,

Hi I have within a quite big function foo1, an internal function foo2. Now, in order to have a cleaner code, I wish to have the internal foo2 as "external". This foo2 was using arguments within the foo1 environment that were not declared as inputs of foo2, which works as long as foo2 is within foo1, but not anymore if foo2 is external, as is the case now. Now, I could add all those

Hi Say I want to add manually an intercept in the function lm. Even if almost all results will be identical, few stats are different as DF counting will be different as intercept will not be included in "automatic" case, while it will be in "manual" case. See: ###usual lm on freeny fr<-lm(freeny.y~freeny.x) ###manual lm on freeny man<-cbind(1,freeny.x)

Hello To update from R 2.6 to 2.8 (on Ubuntu 8.04 both) I had to install new tcl and liblapack packages (excuse me it is in french): > sudo apt-get install r-base-dev > Lecture des listes de paquets... Fait > Construction de l'arbre des d?pendances > Lecture des informations d'?tat... Fait > Les paquets suppl?mentaires suivants seront install?s : >

Hello I'm sorry if this suggestion/correction was already made but after a search in devel list I did not find any mention of it. I would just suggest to add a warning or an exemple for the help-file of the function unique() like "Note that unique() compares only identical values. Values which, are printed equally but in facts are not identical will be treated as different."

Please start a new thread for a new topic and use a meaningful subject for sake of the archives and everyone trying to follow. Using your Int1, correcting your Int2 and using your Y (which goes to Jul not Aug as per the comment): library(zoo) > as.ts(with(merge.zoo(Int1, Int2, Y, fill = 0), pmax(Int1, Int2))) Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 1992

Dear all, I'm trying to plot, in the same window, two different series, using barplot() for the first one and plot() for the second. What happens is that the second chart has a different axes origin, therefore the final plot is wrong. This piece of code shows the differences between the values of par()$usr: barplot(coredata(Z0), beside=TRUE, ylim=c(0,100)); par()$usr;

Thought I should copy the list with Matthieu's response. H -----Original Message----- From: Matthieu Stigler [mailto:matthieu.stigler at gmail.com] Sent: Wednesday, November 05, 2008 8:29 PM To: Horace Tso; rkward-devel at lists.sourceforge.net Subject: Re: [rkward-devel] questions on RKWard some answer only for the third question: do you have the package r-doc-html installed? (on ubuntu

Hi I need for a specific application to change the locale of R 2.9.2 in Ubuntu 9.04. Trying the example in ?Sys.setlocale: Sys.setlocale("LC_TIME", "de_DE.utf8") [1] "" Warning message: In Sys.setlocale("LC_TIME", "de_DE.utf8") : la requ?te OS pour sp?cifier la localisation ? "de_DE.utf8" n'a pas pu ?tre honor?e I tried the code

Hi I have a problem reading files from Windows... these files have, instead of NA on last column. special ending '<a0>' which makes problem... This problem does not appear while reading the same file in Windows! Try: read.csv("http://dl.dropbox.com/u/6113358/prob.csv") Could you please tell me if you also have this problem? I have tried either by cleaning the file on

Hello I think there is a conflict between setwd() and \SweaveOpts{prefix.string=}. In the same document, those both command get Sweave confuse the files and directories. See: say my .Rnw document is in File1 If one inserts some setwd() for another file: -setwd(File2) then the command \SweaveOpts{prefix.string=graphics/Rplots} will search the "graphics" folder in File2 because of

Hi I went through some extensive search to find suitable method (package, function) to fit multiple peaks. The best I found is ALS package but it requires rather complicated input structure probably resulting from GC-MS experimental data and seems to be an overkill to my problem. I have basically simple two column data frame with columns time and sig dput(temp) structure(list(time = c(33, 34,

Greetings! I am looking for advice regarding the best way to compare GAMMs. I know other model outputs return enough information for R's AIC, ANOVA, etc. commands to function, but this is not the case with GAMM unless one specifies the gam or lme portion. I know these parts of the gamm contain items that will facilitate comparisons between gamms. Is it correct to simply use these values