Análise de Dados de Árvores

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)