Boa tarde senhores!
Gostaria de saber como posso fazer uma função ler uma variável a partir de um modelo, nls por exemplo. Vamos supor um CMR.
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)))
fun <-
function(modelo){
mean(all.vars(summary(modelo)$formula)[1])
}
Percebam que eu quero que a função calcule a média de uma variável, cuja variável foi declarada num determinado modelo. A função acima não dá certo, pois a função entende o seguinte:
mean("VARIÁVEL")
Para dar certo, eu acho, a variável não poderia estar entre aspas. Enfim, alguém tem alguma idéia de como eu posso fazer isso?
Desde já agradeço a todos.
(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@yahoo.com.br/ivanalaman@gmail.com
@
\end{signature}