Regressão Exponencial
Boa Tarde! Tenho os seguintes dados: x = c(1,1,2,2,1,3,10,15,9,17,13,18,21,27) tempo=c(2001,2003,2004,2006,2007,2008,2009,2010,2011,2012,2013,2014,2015,2016) Ao plotar (x,y), verifico um crescimento exponencial do x ao longo do tempo. Alguém teria o script de regressão exponencial para me ajudar? Grata Fabiana
Olá! O ajuste fica bom, mas o parâmetro "a" tem problemas (alguém da lista mais especializado poderia falar algo sobre). Att. y = c(1,1,2,2,1,3,10,15,9,17,13,18,21,27) tempo=c(2001,2003,2004,2006,2007,2008,2009,2010,2011,2012,2013,2014,2015,2016) library(minpack.lm) exp<-nlsLM(y~a*exp(b*tempo),start = c(a=0,b=0.25),control=list(maxiter=500)) summary(exp) a<-coef(exp)[1] b<-coef(exp)[2] par(mar=c(4,4,1,1)) plot(tempo,y,las=1) curve(a*exp(b*tempo),xname="tempo",add = T) ________________________________ De: R-br <r-br-bounces@listas.c3sl.ufpr.br> em nome de Fabiana Goncalves Barbosa via R-br <r-br@listas.c3sl.ufpr.br> Enviado: domingo, 17 de setembro de 2017 13:53 Para: r-br@listas.c3sl.ufpr.br Assunto: [R-br] Regressão Exponencial Boa Tarde! Tenho os seguintes dados: x = c(1,1,2,2,1,3,10,15,9,17,13,18,21,27) tempo=c(2001,2003,2004,2006,2007,2008,2009,2010,2011,2012,2013,2014,2015,2016) Ao plotar (x,y), verifico um crescimento exponencial do x ao longo do tempo. Alguém teria o script de regressão exponencial para me ajudar? Grata Fabiana
Fabiana, Se você não se importar de fazer um passo da matemática dessa solução "na mão", a regressão exponencial é normalmente feita sem se necessitar de pacotes mais sofisticados (nada contra o uso deles) empregando a transformação da variável pra tornar regressão num caso de regressão linear:
regr <- lm(log(x) ~ tempo) summary(regr) summary(lm(log(x) ~ tempo))
Call: lm(formula = log(x) ~ tempo) Residuals: Min 1Q Median 3Q Max -1.18204 -0.17017 -0.08422 0.30765 0.77378 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -502.05994 59.90223 -8.381 2.33e-06 *** tempo 0.25074 0.02981 8.410 2.24e-06 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 0.4974 on 12 degrees of freedom Multiple R-squared: 0.855. Adjusted R-squared: 0.8429 F-statistic: 70.73 on 1 and 12 DF. p-value: 2.245e-06 HTH -- Cesar Rabak 2017-09-17 13:53 GMT-03:00 Fabiana Goncalves Barbosa via R-br < r-br@listas.c3sl.ufpr.br>:
Boa Tarde!
Tenho os seguintes dados:
x = c(1,1,2,2,1,3,10,15,9,17,13,18,21,27) tempo=c(2001,2003,2004,2006,2007,2008,2009,2010,2011,2012, 2013,2014,2015,2016)
Ao plotar (x,y), verifico um crescimento exponencial do x ao longo do tempo.
Alguém teria o script de regressão exponencial para me ajudar?
Grata Fabiana
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
participantes (3)
-
Cesar Rabak -
Delcio Rudinei Bortolanza -
Fabiana Goncalves Barbosa