similar to: Plot Confidence Interval for Regression

Dear all, I have a dataset which I need to estimate the regression model and plot the estimated curve two other curves with low and high confidence interval (CI=95%). How can I do that? x<-1:100 y<-x^0.2+rnorm(100,0.1,0.1) mod<-glm(y~log(x)) plot(y~x) lines(predict(mod)~x,col=2) Kind regards, miltinho Brazil para armazenamento! [[alternative HTML version deleted]]

Hi there, I´m a newbie in R and I´m looking for some advice how to teste Piecewise (Broken Stick) models. I´have another simple question: how can I compute descriptive statistics for data grouped for two or more variables? See below, please. Year Site Repetition rainfall 1980 1 1 ... 1980 1 2 ... 1980 1 3 ... 1980 2

Dear R Gurus, I have a matrix dim(1000x1000) and I need create a second matrix with dim(1002x1002) and insert my first matrix at position col=2,line=2. Please, see an example below: 0050055050 5550000000 5000505005 5005000500 0055550555 and I need 333333333333 300500550503 355500000003 350005050053 350050005003 300555505553 333333333333 Thanks a

I am using the quantreg package to build a quantile regression model and wish to generate confidence intervals for the fitted values. After fitting the model, I have tried running predict() and predict.rq(), but in each case I obtain a vector of the fitted values only. For example: library(quantreg) y<-rnorm(50,10,2) x<-seq(1,50,1)

Hi there, I have a binary matrix (dim 100x100) filled with values 0 and 1. I need select a record "n" positions of that matrix when values are 1. How can I do that? Thanks for all, Miltinho Brazil --------------------------------- [[alternative HTML version deleted]]

Hi there, I´m trying to use "simecol" package but I got the error showed below. I´m runnig R version 2.4.0 (2006-10-03). Kind regards, miltinho Brazil --- > require(simecol) Loading required package: simecol Error in loadNamespace(package, c(which.lib.loc, lib.loc), keep.source = keep.source) : in 'simecol' methods specified for export, but

Hi, Is it possible to calculate confidence intervals for support vector regression? In the "ksvm{kernlab}" manual, it says that it supports confidence intervals for regression. However, i couldn't find any information about how to calculate confidence interval. Do you have any documents or examples about this or any other suggestions? Thanks in advance -- [[alternative HTML

hi, at first; thanks for the help on getting confidence intervals in R. now I have a pure statistical question. I hope you don't mind if I ask ... I have an expectation of how large my beta-weight in a regression should be - so I have an "ideal" or expected regression line. Now the real beta-weight is less then the expected and when I draw the confidence interval lines

Dear Users, I wish to know at what confidence level is the confidence interval provided in the Spectrum function (plot.spec) plots. The only information provided in the help regarding this is : a confidence interval will be plotted by plot.spec: this is asymmetric, and the width of the centre mark indicates the equivalent bandwidth. Attached is an example plot with the BLUE confidence

Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060617/d28b0fb4/attachment.pl

List, Would someone be so kind as to point me to a function that will calculate simultaneous confidence intervals of regression coefficients based upon their distribution as (under the assumption of normal errors, with \mathbf{X} as the design matrix): $\hat{\mathbf{\beta}} \sim N(\mathbf{\beta}, \sigma^2(\mathbf{X}^T\mathbf{X})^{-1})$. 'confint' calculates individual coefficients so

Hello, I drew a plot of weight and height of people and fitted it with a polynominal regression x^2. (using curve()) Now I would like to draw the confidence interval line for the fitted curve. Please kindly advise the code for the purpose. Thank you. Elaine [[alternative HTML version deleted]]

I?m trying to run a logistic joinpoint regression utilising the ljr package. I?ve been using the forward selection technique to get the number of knots for the analysis, but I?m uncertain as to my results and the interpretation. The documentation is rather brief ( in the package and the stats in medicine article is quite technical) and without any good examples. At the moment I?m thinking 1)find

Hi, I want to create upper and lower 95% confidence intervals for a p-p plot of an empirical distribution with a theoretical gamma distribution. This is my code: x<-rgamma(100,shape=2, rate=1) # empirical data fitdistr(x,"gamma") # fit a gamma distribution dist<-pgamma(x,shape=1.9884256 ,rate=0.8765314 ) # fitted distribution, using the loglikelihood estimated parameters

Dear all, Every time that I need to generate I plot output as figure I use something like png("myoutput.png") plot(rnorm(100)*1.0,col=4) points(rnorm(100)*0.5,col=2,pch=0) dev.off() But now I need to generate a lot of outputs and it could not be done in a for looping where I can change the names. So I would like to know if is there a way of I automatically redirect the output from my

Dear all, I have a so large image (43,000 x 18,000 pixels) and I need clip this image with a smallest one (1000x1000 pixels). I can read the second image using rgdal package. But the first image can´t be read on my system because if memory limitation (I have about 2GB availabe). So I would like hear from you if anyone have some suggestin in this regards. I have also ArcGis 8.6 ans Erdas 8.3

Dear all, I am interested to calculate confidence interval for fitted values in general for non-linear regressions. Lets say we have y=f(x1,x2,..xN) where f() is a non-linear regression. I would like to calculate a confidence interval for new prediction f(a1,..,aN). I am aware of techniques for calculating confidence intervals for coeffiecients in specific non-linear regressions and with them

Hi, May I know how to obtain the confidence interval of the survival curve of weighted Cox regression model? I tried coxph, cph, and coxphw, but they did not work. Any help would be much appreciated. Thanks, Jeanne -- View this message in context: http://r.789695.n4.nabble.com/Confidence-interval-of-Survival-Curve-of-Weighted-Cox-Regression-tp3308085p3308085.html Sent from the R help mailing

Dear all, Is there any package or function can do the pointwise confidence interval and confidence band for the local polynomial regression? Maybe the local linear regression is enough. Thanks. Regards, Zhen

sorry to bother everyone. i have question to get the coefficient and confidence interval in Logic Regression with Logistic model. below i list the R code X <- matrix(as.numeric(runif(400) < 0.5), 50,8) colnames(X) <- paste("X", 1:ncol(X), sep="") rownames(X) <- paste("case", 1:nrow(X), sep="") # Define expected result: Y = (NOT X1) AND X5 Y