similar to: Confidence bands for model estimates using ns() spline basis

Hi all, I'm looking for smooth monotone estimation packages, preferably using splines. I downloaded the 'cobs' package and intend to use it, but since it offers only quadratic splines based on L1 minimization, I'd like to compare its performance to that of a more 'mainstream' cubic-spline, L2-norm minimizing spline. Preferably a smoothing spline. Does anyone know of such

Hello, I am running a model with four categories and want predicted probabilities in each category. Now for this example I wont give a counterfactual just the training data is fine but is there anyway to get a confidence interval around the predicted probabilities in each group? I have tried but it gives me probabilities and I have used interval="confidence", level=.095 and then interval

Hi all, I want to create a rather standard overlaid qqnorm plot on a single variable, with different subgroups of the same dataset plotted using different colors/symbols/etc. (I don't want side-by-side, rather different-colored curves on the same graph) I managed to do it rather tediously using "lapply" and "split", and wondered whether there is any single-command

Hi all, Could anybody recommend some easy-to-understand and example based notes/tutorials on how to use cubic splines to do filtering on multi-dimension data? I am confused by the 1-dimensional case, and more confused by multi-dimensional case. I found all the books suddenly become very abstract when it comes to this subject. They don't provide examples in R or Splus at all. Specifically,

I'm using nonparametric regression (packeges ksmooth and ks). My question: is there any way to compute confidence bands (or variability bands) with R. Confidence bands for functions are intervals [CLO(x);CUP(x)] such that with probability 1-alpha the true curve is covered by the band [CLO(x);CUP(x)]. Thanks very much for any help you can offer. Michael G??lger

What's the best / simplest way to create 95% confidence bands for a model created with lm() that can be plotted around teh regression line? I've looked everywhere for this - I guess I must be missing something. - Jason

Is there a built-in function in R that will generate simultaneous confidence and prediction bands for linear regression? Tom -- View this message in context: http://www.nabble.com/Simultaneous-Confidence-Prediction-Bands-tp17941537p17941537.html Sent from the R help mailing list archive at Nabble.com.

Dear all, I am trying to plot KM curves with confidence bands with function survplot under package rms. However, the following codes do not seem to work. The KM curves are produced, but the confidence bands are not there. Any insights? Thanks in advance. library(rms) ########data generation############ n <- 1000 set.seed(731) age <- 50 + 12*rnorm(n) label(age) <- "Age"

Dear all, Is there a way to plot / calculate pointwise confidence bands or interval values for a non parametric regression like sm.regression? Thank you in advance. Regards, Martin

Dear R-list, I would like to show the implications of estimating a linear trend to time series, which contain significant serial correlation. I want to demonstrate this, comparing lm() and an gls() fits, using the LakeHuron data set, available in R. Now in my particular case I would like to draw confidence bands on the plot and show that there are differences. Unfortunately, I do not know how to

Hello R users, Given a linear (in the parameters) regression model where one predictor x interacts with time and time*time (ie, a quadratic effect of time t): y = b0 + b1(x) + b2(t) + b3(t^2) + b4(x*t) + b5(x*t^2) + e, I would like to construct 95% confidence bands (optimally, shaded) around this function: *dy* = b1 + b4(t) + b5(t^2) *dx* That is, the partial effect of x on y changing over

Hello, I am trying to plot confidence bands on the mean and prediction bands for the following nonlinear regression, using maximum likelihood via optim. A toy example with data and code of what I am trying to accomplish is: VOL<-c(0.01591475, 1.19147935 ,6.34102460, 53.68809287, 91.90143074, 116.21397007, 146.41843056, 215.64535337, 256.53149673, 315.73609232) Age <-c(1.622222, 2.833333

I adjusted an exponential regression to the following data and wish to plot confidence bands as well. Is this possible? Any help greatly appreciated. Claudia x <- c(1989,1990,1991,1992,1993,1994,1995,1996,1997,1998,1999,2000,2001, 2002,2003,2004,2005,2006,2007,2008,2009) y <- c(987,937,810,749,1087,807,1050,1294,1455,1022,927,403,698,1191,1078, 1176,1125,936,1263,647,868) dat <-

Is there a easy way to put confidence bands on a regression plot or do I have to explictly construct the values and send the lines to the plot? Thanks! Jeff jmorris at ocdus.jnj.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

Dear R experts: I need to get this plot, but also with 95% confidence interval bands: hour <- c(1, 2, 3, 4, 5, 6) millivolts <- c(3.5, 5, 7.5, 13, 40, 58) plot(hour, millivolts, xlim=c(1,10), ylim=c(0,1000)) pm <- lm(millivolts ~ poly(hour, 3)) curve(predict(pm, data.frame(hour=x)), add=TRUE) How can the 95% confidence interval band curves be plotted too? Sincerely,

To my understanding, a confidence interval typically covers a single valued parameter. In contrast, a confidence band covers an entire line with a band. In regression, it is quite common to construct confidence and prediction bands. I have found that many people are connecting individual confidence/prediction interval values produced with predict(object,sd.fit=T,type="conf/pred") and

I have two questions related to plotting predicted values for a linear mixed model using xyplot: 1: With a groups= argument, I can't seem to get the key to appear inside the xyplot. (I have the Lattice book, but don't find an example that actually does this.) 2: With lme(), how can I generate confidence bands or prediction intervals around the fitted values? Once I get them, I'd

Dear R experts, beginners and everyone else, I'm calculating "cumulative periodogram" using the command "cpgram" [1] from the MASS library. Here is a short example with the "lh" (hormone level) dataset: library(MASS) plot(lh,type="l",ylab="value",xlab="time", main="Hormone Levels (lh)") spectrum(lh,

Dear list, I'm a bit perplexed why the 95% confidence bands for the predicted probabilities for units where x=0 and x=1 overlap in the following instance. I've simulated binary data to which I've then fitted a simple logistic regression model, with one covariate, and the coefficient on x is statistically significant at the 0.05 level. I've then used two different methods to

Hi, I have the following data: > est sch190 sch107 sch290 sch256 sch287 sch130 sch139 4.16656026 2.64306071 4.22579866 6.12024789 4.49624748 11.12799127 1.17353917 sch140 sch282 sch161 sch193 sch156 sch288 sch352 3.48197696 -0.29659410 -1.99194986 10.23489859 7.77342138 6.77624539 9.66795001 sch368