Uma alternativa para obter as informações do ajuste do modelo em um data frame. model<-nls(R~k*M^n, start=list(k=0.5,n=1), *data = df*) stats <- tidy(model) # term estimate std.error statistic p.value #1 k 11.53004 0.58948562 19.55949 1.808707e-10 #2 n 0.31174 0.01158954 26.89839 4.284298e-12 class(stats) "data.frame" stats[, 1:2] # term estimate #1 k 11.53004 #2 n 0.31174 2016-07-05 15:21 GMT-03:00 Pedro Campelo via R-br <r-br@listas.c3sl.ufpr.br> :
Pessoal,
Estou rodando um script e preciso extrair os dados dos coeficientes ajustados para um modelo. No meu exemplo, quero extrair e separar os valores "k" e "n".
*> data.frame(A)* * M R* *1 10 24.61* *2 20 29.41* *3 30 33.07* *4 40 37.37* *5 50 38.22* *6 60 41.14* *7 70 41.61* *8 80 43.41* *9 90 47.18* *10 100 48.77* *11 110 51.02* *12 120 52.95* *13 130 52.98* *14 140 52.92*
*> model<-nls(R~k*M^n, start=list(k=0.5,n=1))* *> summary(model)*
Obrigado
*----*
*Pedro Henrique Campelo Felix* *Bachelor of Food Engineering (UNI-BH)* *Msc. Chemical Engineering (UFMG)* *PhD Student of Food Science (UFLA)*
*+55 (31) 97560-1545 <%2B55%20%2831%29%2097560-1545> (TIM)* *"**Nós somos aquilo que fazemos repetidamente. Excelência, então, não é um modo de agir, mas um hábito." Aristóteles*
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- ############################################################### ## Jônatan Dupont Tatsch ## Professor do Departamento de Física ## Centro de Ciências Exatas e Naturais (CCNE) ## Universidade Federal de Santa Maria - UFSM ## Faixa de Camobi, Prédio 13 - Campus UFSM - Santa Maria, RS, Brasil - 97105-900 ## Telefone: +55(55)33012083 ## www.ufsm.br/meteorologia ###############################################################