similar to: gnls : Rho

Dear R-devel members, While writing a new correlation structure similar to 'corCompSymm' and intended to be used with 'gnls', I got puzzled with the 'Initialize' method. Using 'Initialize' before 'gnls' may be regarded as a mean to set an initial value for the corStruct parameter. However 'gnls' does not work properly with a

Hi, What is the reason for getting two different max of a vector from two functions max() and summary()? Here is an example, > set.seed(2222) > x <- sample(x=1:100000, size=10000, replace=T) > max(x) [1] 99992 > summary(x) Min. 1st Qu. Median Mean 3rd Qu. Max. 6 24510 49950 50020 75040 99990 Which value is the correct one? I had this problem while I

Hello List, I am trying to install RODBC package in a debian linux box but getting the following message. Can anyone help me to find what I am doing wrong here: $ R CMD INSTALL RODBC_1.0-4.tar.gz ### * Installing *source* package 'RODBC' ... checking for gcc... gcc checking for C compiler default output... a.out checking whether the C compiler works... yes checking whether we are cross

Dear all, How can I extract the total and residual d.f. from a gnls object? I have tried str(summary(gnls.model)) and str(gnls.model) as well as gnls(), but couldn?t find the entry in the resulting lists. Many thanks! Best wishes Christoph -- Dr. rer.nat. Christoph Scherber University of Goettingen DNPW, Agroecology Waldweg 26 D-37073 Goettingen Germany phone +49 (0)551 39 8807 fax +49

>>>>> "VW" == Viechtbauer Wolfgang (STAT) <Wolfgang.Viechtbauer at STAT.unimaas.nl> >>>>> on Fri, 26 Sep 2008 18:00:19 +0200 writes: VW> Hi all, I'm trying to fit a marginal (longitudinal) VW> model with an exponential serial correlation function to VW> the Orange tree data set. However, R crashes frequently VW>

Hi there, I have four panels in a lattice bwplot. I want to have two different ylim for the panels, for example panels [1,1] and [1,2] with ylim=c(0,200) and panels [2,1] and [2,2] with ylim=c(0,100). Thanks for help in advance. Mahbub.

Hi, I want to use (similar to) las options in lattice (bwplot) plot. Actually I want to have x-axis labels as vertical instead of default horizontal. Thanks in advance for your help. Mahbub

Hi, How can I use portrait/landscape option in lattice bwplot? Is there any option in trellis.device where I can define this? Thanks in advance, Mahbub.

Hi, I was trying to install XML package in a linux (dabian) machine and got the following error message. I am not sure whether there is any error in my linux installion. I appreciate suggestions to install XML properly in this machine. Thanks in advance. Mahbub. ######Errors#### $ R CMD INSTALL -l /usr/local/lib/R/library/ XML_0.93-4.tar.gz * Installing *source* package 'XML' ...

Dear all, I am trying to gain some experience with the function gnls from the nlme package. I tried to model the Theophyline data by trying to model the presumed dependency of the clearance on the body weight. This is my function call of gnls: gnls(conc~SSfol(Dose,Time,lKe,lKa,lCl),data=Theoph, params=list(lKe~1,lKa~1,lCl~Wt),start=c(-2.4,0.46,-3.22,0.01)) That's been the result: Error

Hi, I am looking for C/C++ codes for computing generalized inverse of a matrix. Can anyone help me in this regard? Thanks, Mahbub.

This works fine: fit42<-nls(Vfs~SSlogis(Months,Asym.Int+Asym.Group*Groupdum,xmid,scal), data=df, start=c(Asym.Int=22,Asym.Group=5,xmid=2,scal=6), na.action=na.omit) But this, identical except using gnls, doesn't converge: fit43<-gnls(Vfs~SSlogis(Months,Asym.Int+Asym.Group*Groupdum,xmid,scal), data=df, start=c(Asym.Int=22,Asym.Group=5,xmid=2,scal=6), na.action=na.omit) Error in gnls(Vfs

Hi all, I'm working with a dataset from 10 treatments, each treatment with 30 subjects, each subject measured 5 times. The plot of the dataset suggests that a 3-parameter logistic could be a reasonable function to describe the data. When I try to fit the model using gnls I got the message 'Step halving factor reduced below minimum in NLS step'. I´m using as the initial values of the

Back on January 16, a message on R-help from Ravi Varadhan described a problem with gnls using weights=varPower(). The problem was that the fit failed with error Error in eval(expr, envir, enclos) : Object "." not found I can reliably get this error in version 2.0.1-patched 2004-12-09 on Windows XP and 2.0.1-Patched 2005-01-26 on Linux. The key feature of that example is that the

I have been using the nls function to fit some simple non-linear regression models for properties of graphite bricks to historical datasets. I have then been using these fits to obtain mean predictions for the properties of the bricks a short time into the future. I have also been calculating approximate prediction intervals. The information I have suggests that the assumption of a normal

Hi All, I have question about speficifying a constant in gnls() from package nlme. Here is a testing code: ############# library(nlme) x = exp( rnorm(100)) y = 1/(1+x) + rnorm(100)/10 plot( y ~ x) fm1 = gnls( y ~ 1/(1+(x/v)^w), start=list( v=1, w=1)) a =1; b=1; fm2 = gnls( y ~ a/(b+(x/v)^w), start=list( v=1, w=1)) #This won't work because I don't know to set $a$ and $b$ as

Dear R-Help, I'm trying to fit a model of the following form using gnls. I've fitted it using nlsList with the following syntax: nlsList(Y~log(exp(a0-a1*X)+exp(b0-b1*X))|K,start=list (a0=6,a1=0.2,b0=4.5,b1=0.001),data=data.frame(Y=y,X=X,K=k))) which works just fine: <snip> Coefficients: a0 a1 b0 b1 1 5.459381 0.5006811 5.137458 -0.0040548687

Last week (Wed 9/10/2003, "regression questions") I posted a question regarding the use of gnls( ) and its dissimilarity to the syntax that nls( ) will accept. No one replied, so I partly answered my own question by constructing indicator variables for use in gnls( ). The code I used to construct the indicators is at the end of this email. I do have a nagging, unanswered

Dear friends There is a list of arrays comprising different no of rows and columns even sometimes NULL, such as [[2]] given below. How can we ignore [[2]] or others like this in the complete list. Any help in this regard is needed. Thanks [[1]] [,1] [,2] [1,] 3 1 [2,] 3 1 [3,] 3 1 [[2]] NULL [[3]] [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 3 1

Dear R users; I was trying to fit a nonlinear model using gnls (nlme version 3.1-80, R 2.5.0, WinXP) and I got the following error and warning message: Error in gnls(ht ~ a1 * hd * (1 - a2 * exp(-a3 * (dbh/dq2))), data = hdat, : Step halving factor reduced below minimum in NLS step In addition: Warning message: $ operator is deprecated for atomic vectors, returning NULL in: