Caros Membros,

           Utilizava com frequência o Matlab para resolver integrais por não usar essas operações matemáticas rotineiramente no R. Mas agora resolvi fazer todos os cálculos em um mesmo script e estou tendo problemas para realizar a integral de um modelo não linear para calculo da área abaixo da curva, sendo meu CRM:

#Banco de dados simulado
N=100#tamanho amostral
vI<-c(0.1,0.7,1.3,2,3,4,5,6,7,13)
RES=NULL
vI<-rep(vI,N)
di<-rnorm(N,15)
ap<-rnorm(N,10)
th<-rnorm(N,50)
RES=rbind(Nsim,vI,di,ap,th)
RES_S<-as.data.frame(t(RES))
#

#Ajuste de modelo não linear com o pacote minpack.lm
library(minpack.lm)
d1_1_S<- nlsLM(di ~ ap*b0*(1+b1*log(1-b2*vI^b3*th^-b3)),
                   start=list(b0=1.3333,
                              b1=0.6667,
                              b2=-0.3333,
                              b3=0.3333),
                   data = RES_S)
summary(d1_1_S) ## Fazendo de conta que os coeficientes são significativos

#Plotando o gráfico
newdata<-RES_S
newdata$ap <- mean(RES_S$ap)
newdata$th <- mean(RES_S$th)
newdata$yhat <- predict(d1_1_S, newdata, type = "response")
newdata <- newdata[order(newdata$vI), ]
plot(di ~ vI, data = RES_S, xlab="vI",ylab="di" )
lines(x = newdata$vI, y = newdata$yhat )
CI<-predict(d1_1_S, interval = "confidence")
#--

Agora eu gostaria de resolver a integral de di ~ ap*b0*(1+b1*log(1-b2*vI^b3*th^-b3),
para calcular a área abaixo da curva para meus limites de vI para cada intervalo de 0,1 a 0.7
até 7 a 13. Alguém poderia me dar uma luz?
Obrigado

--
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO)
e-mails:alexandresantosbr@yahoo.com.br
        alexandre.santos@cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
OrcID: orcid.org/0000-0001-8232-6722   -   ResearcherID: A-5790-2016
Researchgate: www.researchgate.net/profile/Alexandre_Santos10
LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635
Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/
======================================================================

_______________________________________________
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 fornea cdigo mnimo reproduzvel.