[R] Plotting Confidence Intervals
Steve Taylor
steve.taylor at aut.ac.nz
Sun May 3 22:36:30 CEST 2015
Have you tried:
library(effects)
plot(allEffects(ines),ylim=c(460,550))
-----Original Message-----
From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Andre Roldao
Sent: Saturday, 2 May 2015 2:50p
To: r-help at r-project.org
Subject: [R] Plotting Confidence Intervals
Hi Guys,
It's the first time i use R-Help and i hope you can help me.
How can i plot conffidence intervals? with the data bellow:
#Package Austria
library(car)
#head(States)
States1=data.frame(States)
ines=lm(SATM ~ log2(pop) + SATV , data=States1)
summary(ines)
NJ=as.data.frame(States["NJ",c(4,2,3)]) #Identificação do estado NJ
p_conf<- predict(ines,interval="confidence",NJ,level=0.95)
p_conf #Intervalo de confiança para o estado NJ e para um nivel de 95%
round(p_conf, digits=3)
p_conf1<- predict(ines,interval="confidence",NJ,level=0.99)
p_conf1 #Intervalo de confiança para o estado NJ e para um nivel de 99%
round(p_conf, digits=3)
p_pred2<- predict(ines,interval="prediction",NJ,level=0.95)
p_pred2 #Intervalo de perdição para o estado NJ e para um nivel de 95%
round(p_pred2,digits=3)
p_pred3<- predict(ines,interval="prediction",NJ,level=0.99)
p_pred3 #Intervalo de perdição para o estado NJ e para um nivel de 99%
round(p_pred3,digits=3)
Thanks
[[alternative HTML version deleted]]
______________________________________________
R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
More information about the R-help
mailing list