<div dir="ltr">Alexandre,<div><br></div><div>Considerando que os dados da sua curva "empírica" já estão dados, é muito mais interessante do ponto de vista matemático e prático (inclusive no R) fazer a integral por meios numéricos, especialmente por você deixar claro que deseja a integral definida (« área abaixo da curva para meus limite. . . »).</div><div><br></div><div>Inclusive há um pacote (caTools [e parece a pracma, que ñ conheço]) que tem a função 'clone' do MATLAB trapz.</div><div><br></div><div>se você fizer a pesquisa "r numerical integration" terá bastante material para resolver esse problema, inclusive com soluções mais sofisticadas.</div><div><br></div><div>HTH</div><div>--</div><div>Cesar Rabak</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-06-14 10:43 GMT-03:00 ASANTOS via R-br <span dir="ltr"><<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Caros Membros,<br>
<br>
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:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
#Banco de dados simulado<br>
N=100#tamanho amostral<br>
vI<-c(0.1,0.7,1.3,2,3,4,5,6,7,<wbr>13)<br>
RES=NULL<br>
vI<-rep(vI,N)<br>
di<-rnorm(N,15)<br>
ap<-rnorm(N,10)<br>
th<-rnorm(N,50)<br>
RES=rbind(Nsim,vI,di,ap,th)<br>
RES_S<-as.data.frame(t(RES))<br>
#<br>
<br>
#Ajuste de modelo não linear com o pacote minpack.lm<br>
library(minpack.lm)<br>
d1_1_S<- nlsLM(di ~ ap*b0*(1+b1*log(1-b2*vI^b3*th^<wbr>-b3)),<br>
start=list(b0=1.3333,<br>
b1=0.6667,<br>
b2=-0.3333,<br>
b3=0.3333),<br>
data = RES_S)<br>
summary(d1_1_S) ## Fazendo de conta que os coeficientes são significativos<br>
<br>
#Plotando o gráfico<br>
newdata<-RES_S<br>
newdata$ap <- mean(RES_S$ap)<br>
newdata$th <- mean(RES_S$th)<br>
newdata$yhat <- predict(d1_1_S, newdata, type = "response")<br>
newdata <- newdata[order(newdata$vI), ]<br>
plot(di ~ vI, data = RES_S, xlab="vI",ylab="di" )<br>
lines(x = newdata$vI, y = newdata$yhat )<br>
CI<-predict(d1_1_S, interval = "confidence")<br>
#--<br>
</blockquote>
<br>
Agora eu gostaria de resolver a integral de di ~ ap*b0*(1+b1*log(1-b2*vI^b3*th^<wbr>-b3),<br>
para calcular a área abaixo da curva para meus limites de vI para cada intervalo de 0,1 a 0.7<br>
até 7 a 13. Alguém poderia me dar uma luz?<br>
Obrigado<br>
<br>
-- <br>
==============================<wbr>==============================<wbr>==========<br>
Alexandre dos Santos<br>
Proteção Florestal<br>
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso<br>
Campus Cáceres<br>
Caixa Postal 244<br>
Avenida dos Ramires, s/n<br>
Bairro: Distrito Industrial<br>
Cáceres - MT CEP: 78.200-000<br>
Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO)<br>
<a href="mailto:e-mails%3Aalexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@yaho<wbr>o.com.br</a><br>
<a href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.<wbr>br</a><br>
Lattes: <a href="http://lattes.cnpq.br/1360403201088680" rel="noreferrer" target="_blank">http://lattes.cnpq.br/13604032<wbr>01088680</a><br>
OrcID: <a href="http://orcid.org/0000-0001-8232-6722" rel="noreferrer" target="_blank">orcid.org/0000-0001-8232-6722</a> - ResearcherID: A-5790-2016<br>
Researchgate: <a href="http://www.researchgate.net/profile/Alexandre_Santos10" rel="noreferrer" target="_blank">www.researchgate.net/profile/A<wbr>lexandre_Santos10</a><br>
LinkedIn: <a href="http://br.linkedin.com/in/alexandre-dos-santos-87961635" rel="noreferrer" target="_blank">br.linkedin.com/in/alexandre-d<wbr>os-santos-87961635</a><br>
Mendeley:<a href="http://www.mendeley.com/profiles/alexandre-dos-santos6/" rel="noreferrer" target="_blank">www.mendeley.com/prof<wbr>iles/alexandre-dos-santos6/</a><br>
==============================<wbr>==============================<wbr>==========<br>
<br>
______________________________<wbr>_________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi<wbr>-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-g<wbr>uia</a>) e fornea cdigo mnimo reproduzvel.</blockquote></div><br></div>