[R-br] RES.: Função sigmóide! Dados!

[Guilherme Moraes Ferraudo]

Walmes, meus parabéns pela sua didática. Impecável!!!

Abs

Em 8 de abril de 2011 17:35, Walmes Zeviani <walmeszeviani em gmail.com> escreveu:
> Thiago,
>
> Segue abaixo um CMR no estilo do que deveria ser (pelo menos para mim) as
> trocas de ajuda dentro da lista. Esclarecimentos estão comentados dentro do
> código.
>
> #------------------------------------------------------------------------------------------
> # exemplo repruzível (exemplo para os demais inscritos da lista)
> #------------------------------------------------------------------------------------------
> # dados hopesdados pelo site http://www.datafilehost.com/
> # link para a página dos dados
> http://www.datafilehost.com/download-a7e94384.html
> # link do download http://www.datafilehost.com/get.php?file=a7e94384
> # o melhor é fazer com que os dados sejam baixados na leitura
> # esse procedimento não funcionou para os seus dados. porque?
>
> da <- read.table("http://www.datafilehost.com/get.php?file=688e5891",
> header=TRUE) # falha
>
> #------------------------------------------------------------------------------------------
> # hospedei os dados no site do LEG, assim eu posso baixar direto da web
> # temos que arrumar algo assim para que todos possam usar
>
> da <- read.table("http://www.leg.ufpr.br/~walmes/docs/umi.txt", header=TRUE)
> # falha
> str(da)
>
> #------------------------------------------------------------------------------------------
> # diagnose gráfica
>
> plot(UR~UEQ, da)
>
> #------------------------------------------------------------------------------------------
> # pelo jeito, você fixou a UR para obter os valores de UEQ. então não seria
> UEQ em função
> # de UR? vou continuar a análise para mostrar os procedimentos, mas pense
> sobre isso.
>
> #------------------------------------------------------------------------------------------
> # modelo logístico não vai se ajustar aos seus dados porque eles não tem
> informação
> # para ajustar uma curva sigmoide. os dados tem padrão de crescimento
> assintótico
> # por isso ajustar um modelo que alçance essa comportamento
>
> A <- 100   # mudar até obter bons chutes
> B <- 0.11  # mudar até obter bons chutes
> plot(UR~UEQ, da)
> curve(A*(1-exp(-B*x)), add=TRUE)
>
> #------------------------------------------------------------------------------------------
> # ajuste do modelo com os chutes
>
> n0 <- nls(UR~A*(1-exp(-B*UEQ)), data=da,
>           start=list(A=100, B=0.11))
> summary(n0)
>
> #------------------------------------------------------------------------------------------
> # curva ajustada
>
> A <- coef(n0)["A"]
> B <- coef(n0)["B"]
> plot(UR~UEQ, da)
> curve(A*(1-exp(-B*x)), add=TRUE, col=2)
>
> #------------------------------------------------------------------------------------------
> # Código mínimo reproduzível: tudo se resolve em questão de segundos!
> #------------------------------------------------------------------------------------------
>
> À disposição.
> Walmes.
>
> ==========================================================================
> Walmes Marques Zeviani
> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
> Departamento de Estatística - Universidade Federal do Paraná
> fone: (+55) 41 3361 3573
> VoIP: (3361 3600) 1053 1173
> e-mail: walmes em ufpr.br
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>



-- 
Guilherme Moraes Ferraudo
Jaboticabal/Campinas - SP
http://lattes.cnpq.br/2096118558794430