similar to: Error in parse(outFile)

Dear all I am tryiing to fit the following lme with an ARMA correlation structure: test <- lme(fixed=fev1f~year, random=~1|id2, data=pheno2, correlation=corARMA(value=0.2, form=~year|id2), na.action=na.omit) But I get the following error message: Error in getGroupsFormula.default(correlation, asList = TRUE) : "Form" argument must be a formula I have used this same form

Thankyou for your help with subsetting within a function. I have now tried to apply the same theory in the framework of an lme as follows: fit1.lme <- eval(substitute(lme(fixed=fixed, data=dataframe, random=random, correlation=corCAR1(form= corr), na.action=na.omit, subset=subset),list(subset=subs))) but I get the following error: Error in switch(mode(object), name = , numeric = , call =

hi I am trying to write a function around a glm or similar function, in which I can pass the subsetting constraint as an argument to the function, but I am having trouble. the following commands (if I wanted SEX==0 to be my subset) in the global environment work fine: subexpr <- expression(SEX==0) subtest <- with(mydata, eval(subexpr)) test.glm <- glm(y~x1+x2, data=mydata,

I have installed Samba2.2.7a on an Ultra5 Solaris 8 Operating System. The smb.conf file has been configured and all appears to work fine. However, along with my shared resource, I am seeing the Printers icon. I only want to use Samba to share resources, and not for printing. Is there a way to make the Printers Icon invisible? I already tried the hide files, and veto files parameters in

Hi I'm trying to wrap predict.lm within a function, but I'm having problems passing arguments into it in this way. Basically I want to create a lm object, then pass it into the predict.lm function and be able to tell predict.lm which variable I want to predict for, by passing the variable name as an argument of the wrapper function. This variable will always be a factor with 3 levels,

Greetings from the day of rest after a week of Where. It was a great time out here - and there is a lot of interest in the community around Mapstraction. The new Sandbox (http://mapstraction.appspot.com) made surprisingly quick rounds and lots of kudos. It will be a great place to demo Mapstraction and encourage developers. We had a very good discussion about Mapstraction at WhereCamp. Pamela

Hi, I have been using lm in R to do a linear regression and find the slope coefficients and value for R-squared. The R-squared value reported by R (R^2 = 0.9558) is very different than the R-squared value when I use the same equation in Exce (R^2 = 0.328). I manually computed R-squared and the Excel value is correct. I show my code for the determination of R^2 in R. When I do not set 0 as the

Google Summer of Code has announced the opening of applications for Mentor organizations. Pamela Fox mentioned awhile ago that Mapstraction may be a good fit. http://socghop.appspot.com/ I would like to brainstorm if we have identifiable projects, and preferably even interested students, that would like to work together on putting in some applications. A couple of initial possibilities: -

I am making a browser-based game using ruby on rails and I am having trouble implementing a ranking system for my game. Each player has five skills: intelligence, creativity, strength, charisma and technological. I want the player to have a rank for each skill. The ranking I have in place right now gives the same rank if two players have the same amount of skill. SELECT COUNT(*) AS [rank] FROM

Hi, I need some help to export a data frame with 83 rows and 1411 colums. I have used the package RODBC until now. But now, I have 1411 colums that I can't send to the old excel. If anybody knows a package to convert my data frame in xlsx tell me! Thanks for help in advance! Pamela Botrel -- View this message in context:

Hi there; I am a neophyte to R though I have been messing around with programming in other languages and environments for some years (my dog's name is punchcard to give you some idea of how many years). I have been trying to make a package and install it, to no avail. The functions I have written all work as expected (by me, that is). But I cannot get the html help pages to work. The

Hello to everyone! I have built CRAN task view for genetics. For now I have not submit it to CRAN yet and it can be accessible from: http://www.bfro.uni-lj.si/MR/ggorjan/software/R/Genetics.html http://www.bfro.uni-lj.si/MR/ggorjan/software/R/Genetics.ctv I have not submitted it to CRAN, since I would like first some opinion about it. Genetics is really so broad field that I belive one person

Hi again; I have placed many examples in the CTFS Rd pages because my audience really really likes lots of examples in order to understand the options available in a functions. At this time I have set them all to \dontrun{}. however, this isn't because the example command line with its function and options won't work, but rather that CTFS specific datasets must exist to have the

I'm not looking to start a discussion of the GPL on R-devel; I'm just giving a reference because the subject of the GPL, the license under which R is issued, has been brought up on this list. If you want to follow some of the legal discussion of the Free Software Foundation's General Public License (GPL), the web site groklaw.net is required reading. One reason is because Pamela

I'm needing to remove the R-squared , and put this in a vector. It's the same that I do with the Beta or alfa in a linear regression: beta <- vector() ## beta for (i in 1:ncol(ret)){ beta[i] <- summary(lm((ret[,i] - ret.selic.ibov[,1]) ~ rmrf[,1] + max[,1]))*$coef[[2]]* } So, I need to know what I'll put after* $*. Thanks. -- View this message in context:

Hi All, I'm trying to plot data that is a time series of flows that are associated with a specific level, and I would like each level to represent a colour in a line plot. Here is some data that approximates what I'm using: date=c(1:300) flow=sin(2*pi/53*c(1:300)) levels=c(rep(c("high","med","low"),100)) data=cbind.data.frame(date, flow, levels)

I've been trying to build a Windows binary of a package of mine without success. It seems that the files under inst/doc throw the script off. I am using the command 'Rcmd INSTALL --build'. ---------- Making package genepi ------------ adding build stamp to DESCRIPTION installing NAMESPACE file and metadata installing R files installing inst files FIND: Parameter format not

Is there a way to append to the outfile when using R CMD BATCH? My code, right now, is: R CMD BATCH --slave --vanilla '--args place .2 -.1 .9 .6' StratificationSimulation example.output Everything else is working the way I'd like it. The first few lines of code of my script file are: options(echo=FALSE) cmd_args = commandArgs() print (cmd_args) #d <-

Any users of the VocTel VOIP service? (Canadian) How have you found the quality (Choppy / smooth audio)? Any problems registering? (I have been unable to register for hours) After reading about the collapse of a big USA VOIP provider, I'm curious Thanks, OCG -------------- next part -------------- An HTML attachment was scrubbed... URL:

Hi All, I have a data set of daily measurements of river flow. I would like to create a "ts" object from this data. Here's a sample data set: date <- as.Date(c(1:300), format="%Y") year=as.numeric(format(date, format = "%Y")) month=as.numeric(format(date, format = "%m")) julianday=as.numeric(format(date, format = "%j"))