[R-br] Função não reconhece variável extraída de um modelo!

[Ivan Bezerra Allaman]

Valeu Walmes! O que importa é que resolveu! Obrigado de novo.

(S,f,P)
Allaman
 

\begin{signature}
<<>>=
Prof. Dr. Ivan Bezerra Allaman
Universidade Estadual de Santa Cruz
Departamento de Ciências Exatas e Tecnológicas
Ilhéus/BA - Brasil
Fone: +55 73 3680-5076
E-mail: ivanalaman em yahoo.com.br/ivanalaman em gmail.com
@
\end{signature}


________________________________
De: Walmes Zeviani <walmeszeviani em gmail.com>
Para: r-br em listas.c3sl.ufpr.br; Ivan Bezerra Allaman <ivanalaman em yahoo.com.br>
Enviadas: Quinta-feira, 29 de Setembro de 2011 15:12
Assunto: Re: [R-br] Função não reconhece variável extraída de um modelo!


Ivan,

Dá pra fazer assim, embora eu acredito que exista um jeito mais simples, 

da <- expand.grid(A=factor(LETTERS[1:3]), x=1:10)
da$y <- 10/(1+exp(-(da$x-5)/1))+rnorm(da$x,0,0.1)

n0 <- nls(y~As[A]/(1+exp(-(x-x0[A])/S[A])), data=da,
          start=list(As=c(10,10,10), x0=c(5,5,5), S=c(1,1,1)))

str(n0)
n0$data
eval(n0$data)
all.vars(summary(n0)$formula)[1]
eval(n0$data)[,all.vars(summary(n0)$formula)[1]]
mean(eval(n0$data)[,all.vars(summary(n0)$formula)[1]])

fun <- function(n0){
  mean(eval(n0$data)[,all.vars(summary(n0)$formula)[1]])
}                         
fun(n0)

À 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
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110929/7e84b71f/attachment.html>