similar to: Need help in nlmeODE

Dear R-helpers, I have some problem with a model I am trying to run with nlmeODE. I have a covariate that I want to include in the model and I can't find in the documentation how to include it in the model so that it can calculate it. >formula(trial) PME ~ time | Ident abcnoniso <- list(DiffEq=list( debdt = ~ -ksol*pH*eb , desdt = ~

Dear all, I have two questions: 1. Does anyone know if nlmeODE package will work once odesolve is discontinued? If i understand well nlmeODE depends on odesolve. 2. Can anyone suggest another package I could use to fit complex PKPD models (described using differential equations) to data? Thanks a lot! Magdalena

Dear All A new non-linear mixed-effects modelling package for R called nlmeODE is now available at www.cran.r-project.org and http://nlmeODE.sourceforge.net. It combines Pinheiro and Bates' NLME library in R with an ordinary differential equations (ODE) solver thereby enabling the use of differential equations in population PK/PD modelling. It includes several PK examples such as one- and

Dear list, I have a degradation model: dX/dt = * I(X2)*( k1*X(t) )/( X(t)+k2 ) where X(t) is concentration at time t, and k1 and k2 are parameters that I want to estimate. I(X) is a known inhibitor function. My questions is whether this is implemented or easily computed in any R package. I have searched the archives but without luck. Any help or comments on this would be appreciated, Klaus

Dear Colleagues, Our group is also working on implementing the use of R for pharmacokinetic compartmental analysis. Perhaps I have missed something, but > fit <- nls(noisy ~ lsoda(xstart, time, one.compartment.model, c(K1=0.5, k2=0.5)), + data=C1.lsoda, + start=list(K1=0.3, k2=0.7), + trace=T + ) Error in eval(as.name(varName), data) : Object

Hello! I would like to know if R can solve Differential Equations... I don't think so because, in my point, I see R like a Statistical System, not a Math System. Am I wrong? Thank you very much. M??rcio de Medeiros Ribeiro Graduando em Ci??ncia da Computa????o Departamento de Tecnologia da Informa????o - TCI Universidade Federal de Alagoas - UFAL Macei?? - Alagoas - Brasil Projeto CoCADa

Hi Great. I did not think that such combination is posssible. Thanks. Petr > -----Original Message----- > From: peter dalgaard [mailto:pdalgd at gmail.com] > Sent: Tuesday, July 11, 2017 1:11 AM > To: PIKAL Petr <petr.pikal at precheza.cz> > Cc: r-help at r-project.org > Subject: Re: [R] fit lognorm to cdf data > > How about > > proc <- c(0.9, 0.84, 0.5,

Hi > -----Original Message----- > From: Robert Baer [mailto:rbaer at atsu.edu] > Sent: Wednesday, August 16, 2017 3:04 PM > To: PIKAL Petr <petr.pikal at precheza.cz>; Duncan Murdoch > <murdoch.duncan at gmail.com> > Cc: r-devel at r-project.org > Subject: Re: [Rd] strange behaviour read.table and clipboard > > You said, "put a name in the cell".

Hi I copied your values to R, here it is > dput(temp) temp <- structure(list(par1 = structure(1:4, .Label = c("x1", "x2", "x3", "x4"), class = "factor"), y1 = c(-0.19, 0.45, -0.09, -0.16), y2 = c(0.4, -0.75, 0.14, -0.01), y3 = c(-0.06, -8.67, 1.42, 2.21), y4 = c(0.13, -0.46, 0.06, 0.06)), .Names = c("par1", "y1",

Hi Peter; The R code you provided works very well. Once again thanks so much for this. The number of variables in my data set that should appear on the y-axis is 733 and they are not numerical (for example the name of one variable is *palmitoyl-arachidonoyl-glycerol (16:0/20:4) [1]**. So, the plot looks very messy in one page. How can I make the plot to print out on multiple pages? Regards,

Hello, I am developing an agent-based model to simulate the spread of infectious diseases in heterogeneous landscapes composed of habitat polygons (or clumps of connected cells). To simplify the model, I consider a habitat grid (or raster) containing the polygon ID of each cell. In addition, I have epidemiological parameters associated with each polygon ID. At each time step, the parameter values

Hi Maybe there are other ways but I would split data to several chunks e.g. in list and use for cycle to fill multipage pdf. With the toy data something like library(reshape2) library(ggplot2) temp <- melt(temp) temp.s<-split(temp, cut(1:nrow(temp), 2)) pdf("temp.pdf") for (i in 1: length(temp.s)) { p <- ggplot(temp.s[[i]], aes(x=par1, y=variable, size=abs(value),

Hi Keep your messages coppied to R helplist, others could give you answers too. See in line From: Rosa Oliveira [mailto:rosita21 at gmail.com] Sent: Tuesday, August 1, 2017 4:38 PM To: PIKAL Petr <petr.pikal at precheza.cz> Subject: Re: [R] Superscript and subscrib R for legend x-axis and y-axis and colour different subjects in longitudinal data with different colours Hi Pikal, I looked

Hi Petr; Thanks so much. Exactly this is what I need. I will play to change color and so on but this backbound is perfect to me. I do appreciate your help and support. Regards, Greg On Thu, Feb 8, 2018 at 1:29 PM, PIKAL Petr <petr.pikal at precheza.cz> wrote: > Hi > > I copied your values to R, here it is > > > > > dput(temp) > > > > temp <-

Hi When you load external file to R, character variables are converted to factors by default and alphabetically sorted. I have limited connection to internet, so I cannot find the answer, you could try it yourself. Maybe you could try not to convert vector with names to factor, which, for plotting issue is not different from factor coding. See ?read.table for details However I am not sure if it

I did the same transition, facing the same issues, so I settled for installation of RKWard. This pretty much reproduces my previous experience with Tinn-R and the Rgui. HTH Ruben ________________________________ Ruben H. Roa-Ureta, Ph. D. Senior Scientist, Center for Environment and Water, Marine Studies Section, King Fahd University of Petroleum and Minerals, KFUPM Box 1927, Dhahran 31261,

How about proc <- c(0.9, 0.84, 0.5, 0.16, 0.1) size <- c(0.144, 0.172, 0.272, 0.481, 0.583) plot(size, proc, xlim=c(0,1), ylim=c(0,1)) fit<-nls(proc~plnorm(size, log(xmid), sdlog, lower=FALSE), start=list(xmid=0.2, sdlog=.1)) summary(fit) lines(fitted(fit)~size) -pd > On 10 Jul 2017, at 16:27 , PIKAL Petr <petr.pikal at precheza.cz> wrote: > > Dear all > > I am

Q1. No. Q2. What do you do? Did you follow the instructions at [1]? If you go to the top level of the R help system you should be able to open the manuals. ?help Q3. Not that I know of. If you have difficulty with R on Linux then you should focus on learning how to use Linux in general... and learn what X-Windows is. But just as with Windows, googling your problem will usually lead to stepwise

Hi try str(varcompobject) to see structure of this object. You can extract parts by standard R means. Cheers Petr > -----Original Message----- > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Sharada > Ramadass > Sent: Tuesday, August 8, 2017 3:33 PM > To: r-help at r-project.org > Subject: [R] how to extract individual values from varcomp? > >

Hi Example, example, example - preferably working. Wild guess - did you try ggplot? Cheers Petr > -----Original Message----- > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of greg holly > Sent: Thursday, February 8, 2018 8:14 AM > To: r-help mailing list <r-help at r-project.org> > Subject: [R] plotting the regression coefficients > > Hi Dear all;