<div><span style="font-size: medium; color: #0000ff;">Não estou conseguindo rodar porque tem um objeto Nsim no último comando abaixo que ele não reconhece!</span></div>
<div dir="ltr"><span style="font-size: 12.8px;">N=100#tamanho amostral</span><br style="font-size: 12.8px;" /><span style="font-size: 12.8px;">vI<-c(0.1,0.7,1.3,2,3,4,5,6,7,13)</span><br style="font-size: 12.8px;" /><span style="font-size: 12.8px;">RES=NULL</span><br style="font-size: 12.8px;" /><span style="font-size: 12.8px;">vI<-rep(vI,N)</span><br style="font-size: 12.8px;" /><span style="font-size: 12.8px;">di<-rnorm(N,15)</span><br style="font-size: 12.8px;" /><span style="font-size: 12.8px;">ap<-rnorm(N,10)</span><br style="font-size: 12.8px;" /><span style="font-size: 12.8px;">th<-rnorm(N,50)</span><br style="font-size: 12.8px;" /><span style="font-size: 12.8px;">RES=rbind(Nsim,vI,di,ap,th)</span></div>
<div dir="ltr"> </div>
<div dir="ltr"> </div>
<div dir="ltr">Alexandre,
<div> </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> </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> </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> </div>
<div>HTH</div>
<div>--</div>
<div>Cesar Rabak</div>
<div> </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="../../../undefined//compose?to=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,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^-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 /> #--</blockquote>
<br /> Agora eu gostaria de resolver a integral de di ~ ap*b0*(1+b1*log(1-b2*vI^b3*th^-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 /> ======================================================================<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="../../../undefined//compose?to=e-mails" target="_blank">e-mails:alexandresantosbr@yahoo.com.br</a><br /> <a href="../../../undefined//compose?to=alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.br</a><br /> Lattes: <a href="http://lattes.cnpq.br/1360403201088680" rel="noreferrer" target="_blank">http://lattes.cnpq.br/1360403201088680</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/Alexandre_Santos10</a><br /> LinkedIn: <a href="http://br.linkedin.com/in/alexandre-dos-santos-87961635" rel="noreferrer" target="_blank">br.linkedin.com/in/alexandre-dos-santos-87961635</a><br /> Mendeley:<a href="http://www.mendeley.com/profiles/alexandre-dos-santos6/" rel="noreferrer" target="_blank">www.mendeley.com/profiles/alexandre-dos-santos6/</a><br /> ======================================================================<br /> <br /> _______________________________________________<br /> R-br mailing list<br /> <a href="../../../undefined//compose?to=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-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-guia</a>) e fornea cdigo mnimo reproduzvel.</blockquote>
</div>
</div>
<div>_______________________________________________<br />R-br mailing list<br />R-br@listas.c3sl.ufpr.br<br /><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br />Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne�a c�digo m�nimo reproduz�vel.</div>