Displaying 1 result from an estimated 1 matches for "ilustracao".
2007 Jun 19
1
help w/ nonlinear regression
...),ylim=c(0,1),pch=16,cex=1.5,
xlab="CO2 (ppm)", ylab="Proporção de respostas correctas")
# Estimaçao (Método Mínimos Quadrados)
ajuste<-nls(Prop~(1/3+exp(B*(T-log(CO2))))/(1+exp(B*(T-log(CO2)))),
data=data.frame(CO2=CO2,Prop=Prop),start=list(B=-10,T=5))
summary(ajuste)
# Ilustracao do ajuste
PropEsp<-predict(ajuste,newdata=list(CO2=seq(0,1000,length=100)),se.fit=T)
lines(PropEsp~seq(0,1000,length=100),lwd=2,col=6)
# IC
upIC<-PropEsp+qt(.975,summary(ajuste)$df[2])*summary(ajuste)$sigma
loIC<-PropEsp-qt(.975,summary(ajuste)$df[2])*summary(ajuste)$sigma
lines(upIC~se...