similar to: cubic spline/smoother with nlme

Greetings, Does anyone know how to get an id number in the little header above each individual plot within a trellis plot? The default seems to be to print the word id and add a line indicating on a linear scale where the current id sits. Thanks in advance for any help you can send, Suzette ================================= Suzette Blanchard, Ph.D. Research Scientist Frontier Science

I am not able to download packages ("xtable") from CRAN not even just a .zip file, I get that the page can not be found. Please could you help? Suzette ================================= Suzette Blanchard, Ph.D. Research Scientist Frontier Science

Dear listers, I am trying to fit a model using nlsList() using alternately a SSfol() selfstart function or its developped equivalent formulae. This preliminary trial works well mydata<-groupedData(Conc~Tps|Organ,data=mydata) mymod1<-nls(Conc~SSfol(Dose,Tps,lKe,lKa,lCl),data=mydata) as well as a developped form: mymod2<-nls(Conc~Dose * exp(lKe+lKa-lCl) *

Hi All, This question was brought up some time ago but I never saw a reply so I'd like to bring it up again. When using ROracle package (version 0.5-5), I am unable to run any queries that are greater than 4000 characters in length. If I do, I get the following message: Error in oraPrepareStatement(con, statement, bind=NULL) : RS-DBI driver: (too long a statement -- it must has less than

Hello. I want to estimate the predicted values and standard errors of Y=f(t) and its first derivative at each unique value of t using the smooth.spline function. However, the data (plant growth as a function of time) show substantial heterogeneity of variance since the variance of plant mass increases over time. What is the consequence of such heterogeneity of variance in terms of bias in the

Hello. I have data (plant weights over time) that are non-linear and in which the variance increases over time. I have to estimate the first derivatives of plant weight given time (i.e. growth rate) and their se, using a regression smoother, and I have been considering cubic spline smoothers. However, I do not know if this can be done given that the error variance would increase over time.

Greetings, \Sexpr{} has worked on MikTex with earlier versions of R, I can not seem to get it to work with R-2.4.0. I run Sweave with MikTex using the following statement in Run under Accessories. latex -include-directory="C:\Program Files\R\R-2.4.0\share\texmf" "C:\Documents and Settings\Suzette\Desktop\MyFile\MyProgram.tex" Any advice would be much appreciated, Suzette

Dear friends of nlme, Running quinModel (Pinheiro/Bates page 380) on R (current release, windows) gives: Nonlinear mixed-effects model fit by maximum likelihood Model: conc ~ quinModel(Subject, time, conc, dose, interval, lV, lKa, lCl) Data: Quinidine Log-likelihood: -497 Fixed: lV + lKa + lCl ~ 1 lV lKa lCl 5.382 -0.273 2.470 Random effects: Formula: list(lV ~ 1, lCl ~ 1)

Greetings, I implemented BRugs to run the EWOC model with a cohort size N=1. I output the simulation data using bugsdata(data), where data is the following list. > data $Dose [1] 140 $Y [1] 0 bugsdata(data) puts out the file data.txt as follows. list(Dose=1.40000E+02, Y=0.00000E+00) after I type the following line modelData("data.txt")

Greetings, I have a simulation of a nonlinear model that is failing. But it does not fail til way into the simulation. I would like to look at the run that is failing and maybe I could if I could capture the seed for the failing run. The help file on set.seed says you can do it but when I tried rs<-.Random.seed print(paste("rs",rs,sep=" ")) I got 626 of them so I

I would like to add two horizontal lines representing acceptible drug levels to a trellis plot. I tried using abline and I get an error that "plot.new has not been called." See below. xyplot(FK~WEEK|Event1/MRN, data=FKdat.o1) abline(h=5) abline(h=10) Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : plot.new has not been called yet Any help that can be provided

Greetings, I used latex type code in my Rd files. The pdf version using R CMD Rd2dvi --output=PKtools.pdf --pdf --title="PKtools" PKtools man came out quite nice. However, my current HTML version does not tex the latex so there is latex code in the files and looks bad. example problem code: AIC$_c$ or \begin{itemize} \item NLME: \begin{itemize} \item population level:

This concerns the "Clinical Study of Quinidine" example on page 380 of the book "Mixed-Effects Models in S and S-PLUS" by Pinheiro and Bates (2000). I have tried to reproduce the example, but get an error: > library(nlme) > fm1Quin.nlme <- nlme(conc ~ quinModel(Subject, time, conc, dose, interval, lV, lKa, lCl), + data=Quinidine, +

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

Dear list, I do have a problem with nls. I use the following data: >test time conc dose 0.50 5.40 1 0.75 11.10 1 1.00 8.40 1 1.25 13.80 1 1.50 15.50 1 1.75 18.00 1 2.00 17.00 1 2.50 13.90 1 3.00 11.20 1 3.50 9.90 1 4.00 4.70 1 5.00 5.00 1 6.00 1.90 1 7.00 1.90 1 9.00 1.10 1 12.00 0.95 1 14.00

Hello everyone,          Dr. wood told me that I can adapting his example to force cubic spline to pass through certain point.          I still have no idea how to achieve this. Suppose we want to force the cubic spline to pass (1,1), how can I achieve this by adapting the following code? # Penalized example: monotonic penalized regression spline ..... # Generate data from a monotonic truth.

Hallo, I'm facing a problem and I would really appreciate your support. I have to translate some Matalb code in R that I don't know very well but I would like to. I have to interpolate 5 point with a cubic spline function and then I expect my function returns the Y value as output a specific X value inside the evaluation range. Let's suppose that: 1- *X = [-10, -5, 0, 5, 10]* 2

Dear all, Is the restricted cubic spline function working properly in the glm model? We used glm(y~rcs(x,5), family=binomial) but it seems that for some theoretical reasons the rcs, restricted cubic spline function can not be fitted by a glm function. Is this correct? Regards, Jan ((Originally, we used lrm(y~ rcs(x,5)) but we couldn't find how to derive the AIC value of the fitted model.

Hi, all I am trying to get the basis matrix and penalty matrix for natural cubic splines. In the "splines" package of R,"ns" can generate the B-spline basis matrix for a natural cubic spline. How can I get the basis matrix and penalty matrix for natural cubic spline. Thanks a lot! Lee [[alternative HTML version deleted]]

Hello everone,            Anyone who knows how to force a cubic smoothing spline to pass through a particular point?            I found on website  someone said that we can use "cobs package" to force the spline pass through certain points or impose shape           constraints (increasing, decreasing). However,  this package is using  B-spline and can only do linear and quadratic