Displaying 20 results from an estimated 4000 matches similar to: "lme: howto specify covariance structure between levels of grouping factors"
2003 Nov 25
1
using pdMAT in the lme function?
Hello. I want to specify a diagonal structure for the covariance matrix
of random effects in the lme() function.
Here is the call before I specify a diagonal structure:
> fit2<-lme(Ln.rgr~I(Ln.nar-log(0.0011)),data=meta.analysis,
+ random=~1+I(Ln.nar-log(0.0011)|STUDY.CODE,na.action=na.omit)
and this works fine. Now, I want to fix the covariance between the
between-groups slopes
2004 Apr 22
1
lme correlation structure error
Hi there fellow R-users,
I am trying to follow an example of modelling a serial correlation structure
in the textbook "Mixed Effects Model in S and Splus".
However, I am getting some very odd results. Here is what I am trying to
run:
library(nlme)
data(Ovary)
fm1<-lme(follicles~sin(2*pi*Time)+cos(2*pi*Time),data=Ovary,random=pdDiag(~s
in(2*pi*Time)))
### The example is fine up
2005 Sep 19
1
How to mimic pdMat of lme under lmer?
Dear members,
I would like to switch from nlme to lme4 and try to translate some of my
models that worked fine with lme.
I have problems with the pdMat classes.
Below a toy dataset with a fixed effect F and a random effect R. I gave
also 2 similar lme models.
The one containing pdLogChol (lme1) is easy to translate (as it is an
explicit notation of the default model)
The more parsimonious
2024 Jan 08
1
how to specify uncorrelated random effects in nlme::lme()
Dear professor,
I'm using package nlme, but I can't find a way to specify two uncorrelated random effects. For example, a random intercept and a random slope. In package lme4, we can specify x + (x ll g) to realize, but how in nlme?
Thanks!
????????????????????????
Zhen Wang
Graduate student, Department of Medical Statistics, School of Public Health, Sun Yat-sen
2007 Feb 28
5
What is happening here?
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
Hi readers of this list,
I have this nasty problem, that looks fixed, but everytime something
else annoying happens...
- From SuSE10.0, there are problems with browsing LAN, due to
susefirewall, was said...and partialy this is true..
All kind of strange things, incompatibilities, or bugs, keep me from
just browsing shares in the network...
I filed a
2004 Aug 06
2
icecast/ices
Marco,
I have some preliminary HowTos and software available at
http://quasi.ksl.com/icecast/
It's not as nice as http://www.6809.org.uk/kja3/ices2-howto.shtml, but
I'm working to make it better as I become more familiar. Hope it helps.
The ices2 program is designed for .ogg streaming while icecast-1.3.12
and related streamers, MuSE, liveice and DarkIce will still stream .mp3
quite nicely.
2003 May 12
1
plot.ranef.lme (PR#2986)
library(nlme)
data(Phenobarb)
na.include <- function(x)x
phe1 <- nlme(conc~phenoModel(Subject, time, dose, lCl, lV),
data = Phenobarb,
fixed = lCl+lV~1,
random= pdDiag(lCl+lV~1),
start = c(-5,0),
na.action = na.include,
naPattern = ~!is.na(conc))
phe.ranef <- ranef(phe1,augFrame=TRUE)
plot(phe.ranef, form=lCl~Wt+ApgarInd)
[Error in max(length(x0),
2000 Jan 06
1
nlme
Among others, datam contains the columns: logconc, tm, dose, subj, bilirubin.
None of these are factor variables.
The following compartment models work (the first still has not
converged after 100 interations):
res1 <- nlme(logconc~p2+p3+log(dose/(exp(p1)-exp(p2))*
(exp(-exp(p2)*tm)-exp(-exp(p1)*tm))),start=list(fixed=c(5,-2,-0.1)),
fixed=list(p1+p2+p3~1),control=list(maxIter=100),
2010 Jul 11
1
How to automatically restart nlme in a user-defined function
Hi, everyone,
I wrote a function, which includes an nlme estimation. The problem is
sometimes nlme may not converge due to too many random effects.
Say a, b are two parameters. if I specify random effects by:
random = a+b~1,
nlme fails to converge.
Then I have to constrain the random effects in a positive definite diagonal
matrix by:
random = list(pdDiag(a+b~1))
My question is how I can
2011 Mar 17
2
fitting gamm with interaction term
Hi all,
I would like to fit a gamm model of the form:
Y~X+X*f(z)
Where f is the smooth function and
With random effects on X and on the intercept.
So, I try to write it like this:
gam.lme<- gamm(Y~ s(z, by=X) +X, random=list(groups=pdDiag(~1+X)) )
but I get the error message :
Error in MEestimate(lmeSt, grps) :
Singularity in backsolve at level 0, block 1
2012 Feb 05
1
Covariate model in nlme
Dear R users,
I am using nlme to fit a pharmacokinetic model. The base model is
parameterized in terms of CL, V1, V2 and Q.
basemodel<-nlme(Conc ~TwoCompModel(CL,Q,V1,V2,Time,ID),
data = data2, fixed=list(CL+Q+V1+V2~1),
random = pdDiag(CL+V1+V2~1),
start=c(CL=log(20),Q=log(252),V1=log(24.9),V2=log(120)),
control=list(returnObject=TRUE,msVerbose=TRUE,
msMaxIter=20,pnlsMaxIter=20,pnlsTol=1),
2010 Apr 14
3
pdMat
Alguien tiene experiencia en escribir una pdMat. Para aquellos que no lo
recuerden son las matrices de covarianzas de los efectos aleatorios que
ajusta la función lme de la librería nlme
Estas matrices tiene especial importancia en aplicaciones de genética de
poblaciones y en particular en mapeo de asociación. Pinheiro y Bates dicen
que el usuario puede crear sus propias pdMat y sugiere como
2010 Jan 04
1
glmer (lme4), glmmPQL (MASS) and xtmepoisson (Stata)
Dear R users,
I'm trying to specify a generalized linear mixed model in R, basically a Poisson model to describe monthly series of counts in different regions.
My aim is to fit subject-specific curves, modelling a non-linear trend for each region through random effects for linear splines components (see Durban et al, Stat Med 2005, or " Semiparametric regression" by Ruppert et al,
2006 Jul 03
1
analogue of group option of SAS MIXED/random in R
Dear list,
I am trying to use lme to build the analogue of the following SAS MIXED
random specification:
random int+Variable1+Variable2 /subject = Subject group=Condition type=vc;
which gives a Condition-blocked heterogeneity in the random effects
variance-covariance matrix.
Needless to say, I have a hard time in specifying Condition-specific
heterogeneities in the variance-covariance
2008 May 27
3
How to test significant differences for non-linear relationships for two locations
Hi List,
I have to compare a relationship between y and x for two locations. I found logistic regression fits both datasets well, but I am not sure how to test if relationships for both sites are significantly different. I searched the r site, however no answers exactly match the question.
I used Tukey's HSD to compare two means, but the relationship in my study was not simply linear. So I
2004 Jan 21
0
intervals in lme() and ill-defined models
There has been some recent discussion on this list about the value of using
intervals with lme() to check for whether a model is ill-defined. My
question is, what else can drive very large confidence intervals for the
variance components (or cause the error message "Error in
intervals.lme(Object) : Cannot get confidence intervals on var-cov
components: Non-positive definite approximate
2012 Mar 09
0
pdMat class in LME to mimic SAS proc mixed group option? Group-specific random slopes
I would like to be able to use lme to fit random effect models In which some but not all of the random effects are constrained to be independent. It seems as thought the pdMat options in lme are a promising avenue. However, none of the existing pdMat classes seem to allow what I want.
As a specific example, I would like to fit a random intercept/slope mixed model to longitudinal observations in
2013 Oct 26
2
Problems with lme random slope+intercept model
Dear all,
I'm trying to fit a model on ecological data in which I have measured a few
biotic and abiotic factors over the course of a few days in several
individuals. Specifically, I'm interested in modelling y ~ x1, with x2, x3,
and 'factor' as independent variables. Because data suggests both slope and
intercept (for y ~x1) might differ between individuals, I'd want to
2008 Aug 22
1
lme questions re: repeated measures & covariance structure
Hello,
We are attempting to use nlme to fit a linear mixed model to explain bird
abundance as a function of habitat:
lme(abundance~habitat-1,data=data,method="ML",random=~1|sampleunit)
The data consist of repeated counts of birds in sample units across multiple
years, and we have two questions:
1) Is it necessary (and, if so, how) to specify the repeated measure
(years)? As written,
2010 Feb 04
1
random slope models with lme --> failured to converge
Dear all,
I am working on a data set in which I have sequentially measured egg
temperatures ("eggtemp") in birds incubating in different ambient
temperatures ("treat", sample data set below), "id" is not replicated within
treatment.
id treat eggtemp
1 79 3 30.90166
2 42 3 34.94044
3 10 3 32.69945
4 206 3 36.64127
5 23 3 31.80055
6