Integral de um modelo não linear para calculo da área abaixo da curva
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/ ======================================================================
Alexandre, 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. . . »). Inclusive há um pacote (caTools [e parece a pracma, que ñ conheço]) que tem a função 'clone' do MATLAB trapz. se você fizer a pesquisa "r numerical integration" terá bastante material para resolver esse problema, inclusive com soluções mais sofisticadas. HTH -- Cesar Rabak 2017-06-14 10:43 GMT-03:00 ASANTOS via R-br <r-br@listas.c3sl.ufpr.br>:
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.
participantes (2)
-
ASANTOS -
Cesar Rabak