Displaying 20 results from an estimated 500 matches similar to: "lme, lmer, convergence"
2011 Nov 09
1
Invalid 'yscale' in viewport
I tried Googling for help on this problem, but what I saw only increased my
puzzlement.
I found a lovely bit of code called esplot() that makes scatterplots with
associated histograms and rug plots. I have been trying to use it on my own
data, but the viewport routine on which it depends (I think) is choking on the
data I'm trying to plot.
$ Error in valid.viewport(...) : Invalid
2012 Mar 31
2
unable to move temporary installation
Hi all,
I'm having a strange error that prevents me from installing new packages,
or updating packages after reinstalling. The error message is
Warning: unable to move temporary installation ‘C:\Program
Files\R\R-2.14.2\library\file15045004ac2\sandwich’ to ‘C:\Program
Files\R\R-2.14.2\library\sandwich’
for one of the packages that is failing to install/update. This error
started happening
2006 Jun 28
3
lme convergence
Dear R-Users,
Is it possible to get the covariance matrix from an lme model that did not converge ?
I am doing a simulation which entails fitting linear mixed models, using a "for loop".
Within each loop, i generate a new data set and analyze it using a mixed model. The loop stops When the "lme function" does not converge for a simulated dataset. I want to
2006 May 26
2
lme, best model without convergence
Dear R-help list readers,
I am fitting mixed models with the lme function of the nlme package.
If I get convergence depends on how the method (ML/REM) and which (and
how much) parameters will depend randomly on the cluster-variable.
How get the bist fit without convergence?
I set the parameters msVerbose and returnObject to TRUE:
lmeControl(maxIter=50000, msMaxIter=200, tolerance=1e-4,
2008 Jun 14
1
"False convergence" in LME
I tried to use LME (on a fairly large dataset, so I am not including it), and I got this error message:
Error in lme.formula(formula(paste(c(toString(TargetName), "as.factor(nodeInd)"), :
nlminb problem, convergence error code = 1
message = false convergence (8)
Is there any way to get more information or to get the potentially wrong estimates from LME?
(Also, the page in the
2003 Jan 18
4
File Creation Dates Question
Hi all,
I searched the archives for a solution to this issue
but didn't find any thing.
I have an office with a mix of PCs running ME and XP
that use Samba 2.2.7 under RH Linux 7.3 as an archive
for mortgage information.
When documents are moved to the Samba share, the file
creation timestamp on the Linux box reflects the time the files
were created on the PCs. What I need is for the
2009 Jan 22
1
convergence problem gamm / lme
Hope one of you could help with the following question/problem:
We would like to explain the spatial
distribution of juvenile fish. We have 2135 records, from 75 vessels
(code_tripnr) and 7 to 39 observations for each vessel, hence the random effect
for code_tripnr. The offset (‘offsetter’) accounts for the haul duration and
sub sampling factor. There are no extreme outliers in lat/lon. The model
2005 Sep 06
3
Samba Over A Wan
Hi All,
I have a customer that I need to give access to one of our servers over a WAN.
Their box is running some windows server version and they use terminal services
to let their employees access their applications.
We use Samba in or local office environment where it works quite well. I need
to know the following:
1) Will Samba work over the internet?
2) Can Samba connections be encrypted or
2007 Feb 28
1
LME without convergence
Dear R-help list readers,
I am fitting a mixed model using the lme function (R V 2.3.1 for
Windows). This is an example:
dep<-c(25,40,33.33,60,70.83,72,71.43,50,40,53.33,64,54.17,60,53.57)
yes<-c(0,1,2,3,4,5,6,0,1,2,3,4,5,6)
treat<-c(1,1,1,1,1,1,1,0,0,0,0,0,0,0) #factor
If I now fit a model with random slopes as well as intercepts:
model1<-lme(dep~yes,random=yes|treat)
R
2011 Jun 22
1
lme convergence failure within a loop
Hi R-users,
I'm attempting to fit a number of mixed models, all with the same
structure, across a spatial grid with data points collected at various
time points within each grid cell. I'm trying to use a 'for' loop to try
the model fit on each grid cell. In some cells lme does not converge,
giving me the error:
Error message: In lme.formula(logarea ~ year + summ_d, data =
2004 Nov 23
2
Convergence problem in GLMM
Dear list members,
In re-running with GLMM() from the lme4 package a generalized-linear mixed
model that I had previously fit with glmmPQL() from MASS, I'm getting a
warning of a convergence failure, even when I set the method argument of
GLMM() to "PQL":
> bang.mod.1 <- glmmPQL(contraception ~ as.factor(children) + cage + urban,
+ random=~as.factor(children) + cage +
2012 Feb 20
1
prediction for linear mixed model
Hi, I am wondering if we can make prediction on a linear mixed model by lmer() from lme4 package? Specifically I am fitting a very simple glmer() with binomial family distribution, and want to see if I can get the predicted probability like that in regular logistic regression?
fit<-glmer(y~x+(1|id),dat,family=binomial)
where y is the response variable (0, 1), and x is a continuous variable.
2006 Feb 15
1
no convergence using lme
Hi. I was wondering if anyone might have some suggestions about how I can
overcome a problem of "iteration limit reached without convergence" when
fitting a mixed effects model.
In this study:
Outcome is a measure of heart action
Age is continuous (in weeks)
Gender is Male or Female (0 or 1)
Genotype is Wild type or knockout (0 or 1)
Animal is the Animal ID as a factor
2011 Jul 25
1
lme convergence error
Hello, I am working from a linux 64 machine on a server with R-2.12 (I can't
update to 2.13). I am iterating through many linear mixed models for
longitudinal data and I occasionally receive the following convergence
error:
> BI.lme <- lme(cd4 ~ time + genBI + genBI:time + C1 + C2 + C11 + C12,
random =~ 1 + time | IID, data = d)
Error in lme.formula(cd4 ~ time + genBI + genBI:time +
2004 Sep 22
1
Kernel Oplocks and FC2
Hi All,
I am running Samba version 3.0.7-2.FC2.
I've been noticing the following messages in /var/log/messages:
Sep 22 09:35:47 www kernel: lease broken - owner pid = 6024
Sep 22 09:36:32 www kernel: lease broken - owner pid = 25134
I've been told that if I place "kernel oplocks = false" in my
smb.conf file, I won't see this message any more.
However, if I do this, does
2006 Jan 03
1
Raw Hangup messages with IAX2?
Hi All,
I am running asterisk 1.2. I have a softphone connecting from a coworkers
home through their router using IAX2 through our router at the office. Both
have port 4569 for TCP and UDP opened and forwarded to the right pc
and server.
I'm seeing Raw Hangup <person's IP address>, src=0. dst=10787
messages show up in the log like 10 every 5 seconds. We can still
make calls and
2004 Jan 14
2
Logging Detail
Hi all,
I am using Samba 2.2.7-security-rollup-fix on a few RH7.3
systems.
I need to configure the logging to show what user/pc is accessing
what file.
I've been through the How-Tos at Samba.org and read the section
on configuring logging detail but either I'm too tired or too stupid
(both viable options ;-) ) to figure out which level I need to get this
information.
Can anyone tell
2012 Feb 17
2
lmer - error message
Hi all,
I am fairly new to mixed effects models and lmer, so bear with me.
Here is a subset of my data, which includes a binary variable (lake (TOM or
JAN)), one other fixed factor (Age) and a random factor (Year).
lake FishID Age Increment Year
1 TOM 1 1 0.304 2007
2 TOM 1 2 0.148 2008
3 TOM 1 3 0.119 2009
4 TOM 1 4 0.053 2010
5
2006 Jul 23
3
Making a patch
Dear R developers,
is there a preferred format or strategy for making a patch to
contribute to a package that is maintained by R-core? Berwin Turlach
and I have written a very minor extension to lmeControl to allow it to
pass an argument to nlminb for the maximum number of evaluations of
the objective function. I've edited the nlme/R/lme.R and
nlme/man/lmeControl.Rd files. I can diff the
2005 Sep 29
2
how to fix the level-1 variances in lme()?
Dear all,
Edmond Ng (http://multilevel.ioe.ac.uk/softrev/reviewsplus.pdf) provides
an example to fit the mixed effects meta-analysis in Splus 6.2. The
syntax is:
lme(fixed=d~wks, data=meta, random=~1|study, weights=varFixed(~Vofd),
control=lmeControl(sigma=1))
where d is the effect size, study is the study number, Vofd is the
variance of the effect size and meta is the data frame.