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: [email protected] twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
