[R] convergence

2007-04-19 Thread rach.s
hie.. how can i write a loop that makes algorithm keeps repeating until a solution is converged?do i use a for loop? i know that we can use for loop to ask for a number of repetitions, but how to use it to ask the algorithm to keep repeating until a solution is converged? Thanks -- View this mess

Re: [R] convergence

2007-04-19 Thread Ted Harding
On 19-Apr-07 12:00:17, rach.s wrote: > > hie.. > how can i write a loop that makes algorithm keeps repeating > until a solution is converged?do i use a for loop? i know > that we can use for loop to ask for a number of repetitions, > but how to use it to ask the algorithm to keep repeating > until

Re: [R] convergence

2007-04-19 Thread Doran, Harold
8:00 AM > To: r-help@stat.math.ethz.ch > Subject: [R] convergence > > > hie.. > how can i write a loop that makes algorithm keeps repeating > until a solution is converged?do i use a for loop? i know > that we can use for loop to ask for a number of repetitions, > but how

Re: [R] convergence

2007-04-19 Thread Jim Lemon
rach.s wrote: > hie.. > how can i write a loop that makes algorithm keeps repeating until a solution > is converged?do i use a for loop? i know that we can use for loop to ask for > a number of repetitions, but how to use it to ask the algorithm to keep > repeating until a solution is converged? g

Re: [R] convergence

