
Walmes, meus parabéns pela sua didática. Impecável!!! Abs Em 8 de abril de 2011 17:35, Walmes Zeviani <walmeszeviani@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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@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