similar to: Confidence Bands on Regression Line

I am submitting the command: calcium<-read.table("calcium.txt",header=TRUE) and with Rw0630 and Bdr0631 I get an application error and R crashes. But the command works fine with rjune. I am running under NT4.0 Anyone have a similar problem? Jeff jmorris at ocdus.jnj.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

> Date: Mon, 11 Jan 1999 15:23:10 +0100 > From: Kjetil Halvorsen <khal at alumni.uv.es> > To: "Morris, Jeffrey (JCD-RL)" <JMorris at OCDUS.JNJ.COM> > CC: "'R-Help'" <r-help at stat.math.ethz.ch> > Subject: [R] R: Regression with t errors? > > Has anybody implemented/knows about ways to do in R regression > based on

I'm running Guido's latest version on NT 4.0 and when I try to copy a plot to the clipboard, I get an information window that says "Impossible to Open the Clipboard". Has anyone else seen this? BTW, it works fine on my laptop running Windows95. Thanks! Jeff jmorris at ocdus.jnj.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing

In the help file for the non-central t, the following appears: ncp: non-centrality parameter delta; currently `ncp <= 37.62'. I assume that this means the ncp cannot exceed 37.62. Is this still the case and is there any plans to increase this restriction? Thanks! Jeff Jeff Morris Design Support Clinical Chemistry R&D Ortho-Clinical Diagnostics email: jmorris6 at ocdus.jnj.com

I am trying to run, in R -Windows95, the example in MASS2, section 17.2. When I run predplot() on cush.lda I get > predplot(cush.lda, "LDA") Error: couldn't find function "max.col" I get the plot OK, but within predict(), it balks on max.col A more general question about using VR. Is there an easy way to use library on VR in general or do I have to move each into

I have downloaded an installed the xgobi library and can plot(laser), but when I try xgobi(laser), I get these messages in the MSDOS window... Out of environment space Out of environment space Bad command or filename and in the RGUI window it comes back with... D:\RW0901/library/xgobi/scripts/xgobi.bat -vtitle 'laser' C:/WINDOWS/TEMP/xgobi24464 Any ideas what I need to do? I am

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"

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

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

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

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

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 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

Hi, I am using ns() to model the effect of time on some outcome y [ specifically, I am using polr() in a model of the form mod1=polr(y~x1+x2*ns(Year,df=3),...) , with x1 and x2 denoting several covariates each ] I understand how to use the spline basis as recorded in the model matrix in order to reproduce the model fit and to generate curves of the point estimates of the time effect,

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,

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

Hello Fellow R Users,I have spent the last week trying to find a work around to this problem and I can't seem to solve it. I simply want to plot my GEE model result with 95% confidence bands. I am using the geepack package to run a basic GEE model involving nestling weights, to a Gaussian distribution, with "exchangeable" error structure. I am examining how nestling weight varies

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

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

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 <-