2007-04-19 Thread Alberto Monteiro
Ted Harding wrote: > > There are various ways round this, but a 'for' loop with > a fixed number of iterations is not usully one of them! > > The simplest is to use while(). A possibly strategy is > > Y.old <- initial.Y > while(TRUE){ > Y <- compute.Y(Y.old, ...) > if(abs(Y - Y.old)

Re: [R] convergence

2007-04-19 Thread Ravi Varadhan
, April 19, 2007 9:24 AM To: [EMAIL PROTECTED]; r-help@stat.math.ethz.ch Subject: Re: [R] convergence Ted Harding wrote: > > There are various ways round this, but a 'for' loop with > a fixed number of iterations is not usully one of them! > > The simplest is to use while(

Re: [R] convergence

2007-04-20 Thread Julien Barnier
Hi, > how can i write a loop that makes algorithm keeps repeating until a solution > is converged?do i use a for loop? i know that we can use for loop to ask for > a number of repetitions, but how to use it to ask the algorithm to keep > repeating until a solution is converged? I don't know what

[R] convergence in polr

2004-02-24 Thread C. Spanou
Hello splus-users, I am trying to fit a regression model for an ordered response factor. So I am using the function polr in library(MASS). My data is a matrix of 1665 rows and 63 columns (one of the column is the dependent variable). The code I use is polr(as.ordered(q23p)~.,data=newdatap) but I

[R] Convergence in geese/gee

2006-08-10 Thread Sharon Kuhlmann-B
We are currently analyzing data on children clustered in day care-centers (DCC). We have tried to use geepack and gee libraries to estimate an overall incidence rate for absences (=number of absences/risk time) by specifying geese(number.absences ~ offset(log(risktime)), id=day.care.id,

[R] Convergence problem in GLMM

2004-11-23 Thread John Fox
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 ~ a

[R] Convergence problems with survreg()

2005-02-01 Thread Stefano Conti
Dear R mailing list, I'm trying to fit a censored regression model to a large (dimension of the design matrix is 2e5 by 7) right truncated data by means of the survreg(Surv()) function, as suggested by Paul Johnson on his "R Tips" web page. Possibly due to the sensitivity to the initial values

Re: [R] convergence in polr

2004-02-24 Thread Prof Brian Ripley
Why have you sent a message about S-PLUS to R-help, one that has already been answered on S-news? There is no function nlminb in R. On 24 Feb 2004, C. Spanou wrote: > Hello splus-users, I am trying to fit a regression model for an ordered > response factor. So I am using the function polr in l

Re: [R] convergence in polr

2004-02-24 Thread C. Spanou
I am really sorry. I was supposed to send it to the Splus users but by mistake I sent to the R-users. Sorry once again On Feb 24 2004, Prof Brian Ripley wrote: Why have you sent a message about S-PLUS to R-help, one that has already been answered on S-news? There is no function nlminb in R.

[R] convergence for proportional odds model

2005-09-05 Thread liu abc
Hey, everyone, I am using proportional odds model for ordinal responses in dose-response experiments. For some samll data, SAS can successfully provide estimators of the parameters, but the built-in function polr() in R fails. Would you like to tell me how to make some change so I can use polr

[R] convergence for proportional odds model

2005-09-06 Thread David Duffy
liu abc <[EMAIL PROTECTED]> wrote: > > I am using proportional odds model for ordinal responses in > dose-response experiments. For some samll data, SAS can successfully > provide estimators of the parameters, but the built-in function polr() > in R fails. Would you like to tell me how to make some

[R] convergence of coxfilter and coxph

2007-05-21 Thread carol white
Hi, coxfilter function in genefilter package uses coxph to fit a model to filter genes. how come that coxfilter could converge to find a solution in cox model fitting using a data matrix of 8000 variables and 600 samples but coxph doesn't converge with the same matrix? regards, carol ---

[R] Convergence code in nlm function

2004-08-21 Thread Victoria Landsman
Dear R users, I am using the nlm function for minimization of the very non-linear function of four parameters. I am running 100 simulations and almost always I get the convergence code =2 (Successive iterates within tolerance. Current iterate is probably solution.) [about 75 times of 100]. Fre

Re: [R] Convergence problem in GLMM

2004-11-23 Thread Deepayan Sarkar
On Tuesday 23 November 2004 11:14, John Fox wrote: > 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

RE: [R] Convergence problem in GLMM

2004-11-23 Thread John Fox
nal Message- > From: Deepayan Sarkar [mailto:[EMAIL PROTECTED] > Sent: Tuesday, November 23, 2004 5:10 PM > To: [EMAIL PROTECTED] > Cc: John Fox; 'R-Help' > Subject: Re: [R] Convergence problem in GLMM > > On Tuesday 23 November 2004 11:14, John Fox wrote: &g

Re: [R] Convergence problem in GLMM

2004-11-23 Thread Richard A. O'Keefe
I was trying to install some more packages and ran into a problem I hadn't seen before. Version: platform sparc-sun-solaris2.9 arch sparc os solaris2.9 system sparc, solaris2.9 status major2

[R] convergence diagnostic test in R

2003-07-30 Thread GUNES Evrim-Didem
Hello, I am trying to perform a diagnostic test using BOA in R. I get the following error message for Heidelberger and Welch test. - Error in if (I[i] < 0.46) { : missing value where logical needed I am new to this environment, but I suspect this might be a bug. Does anyone have an idea? T

[R] Convergence warnings from zeroinfl (package pscl)

2006-07-20 Thread Daniel E. Bunker
Dear R-Helpers, Can anyone please help me to interpret warning messages from zeroinfl (package pscl) while fitting a zero inflated negative binomial model? The console reports convergence and the parameters seam reasonable, but these <> suggest that it did not converge. (See full output belo

Re: [R] convergence of coxfilter and coxph

2007-05-22 Thread Thomas Lumley
On Mon, 21 May 2007, carol white wrote: > Hi, coxfilter function in genefilter package uses coxph to fit a model > to filter genes. how come that coxfilter could converge to find a > solution in cox model fitting using a data matrix of 8000 variables and > 600 samples but coxph doesn't converge

[R] convergence error (lme) which depends on the version of nlme (?)

2005-12-12 Thread Leo Gürtler
Dear list members, the following hlm was constructed: hlm <- groupedData(laut ~ design | grpzugeh, data = imp.not.I) the grouped data object is located at and can be downloaded: www.anicca-vijja.de/lg/hlm_example.Rdata The following works: library(nlme) summary( fitlme <- lme(hlm) ) with out

Re: [R] convergence error (lme) which depends on the version of nlme (?)

2005-12-12 Thread Dieter Menne
Leo Gürtler anicca-vijja.de> writes: > > > hlm <- groupedData(laut ~ design | grpzugeh, data = imp.not.I) > > the grouped data object is located at and can be downloaded: > www.anicca-vijja.de/lg/hlm_example.Rdata ... > > > 2) What is the cause of the convergence error which seems to d

Re: [R] convergence error (lme) which depends on the version of nlme (?)

2005-12-13 Thread Douglas Bates
On 12/12/05, Leo Gürtler <[EMAIL PROTECTED]> wrote: > Dear list members, > > the following hlm was constructed: > > hlm <- groupedData(laut ~ design | grpzugeh, data = imp.not.I) > > the grouped data object is located at and can be downloaded: > > www.anicca-vijja.de/lg/hlm_example.Rdata > > The fo