Produção de gás


		

# para exemplificar vamos utilizar os dados simulados abaixo  (produção de gas em função tempo)


		

tempo=c(0,12,24,36,48,60,72,84,96,108,120,144,168,192)
gas=c(0.002,3.8,8,14.5,16,16.5,17,17.4,17.9,18.1,18.8,19,19.2,19.3)

 

# criando uma tabela

dados=data.frame(tempo, gas)

 

# ajustando o modelo com a função nlsfit() do pacote easynls

# note na programação abaixo que o argumento model define o modelo e start define os valores iniciais para o processo de estimação


nlsfit(dados, model=13, start=c(19,4,0.025,0.004,5))

 

 

# para observar graficamente o modelo é só utilizar a função nlsplot() do pacote easynls

nlsplot(dados, model=13, start=c(19,4,0.025,0.004,5))

 

# veja abaixo o gráfico