[R] reading the seed from a simulation
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 don't know how to identify the right one. Please can you help? Thank you, Suzette = Suzette Blanchard, Ph.D. UCSD-PPRU __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] cubic spline/smoother with nlme
Greetings, I would like to use a cubic spline or smoother to model the fixed effects within nlme. So far the only smoother I have been able to get to run successfully in nlme is smooth(). I tried smooth.spline: fixed=list(lKa~1,lCL~smooth.spline(BSA, df=3)) the error I got was the following. Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : invalid variable type Can anyone suggest a cubic spline that would work within this context? Thank you for any help you can send, Suzette = Suzette Blanchard, Ph.D. UCSD-PPRU __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] nlme
Greetings, I am trying to use nlme to model a large data set of pharmacokinetic concentrations with sparse data per cluster. The base model is a one compartment model with first order absorption. The program runs but the PNLS step never converges. The parameter estimates are reasonable. The plot of the individual level residuals is pretty reasonable and the plot of the residuals from the population model is ok but not great. My questions are: 1) I know the PNLS step starts and updates the fixed parameters and random effects. Step 2 is the LME step and it updates the interpatient variance. Does the LME step include both the EM algorithm and a Newton- Raphson ms step? 2) I set Maxiter=500 interations and I get convergence on what I think is the LME step but PNLS step never converges. Any suggestions on what I might do to achieve convergence would be very much appreciated. Thank you, Suzette = Suzette Blanchard, Ph.D. Research Scientist Frontier Science Foundation 1244 Boylston St. Suite 303 Chestnut Hill, MA 02467 Email: [EMAIL PROTECTED] Phone: (617) 632-2007 Fax:(617) 632-2001 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] trellis plot question
Yes, that solved the problem. Thank you, Suzette On Fri, 25 Jul 2003, Petr Pikal wrote: Hi You probably use numeric id variable. I think you need to make an id variable to become a factor. e.g. use as.factor(id) Cheers Petr On 24 Jul 2003 at 10:39, Suzette Blanchard wrote: 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 Foundation 1244 Boylston St. Suite 303 Chestnut Hill, MA 02467 Email: [EMAIL PROTECTED] Phone: (617) 632-2007 Fax:(617) 632-2001 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help Petr Pikal [EMAIL PROTECTED] [EMAIL PROTECTED] = Suzette Blanchard, Ph.D. Research Scientist Frontier Science Foundation 1244 Boylston St. Suite 303 Chestnut Hill, MA 02467 Email: [EMAIL PROTECTED] Phone: (617) 632-2007 Fax:(617) 632-2001 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] trellis plot question
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 Foundation 1244 Boylston St. Suite 303 Chestnut Hill, MA 02467 Email: [EMAIL PROTECTED] Phone: (617) 632-2007 Fax:(617) 632-2001 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help