similar to: multiple fitted curves plot

Hi, I have a data set of observations by either one person or a pair of people. I want to only keep the pair observations, and was using the code below until it gave me the error " $ operator is invalid for atomic vectors". I am just beginning to learn R, so I apologize if the code is really rough. Basically I want to keep all the rows in the data set for which the value of

Dear R-list members, Cynthia M. Jones wrote a paper (Fitting growth curves to retrospective size-at-age data, Fisheries Research 46(2000):123-129; abstract at http://www.elsevier.nl/gej-ng/10/19/44/70/24/37/abstract.html)where the SAS procedure MIXED, Macro NLINMIX (Littell et. al., 1996)was used to estimate the von Bertalanffy growth function parameters assuming that data from the same fish are

Dear list, I am faced with the following model: y=E+P+M+H+PxE+Error y is a response, E and P are factors with fixed effects. M is a random effect nested in P and H is a random effect nested in M. PxE is interaction of P and E. It seems that I should fit such a model with the lme function of library nlme, but I was not able to figure out from the help page how to specify the formula. In the

Greetings- I have some data on which I need to do something like a Hierarchical Linear Model (please bear with me, I'm only learning the technique so I don't know yet if my language is correct). Essentially I'm analyzing data at two levels simultaneously; data are about individuals an organizations of which they are members. Can someone point me toward an appropriate package in R?

Hi All. I'm trying to run nlsList an getting an error that makes no sense to me. I have accuracy and reaction time data over many trials for each person (id) When I use nlsList code that is virtually identical to the example in the doc file I get the following error. I've tried everything I could think of and can't get around it. Any ideas what I'm doing wrong? **************

Hi, I have a program written in R which is good on the version 1.2, but for the fallowing versions of R, an error always is at the same place. That is at the level of the fallowing line: Sur<- getInitial(res2[m:M,2]~SSasymp(res2[m:M,1],Asymp,resp0,lrc),data=res2) Error in eval(expr,envir,enclos):numeric envir arg not of length one I don't know at all this langage for the instant.

Hola a todos Les consulto por un problema con xtable y nlme, tomando un ejemplo del manual de nlme para obtener los resultados en latex utilizando xtable, se puede utilizar el siguiente código, pero hay un problema y causa error. library(nlme) library(xtable) fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc), data = Loblolly, fixed = Asym + R0 + lrc ~ 1,

I ran into a problem trying to make a plot from a file that's read using source. Basically, I have the following code in a file "plot.R" : library(nlme) data(Loblolly) fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc), data = Loblolly, fixed = Asym + R0 + lrc ~ 1, random = Asym ~ 1, start = c(Asym = 103, R0 = -8.5, lrc = -3.3))

Better posted on r-sig-mixed-models , no? Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Aug 23, 2017 at 5:17 AM, PIKAL Petr <petr.pikal at precheza.cz> wrote: > Dear all > > I encountered strange

Dear all I encountered strange behaviour of augPred with virtually the same data First I made groupedData object. > mar.g<-groupedData(rutilizace~doba|int, data=mar) When I perform nlme on complete dataset I get an error with augPred > fit<-nlsList(rutilizace~SSasymp(doba, Asym, R0, lrc), data=mar.g) Warning message: c("1 error caught in nls(y ~ cbind(1 - exp(-exp(lrc) * x),

> On Aug 23, 2017, at 8:08 AM, PIKAL Petr <petr.pikal at precheza.cz> wrote: > > Hi > > Well, yes I tried it about two weeks ago but my post did not get through as it still awaits moderator approval. It got through just fine. It appeared on Aug 15. It just didn't get any replies. As I read your original question in this thread, it was not clear to me that you had

Hi Well, yes I tried it about two weeks ago but my post did not get through as it still awaits moderator approval. I could check which column is offending but actually it is only minor nuisance, I can live with selection of columns before fitting a model. What seems to me strange is that both full dataset and only selected colums gave me identical fit results but only one works within augPred.

I need your help to make R works on my linux box, runing fedora 11. Few things on the machine and steps I did: 1. uname -a shows: Linux bagvapp 2.6.29.4-167.fc11.i586 #1 SMP Wed May 27 17:14:37 EDT 2009 i686 athlon i386 GNU/Linux 2. it has a gcc (version 4.4) RedHat's build 3. the g77 comes with the box does not work (R's configure tells it could not compile simple fortran code); so

Hi I would like to fit a non-linear regression to the follwoing data: quantiles<-c(seq(.05,.95,0.05)) slopes<-c( 0.000000e+00, 1.622074e-04 , 3.103918e-03 , 2.169135e-03 , 9.585523e-04 ,1.412327e-03 , 4.288103e-05, -1.351171e-04 , 2.885810e-04 ,-4.574773e-04 , -2.368968e-03, -3.104634e-03, -5.833970e-03, -6.011945e-03, -7.737697e-03 , -8.203058e-03, -7.809603e-03, -6.623985e-03,

Am Dienstag, 25. April 2017, 11:21:31 schrieb Dirk Eddelbuettel: > On 25 April 2017 at 16:11, Johannes Ranke wrote: > | This looks similar to what I got this morning when I tested my > | (unreleased) > | backport of R 3.4.0 to Debian jessie. My test was > | > | library(MASS) > | example(rlm) > | > | and there was an object that was not found. I am on a train on the way

--- nva/set_timings.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/nva/set_timings.c b/nva/set_timings.c index 7376486..985a707 100644 --- a/nva/set_timings.c +++ b/nva/set_timings.c @@ -506,7 +506,7 @@ shallow_dump(struct nvamemtiming_conf *conf) if (conf->range.start == (unsigned char) -1) conf->range.start = 0; if (conf->range.end == (unsigned char) -1) -

Hello all. I have been working on a (fairly simple) function for way too long. I’ve really hit a wall and I was hoping someone might be able to point me in the right direction. I have (attempted) to create a function that has an embedded while loop. The while loop works fine by itself, however, when the while loop is embedded in the function, the function fails. I’m not sure why this

Dear R Help, At the risk of sounding na?ve, is there a way for obtaining something resembling an r^2 value for the fit obtained using SSasymp? I have read ?SSasymp and looked through the archives but to no avail. Many thanks, Robert Terwilliger

Guys, apache cpus usage is hitting 100% sometimes ( to such an extent that its very noticeable) on a box with just 8 users or so. i m getting this when i run 'top'. The worrying thing is seeing the work 'atack' under command PID USER PR NI VIRT RES SHR S %CPU %MEM TIME+ COMMAND 23119 apache 15 0 964 556 472 S 0.7 0.0 0:03.68 atack 23479 apache 15

Hi, I'm a new R user, with a lot of questions. At the moment I'm stoped on an error traying to fit a model: > x <- sandeel ## numeric data (2500-60000) > y <- Noss ## numeric data (0-1.2) > A <- 0.8 > B <- 0.6 > C <- 1/40000 > nls( y ~ A-B*exp(-C*x)) Error in match.call(definition, call, expand.dots) : .Primitive... is not a function I'm