similar to: Error in NLME (nonlinear mixed effects model)

Dear R-users, I am trying to create a model using the NLS function, such that: Y = f(X) + q + e Where f is a nonlinear (Weibull: a*(1-exp(-b*X^c)) function of X and q is a covariate (continous variable) and e is an error term. I know that you can create multiple nonlinear regressions where x is polynomial for example, but is it possible to do this kind of thing when x is a function with unknown

Hola! I am encountering the following problem, in a multilevel analysis, using glmmPQL from MASS. This occurs with bothj rw1062 and r-devel, respectively with nlme versions 3.1-38 and 3.1-39 (windows XP). > S817.mod1 <- glmmPQL( S817 ~ MIEMBROScat+S901+S902A+S923+URBRUR+REGION+ + S102+S103+S106A+S108+S110A+S109A+S202+S401+S557A+S557B+ + YHOGFcat,

Hi there, I am looking to compare nonlinear mixed effects models that have different nonlinear functions (different types of growth curve)embedded. Most of the literature I can find focuses on comparing nested models with likelihood ratios and AIC. Is there a way to compare model fits when models are not nested, i.e. when the nonlinear functions are not the same? Many thanks in advance! Lindsay

Dear all, I am trying to fit a mixed-effects non linear regression, but I have some trouble with it. My data are a balanced panel of 904 subjects with 8 observations (at regular periods) per subject. The functional form of my model is Y=Aexp(-BX1)X2 +e. I want to allow parameters A and B to vary among subjects and also include an autocorrelation term. I have already fitted a standard nonlinear

Subject: nls, step factor 0.000488281 reduced below 'minFactor' of 0.000976563 Hi there, I'm trying to conduct nls regression using roughly the below code: nls1 <- nls(y ~ a*(1-exp(-b*x^c)), start=list(a=a1,b=b1,c=c1)) I checked my start values by plotting the relationship etc. but I kept getting an error message saying maximum iterations exceeded. I have tried changing these

Hello Timo, With dovecot 2.1.rc1 (056934abd2ef) there seems to be something wrong with virtual plugin mailbox search patterns. I'm using a virtual mailbox 'unread' with the following dovecot-virtual file $ cat dovecot-virtual * unseen For testing propose I created the following folders with each containing one unread message INBOX, INBOX/level1 and INBOX/level1/level2

>> One clear moral seems to be don't do anything more inside a >> generic function than you really need to do. Keep it *very* >> simple indeed. >> > I recall JMC saying something like, all generic functions > should be one line long; a call to the appropriate UseMethod. It certainly is encouraging to know that others also have been confused by aspects of

Hello all, In lme4 if you want to model two non-nested random effects you code it like this: mod1 <- lmer(y~x + (1|randomvar1) + (1|randomvar2)) How would you go about to model something similar in nlme? In my database I have two variables for which I have repeated measures, lets call them "individual" and "year". But none of the "individuals" were measured in

Hi R Users, When I use package nlme for linear model with random effects, there exists errors and I don't know the data structure of lme. Here is my data: Rice<-data.frame(Yield=c(8,7,4,9,7,6,9,8,8,8,7,5,9,9,5,7,7,8,8,8,4,8,6,4,8,8,9), Variety=rep(rep(c("A1","A2","A3"),each=3),3),

Hello, I have an unbalanced mixed model design with two fixed effects "site" (2 levels) and "timeOfDay" (4 levels) and two random effects "day" (3 consecutive days) and "trap" (6 unique traps, 3 per site). The dependent variable is the body length ("BL") of insect larvae from 7 to 29 individuals per trap (104 individuals in total). To account

Hello list, A change in smbclient between 3.0.24 and 3.0.25c breaks [1] the gollem filemanager [2] and maybe others. In particular folder creation is broken unless the folder is a top level folder on the share. To create a subfolder "level2" in the folder /level1 of a share gollem executes this command: /usr/bin/smbclient "//1.2.3.4/sharename" "-p139"

I have some matrices stored as elements in a list that I am working with. On example is provided below as TP[[18]] > TP[[18]] level2 level1 1 2 3 4 1 79 0 0 0 2 0 0 0 0 3 0 0 0 0 4 0 0 0 0 Now, using prop.table on this gives > prop.table(TP[[18]],1) level2 level1 1 2 3 4 1 1 0 0 0 2 3

Hi experts I have a trouble in access rights I am running Samba 3.0.31 on Solaris 10 x86 64 bits as member server of an Active Directory 2003 R2 domain (MYDOMAIN) using Identity Management for Unix I set rights to access a sub folder of a Samba share. On Solaris the user "toto" jdoe can write a new file. From Windows, the same user can't. Itlooks like OK when the primary group

On 26/01/13 19:20, Torsten Rüger wrote: [...] > Just at the end it mentions closures and I was wondering how those are done in llvm. > The link was to wikipedia, and i do know what closures/blocks/continuations are, (i think) but maybe someone could point me to where to read about how to do them in llvm. I've implemented closures in the past --- it's fiddly and a surprising amount

Dear R experts, I wonder how can one input a string variable in multiple lines in a R script. I've seen solution to the command line continuation in the non-string situation (just continute at the new line). However, that doesn't work for strings because the new line character, "\n", is included in the resultant string. For example, if I type the following in a script:

>From my investigations, and searching the archives, it seems that smbfs can only mount shares, and not directories under the share i.e. smbmount //host/share /mnt works smbmount //host/share/directory /mnt doesn't work Our W2K3 server has user home directories in a tree where there are no user rights to the next node up the tree (level1), which is the directory on the share, i.e.

Hello, I''ve done potencial replacement for RoR''s routing mechanism. It just works as an independent class so you can try bundled examples withou installing anything else. I didn''t integrated it into Rails (as a plugin?) because I don''t know if anyone would appreciate it etc. SiteMapper features: -------------------- 1. uses regexp mask to get variables from

Hi, the following R lines work fine in R 2.4.0, but not in R 2.4.1 or any devel versions of R 2.5.0 (see below for details). library(drc) # to load the dataset 'PestSci' library(nlme) ## Setting starting values sv <- c(0.43355869, 2.49963220, 0.05861799, 1.73290589, 0.38153146, 0.24316978) ## No error m1 <- nlme(SLOPE ~ c + (d-c)/(1+exp(b*(log(DOSE)-log(e)))), fixed =

Dear R users, Does anyone knows how to run a glmm with one fixed factor and 2 random numeric variables (indices)? Is there any way to force in the model a separate interaction of those random variables with the fixed one? I hope you can help me. #eg. Reserve <- rep(c("In","Out"), 100) fReserve <- factor(Reserve) DivBoulders <- rep

Hi Guys, is there a way to have a absolute file paths in a file - and use --files-from together with --delete option? what i mean is - lets say i have a daily generated file a.txt from an external script which contains something like this: /level1/level2/level/.../.../leveln/file.txt can i ask rsync to copy all files from that kind - /level1/level2/level/.../.../leveln/file.txt which my file