<div>Desculpem-me pela mancada, segue email anterior, saída do dput() e CMR... </div><div><br></div><div><span style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">Prezados colegas, tenho alguns dados de produção de metano ao longo do tempo em intervalos de 12 minutos. Gostaria da opinião de vocês sobre a melhor forma de integrar a área sob a curva formada. Já pesquisei algumas coisas e acredito que usar a função integrate() sobre o vetor gerado com as funções approxfun() ou splinefun() podem resolver o problema ao ligarem os pontos, embora redundem em resultados diferentes e não contemplem exatamente todos os pontos. Segue o banco de dados e as funções usadas por mim... </span></div>
<div><span style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)"><br></span></div><div>### Banco de dados</div><div><div>data<-structure(list(Sample = 1:100, CH4A = c(0.02492, 0.0272, 0.03284, </div>
<div>0.04048, 0.04256, 0.0474, 0.05128, 0.05348, 0.05292, 0.05156, </div><div>0.05532, 0.05608, 0.05152, 0.05234, 0.05222, 0.05266, 0.056, </div><div>0.05008, 0.05038, 0.05148, 0.0485, 0.04748, 0.04272, 0.04294, </div><div>
0.04132, 0.03824, 0.03672, 0.04056, 0.03592, 0.03492, 0.03818, </div><div>0.03628, 0.03492, 0.03588, 0.03196, 0.0338, 0.03456, 0.03284, </div><div>0.03684, 0.03282, 0.03116, 0.03082, 0.03072, 0.03488, 0.0349, </div><div>0.035, 0.03516, 0.03256, 0.02988, 0.03408, 0.02854, 0.03392, </div>
<div>0.02844, 0.0327, 0.031, 0.03074, 0.0308, 0.02808, 0.02796, 0.02822, </div><div>0.04016, 0.03704, 0.04814, 0.04636, 0.04824, 0.0514, 0.0506, </div><div>0.05602, 0.0623, 0.06436, 0.06132, 0.06396, 0.0612, 0.05466, </div>
<div>0.05832, 0.06704, 0.06776, 0.06076, 0.05698, 0.06096, 0.06566, </div><div>0.06282, 0.05626, 0.05456, 0.05716, 0.05634, 0.05258, 0.05564, </div><div>0.05608, 0.05284, 0.05228, 0.05008, 0.05322, 0.05018, 0.04752, </div>
<div>0.04478, 0.04732, 0.04524, 0.03788, 0.04192)), .Names = c("Sample", </div><div>"CH4A"), row.names = c(NA, 100L), class = "data.frame")</div></div><div><br></div><div><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
#### Plotagem -----</div><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">plot(CH4A~Sample,data,type='l')</div><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
#### integrando com approxfun()</div><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">f2 <- approxfun(data$Sample, data$CH4A)</div><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
integrate(f2, min(data$Sample), max(data$Sample),subdivisions=1000)</div><span style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)"></span><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
#### integrando com splinefun()</div><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">f3 <- splinefun(data$Sample, data$CH4A) </div><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
integrate(f3, min(data$Sample), max(data$Sample),subdivisions=1000)</div><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">#### plotagem final --------</div>
<div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">plot(data$CH4A~data$Sample, type = "n") </div><span style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)"></span><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
curve(f2(x), add = TRUE, col=1) </div><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">curve(f3(x), add = TRUE, col=3) </div><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
f4 <- splinefun(data$Sample, data$CH4A, method ="natural") </div><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">curve(f4(x), add = TRUE, col=4)</div>
<div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">### Percebam que algumas variações bruscas não são contempladas pelas curvas ajustadas, o que provavelmente afeta a área.</div>
</div><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)"><br></div><div style="font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
<br></div>