Displaying 20 results from an estimated 20000 matches similar to: "Plot Confidence Interval for Regression"
2008 Jan 06
1
Confidence Inteval for regression curve.
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]]
2006 May 30
1
Piecewise (broken stick) models in R
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
2006 Sep 27
1
matrix with additional upper, botton, left and right cells
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
2011 Jan 11
1
Confidence interval on quantile regression predictions
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)
2006 Dec 10
4
sample "n" random positions from a matrix
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]]
2006 Dec 06
2
require(simecol) error
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
2011 Oct 30
1
Confidence interval for Support Vector Regression
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
2009 Mar 08
0
statistical question: confidence interval of regression weight - significance
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
2011 May 21
1
Confidence Interval of Specral Density Plot
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
2006 Jun 17
2
substr or split help needed
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
2008 Nov 29
1
function for simultaneous confidence interval of regression coefficients
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
2013 Feb 07
2
how to draw confidence interval lines of a fitted curve of polynominal regression
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]]
2010 Nov 30
1
confidence interval for logistic joinpoint regression from package ljr
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
2013 Apr 01
1
95% Confidence Interval for a p-p plot
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
2008 Jan 04
1
redirect plot device automatically
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
2008 Jan 11
2
clipping a large image on R
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
2010 Mar 14
1
confidence intervals for non-linear regression
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
2011 Feb 16
2
Confidence interval of Survival Curve of Weighted Cox Regression
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
2004 Mar 12
1
confidence interval in local polynomial regression
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
2008 May 19
0
How to get confidence interval and coefficient in Logic Regression
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