Baixe o app para aproveitar ainda mais
Prévia do material em texto
Rotina – aulasR aula2 <- read.table(file="C:/Users/Lorena/Documents/Arquivos do R/rh213.txt",header=T) aula2 names(aula2) #criando variaveis lnh<- log(aula2$h) lnh invd<- (1/aula2$d) invd lnhd<- log(aula2$hd) lnhd ddq<- (aula2$d/aula2$dq) ddq lnd<- log(aula2$d) lnd d2<- (aula2$d^2) d2 lninvd<- log(invd) lninvd hhd<- (aula2$h/aula2$hd) hhd dhd<- (aula2$d/aula2$hd) dhd dhd2<- (aula2$d/aula2$hd)^2 dhd2 mod1<- lm(lnh~invd+lnhd,data=aula2) mod1 names(mod1) residuals(mod1) mod1est<- exp(fitted.values(mod1)) mod1est erro1<- aula2$h-mod1est erro1 epmod1<- (erro1/aula2$h)*100 epmod1 media1<- mean(epmod1) media1 #syx variaveis somaerro²<- sum(erro1^2) somaerro² glres<- (69-3) glres syx<-sqrt(somaerro²/glres) syx syxporc<- (syx/mean(aula2$h))*100 syxporc #r² variaveis cmod1<- (sum(aula2$h))^2/69 cmod1 h2<- aula2$h^2 h2 n1<- (69-1) n1 sqtotal1<- (sum(h2))-cmod1 sqtotal1 r²<- 1-(somaerro²/sqtotal1) r² r²ajus<- 1-((1-r²)*(n1/glres)) r²ajus mod2<- lm(aula2$h~ddq) mod2 erro2<- residuals(mod2) erro2 mod2est<- fitted.values(mod2) mod2est epmod2<- (erro2/aula2$h)*100 epmod2 media2<- mean(epmod2) media2 anova(mod2) summary(mod2) #syx2 variaveis somaerro²2<- sum(erro2^2) somaerro²2 glres2<- (69-2) glres2 syx2<- sqrt(somaerro²2/glres2) syx2 syxporc2<- (syx2/mean(aula2$h))*100 syxporc2 mod3<- lm(lnh~lnd) mod3 residuals(mod3) mod3est<- exp(fitted.values(mod3)) mod3est erro3<- aula2$h-mod3est erro3 epmod3<- (erro3/aula2$h)*100 epmod3 media3<- mean(epmod3) media3 #syx3 variaveis somaerro²3<- sum(erro3^2) somaerro²3 glres3<- (69-2) glres3 syx3<- sqrt(somaerro²3/glres3) syx3 syx3porc<- (syx3/mean(aula2$h))*100 syx3porc #r² - variaveis cmod3<- (sum(aula2$h))^2/69 cmod3 h3<- aula2$h^2 h3 n3<- 69-1 n3 sqtotal3<- (sum(h3))-cmod3 sqtotal3 r²<- 1-(somaerro²3/sqtotal3) r² r²ajus<- 1-((1-r²)*(n3/glres3)) r²ajus mod4<- lm(aula2$h~aula2$d+d2+aula2$hd,data=aula2) mod4 erro4<- residuals(mod4) erro4 mod4est<- fitted.values(mod4) mod4est epmod4<- (erro4/aula2$h)*100 epmod4 media4<- mean(epmod4) media4 summary(mod4) syx4porc<- (1.772/mean(aula2$h))*100 syx4porc mod5<- lm(lnh~lninvd) mod5 residuals(mod5) mod5est<- exp(fitted.values(mod5)) mod5est erro5<- aula2$h-mod5est erro5 epmod5<- (erro5/aula2$h)*100 epmod5 media5<- mean(epmod5) media5 #syx variaveis somaerro²5<- sum(erro5^2) somaerro²5 glres<- (69-2) glres syx5<- sqrt(somaerro²5/glres) syx5 syx5porc<- (syx5/mean(aula2$h))*100 syx5porc #r² - criando variaveis cmod5<- (sum(aula2$h)^2/69) cmod5 h5<- aula2$h^2 h5 n5<- 69-1 n5 sqtotal5<- (sum(h5)-cmod5) sqtotal5 r²<- (1-(somaerro²5/sqtotal5)) r² r²ajus<- 1-((1-r²)*(n5/glres5)) r²ajus mod6<- lm(hhd~dhd+dhd2) mod6 names(mod6) residuals(mod6) mod6est<- aula2$hd*(fitted.values(mod6)) mod6est erro6<- aula2$h-mod6est erro6 epmod6<- (erro6/aula2$h)*100 epmod6 media6<- mean(epmod6) media6 #syx variaveis somaerro²6<- sum(erro6^2) somaerro²6 glres6<- (69-3) glres6 syx6<- sqrt(somaerro²6/glres) syx6 syx6porc<- (syx6/mean(aula2$h))*100 syx6porc #r² - criando variaves cmod6<- (sum(aula2$h)^2)/69 cmod6 h6<- (aula2$h^2) h6 n6<- 69-1 n6 sqtotal6<- (sum(h6)-cmod6) sqtotal6 r²<- (1-(somaerro²6/sqtotal6)) r² r²ajus<- 1-((1-r²)*(n6/glres6)) r²ajus #Gráficos: Erro/DAP par(mfrow=c(2,3)) plot(aula2$d,epmod1,xlim=c(0,20),ylim=c(-5,5),main="modelo 1",xlab="DAP",ylab="erro %") abline(h=0) plot(aula2$d,epmod2) plot(aula2$d,epmod3) plot(aula2$d,epmod4) plot(aula2$d,epmod5) plot(aula2$d,epmod6) #Graficos: Erro/H estimada plot(mod1est,erro1,xlim=c(0,30), ylim=c(-5,5),main="modelo1", xlab="H estimada", ylab="Erro (m)") abline(h=0) plot(mod2est, erro2) plot(mod3est,erro3) plot(mod4est,erro4) plot(mod5est,erro5) plot(mod6est,erro6)
Compartilhar