Script_Demaerschalk_(1971)

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

#1) Importando arquivo de dados:
T <- read.table(file="C:/User_R/dadostaper_edit.txt",head=T)
T
names(T)
#2) Ajustando modelos Taper:
di <- (T$di)
d <- (T$d)
did2 <- (di/d)^2
hih <- (T$hi/T$h)
did <- (di/d)
ldi <- log10(T$di)
ld <- log10(T$d)
lhhi <- log10(T$h-T$hi)
lh <- log10(T$h)
#a) Modelo 1 (md1): Demaerschalk (1971)
md1 <- nls(did2 ~ 10^(2*b0)*d^(2*b1-2)*(h-hi)^(2*b2)*h^(2*b3), data=T, start=list (b0=0.1, b1=0.1, b2=0.1, b3=0.1))
md1
summary(md1)
names(md1)
residuals(md1)
md1est <- (10^(0.1450)*d^(0.9114)*(T$h-T$hi)^(0.8386)*(T$h)^(-0.8356))
md1est
Emd1 <- (T$di) - md1est
Emd1
summary(Emd1)
EPmd1 <- (Emd1/(T$di))*100
EPmd1
summary(EPmd1)
mdpmd1 <- mean(EPmd1)
mdpmd1
c <- (sum(T$di))^2/282
c
sqtotal<- (sum(T$di^2))-c
sqtotal
sqresmd1<- sum(Emd1^2)
sqresmd1
glresmd1<- 278
glresmd1
n1<-282-1
r2md1<- 1-(sqresmd1/sqtotal)
r2md1
r2ajust<- 1-((1-r2md1)*(n1/glresmd1))
r2ajust
syxmd1<- sqrt(sqresmd1/glresmd1)
syxmd1
syxpmd1<- (syxmd1/mean(T$di))*100
syxpmd1