<font face="trebuchet ms,sans-serif">Ivan,<br><br>Dá pra fazer assim, embora eu acredito que exista um jeito mais simples, <br><br><span style="font-family: courier new,monospace;">da <- expand.grid(A=factor(LETTERS[1:3]), x=1:10)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">da$y <- 10/(1+exp(-(da$x-5)/1))+rnorm(da$x,0,0.1)</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">n0 <- nls(y~As[A]/(1+exp(-(x-x0[A])/S[A])), data=da,</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">          start=list(As=c(10,10,10), x0=c(5,5,5), S=c(1,1,1)))</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">str(n0)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">n0$data</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">eval(n0$data)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">all.vars(summary(n0)$formula)[1]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">eval(n0$data)[,all.vars(summary(n0)$formula)[1]]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">mean(eval(n0$data)[,all.vars(summary(n0)$formula)[1]])</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">fun <- function(n0){</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">  mean(eval(n0$data)[,all.vars(summary(n0)$formula)[1]])</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">}                         </span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">fun(n0)</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;">À disposição.<br>Walmes.<br><br clear="all"></font><span style="font-family: trebuchet ms,sans-serif;">==========================================================================</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">Walmes Marques Zeviani</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">fone: (+55) 41 3361 3573</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">VoIP: (3361 3600) 1053 1173</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a></span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">twitter: @walmeszeviani</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">linux user number: 531218</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">==========================================================================</span><br>