Displaying 6 results from an estimated 6 matches for "getgroupsformula".
2006 May 17
1
nlme model specification
Hi folks,
I am tearing my hair out on this one.
I am using an example from Pinheiro and Bates.
### this works
data(Orange)
mod.lis <- nlsList(circumference ~ SSlogis(age, Asymp, xmid, scal),
data=Orange )
### This works
mod <- nlme(circumference ~ SSlogis(age, Asymp, xmid, scal),
data=Orange,
fixed = Asymp + xmid + scal ~ 1,
start =
2003 Apr 14
1
Problem with nlme or glmmPQL (MASS)
...+REGION+
+
S102+S103+S106A+S108+S110A+S109A+S202+S401+S557A+S557B+
+ YHOGFcat, data=Hogares, random= ~1,
family=binomial,
+ na.action=na.omit)
iteration 1
Error in getGroups.data.frame(dataMix, groups) :
Invalid formula for groups
> getGroupsFormula(Hogares)
~ID08
Hogares is a groupedData object. The only non-default thing I have
done is to use order.groups=FALSE in groupedData, since the response
is a factor and the max ordering function does not apply.
What is wring here?
Kjetil Halvorsen
2005 Jun 10
1
Problems with corARMA
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 argument with differerent correlation
structures and it has worked fine. Does anyone know why it won't
recognise ~year | id2 (or even ~ 1 | id2) as a formula?
Any help would...
2006 Mar 07
1
lme and gls : accessing values from correlation structure and variance functions
Dear R-users
I am relatively new to R, i hope my many novice questions are welcome.
I have problems accessing some objects (specifically the random effects, correlation structure and variance function) from an object of class gls and lme.
I used the following models:
yah <- gls (outcome~ -1 + as.factor(Trial):as.factor(endpoint)+
1997 Apr 28
0
R-beta: a recursive anonymous function
...;call") {
if (x[[1]] == as.name("/"))
return(c(sys.function()(x[[2]]), sys.function()(x[[3]])))
if (x[[1]] == as.name("(")) # for R
return(sys.function()(x[[2]]))
}
if (mode(x) == "(") return(sys.function()(x[[2]])) # for S
list(x)
})(getGroupsFormula(data, form, ...)[[2]])
## I know it's horribly obscure. Blame Bill Venables for teaching me this.
Regretably, it doesn't work in R. Using the debugger one finds that
sys.function() returns the function being called the first time
through but the second time through it returns NULL. Is...
2002 Feb 12
1
problem plotting nls objects: couldn't find function "..."
...== cF))
ylab <- cL
else if (!is.null(rF) && (ylab == rF))
ylab <- rL
}
argForm <- .y ~ .x
argData[, ".y"] <- .y
if (is.null(args$ylab))
args$ylab <- ylab
}
grpsF <- getGroupsFormula(form)
if (!is.null(grpsF)) {
gr <- splitFormula(grpsF, sep = "*")
for (i in 1:length(gr)) {
auxGr <- all.vars(gr[[i]])
for (j in auxGr) {
argData[[j]] <- eval(as.name(j), data)
}
}
if (length(...