similar to: predictive interval in nlme

How can I get prediction intervals from a mixed-effects model? Consider the following example: library(nlme) fm3 <- lme(distance ~ age*Sex, data = Orthodont, random = ~ 1) df3.1 <- with(Orthodont, data.frame(age=seq(5, 20, 5), Subject=rep(Subject[1], 4), Sex=rep(Sex[1], 4))) predict(fm3, df3.1, interval='prediction') # M01 M01

Hello- It has been several years since anyone has asked, so i am asking again- has anyone created a routine to estimate confidence intervals for predictions from nls models (ala Bates and Watts 1988)? Thanks - Alice Shelly [[alternative HTML version deleted]]

something like "sprintf" in C? so I can do: print(sprintf("the correct result is %3.4f\n", myresult)); ------- Also, I am desperately looking for a "clear console screen" function in R... thanks a lot! [[alternative HTML version deleted]]

Is there a function in R for doing Greenhouse-Geiser correction in ANOVA models? Is it already available in the aov function? How do we use it? Best, Darren [[alternative HTML version deleted]]

Hello. Some months age an e-mail was posted in which a comparison between Glimmix and glm was discussed. I have not been able to find that e-mail on the R archive. Does anyone recall the date of the above e-mail? Thank you very much. ******************************************* Antonio Paredes USDA- Center for Veterinary Biologics Biometrics Unit 510 South 17th Street, Suite 104 Ames, IA 50010

Hi all I have a question regarding differences in the way gregpexr works in R 2.3.0 and R 2.4.0. In R 2.3.0, this is what happens: > gregexpr(" [a-z] [a-z] ", " a b c d e f ", perl=T) [[1]] [1] 1 3 5 7 9 attr(,"match.length") [1] 5 5 5 5 5 ... while in R 2.4.0, this is what happens: > gregexpr(" [a-z] [a-z] ", " a b c d e f ", perl=T)

Hello: Could anyone point me to a nice example where someone has created methods for "[" on a user defined Class? I looked at the package Matrix but that was a little daunting. I'm looking for someone a little more introductory. I've tried to search the help section and the web but its difficult since "[" isn't searchable. Thanks in advance! Greg [[alternative

Inspired by the functions that Barry Rawlingson and Dave Forrest posted for searching Rwiki and R-help archive, I've made up a function that does the search on Prof. Baron's site (Thanks to Prof. Baron's help on setting up the query string!): RSiteSearch <- function(string, restrict="Rhelp", format="long", sortby="score",

Based on some recent r-help discussions, I have been trying out plotting confidence intervals using predict and matplot. Matplot appeared to not be plotting the linear regression when using the default column names generated by read.table (V1, V2, etc). On further inspection, the error seemed to be with predict and vector names (V1, V2) rather than with matplot. I was using some textbook

Hi, I´m trying to plot a nonlinear regresion with the confidence bands for the curve obtained, similar to what nlintool or nlpredci functions in Matlab does, but I no figure how to. In nls the option is there but not implemented yet. Is there a plan to implement the in a relative near future? Thanks in advance, Florencio La información contenida en este e-mail y sus ficheros adjuntos es

Hello everybody, I have two series of intervals, and I'd like to output the shared regions. For example: series1<-cbind(Start=c(10,21,40,300),End=c(20,26,70,350)) series2<-cbind(Start=c(25,60,210,500),End=c(40,100,400,1000)) > series1 Start End [1,] 10 20 [2,] 21 26 [3,] 40 70 [4,] 300 350 > series2 Start End [1,] 25 40 [2,] 60 100 [3,]

Hello useRs and helpRs, Some time ago, in a gallaxy far away (here is the thread : http://finzi.psych.upenn.edu/R/Rhelp02a/archive/46532.html ) we discussed about a graph gallery showing the power of R in that domain. I did some work around that, and there is a (pretty advanced) draft here : http://addictedtor.free.fr/graphiques/displayGallery.php For instance, there are some of my graphs,

This below is not solvable with uniroot to find "a": fn=function(a){ b=(0.7/a)-a (1/(a+b+1))-0.0025 } uniroot(fn,c(-500,500)) gives "Error in uniroot(fn, c(-500, 500)) : f() values at end points not of opposite sign" I read R-help posts and someone wrote a function: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/92407.html but it is not very precise. Is there any

hie all i am trying to carry out a categorical data analysis but my problem is that when in i use the chi squared test some of my expected values are less than 5. is there a test that can handle this situation. the data is not a 2*2 table. its more from the social sciences where you have from strongly agree to strongly disagree. i know i can collapse vthe tables but there is a loss of

Hi all, Is it possoble to include an R function in a Perl program? I would like to use the "limma" library for microarray analysis, and have no clue how this can be done. Thanks in Advance -Tarun

The CRAN host in Tampa, FL (cran.hostingzero.net) isn't responding, and hasn't responded in quite a while -- at least problems were reported more than a year ago (Oct 2007) http://finzi.psych.upenn.edu/R/Rhelp02a/archive/116706.html although someone apparently succeeded in July 2008 http://finzi.psych.upenn.edu/R/Rhelp02a/archive/137451.html Perhaps it's just flaky, and not

Hello, I have a deriv function that I am feeding to nlme. It works, and I can use it in nls, but when I try to use it in nlme I get Error: subscript out of bounds. I can fit the model using SSasympOrig, instead of the deriv function, but I am trying to reproduce an earlier analysis (done in Splus) and I get slightly different results with SSasympOrig. Here are my calls: This does fit

Dear Sirs, I would like to kindly ask for your assistance on the folllowing issue. After uploading the most updated the versions of some libraries (namely: survival and nlme) I became unable to open R. Every time I try to open R a dialogue box pops up with the following message: --------------------------- Information --------------------------- Fatal error: unable to restore saved data in

Hi to all, I found in the R-help archive how to calculate the p-value for a gee result: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/74150.html but there are two questions (I am afraid they are basic questions ...) 1. why is the result multiplicated with 2 2. how could I decide between lower.tail =TRUE and FALSE: example:

Hi, I' trying to learn how to use lme for Linear Mixed Models and I have a problem when I have to include more than one random effect in my model. I know that this could be a stupid question to ask, but I'm not able to solve it by myself... One example: if my model is response = operator + block + day with operator and block as fixed effects and day as random effect, I use res.lme