Extrair dados de ajuste
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 (TIM)* *"**Nós somos aquilo que fazemos repetidamente. Excelência, então, não é um modo de agir, mas um hábito." Aristóteles*
É isso que quer: m.1 <- summary(model) names(m.1) m.1$parameters[,1] Em 05-07-2016 15:21, Pedro Campelo via R-br escreveu:
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 (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.
-- Atenciosamente Felipe E. Barletta Mendes Estatístico(UFPR) - Conre3 9766-A Mestrando em Bioestatística(UEM) +55 (41)-92077191 +55 (41)-33287216
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 ###############################################################
Opa, faltou carregar o pacote: library(broom). 2016-07-05 15:50 GMT-03:00 Jônatan <jdtatsch@gmail.com>:
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 ###############################################################
-- ############################################################### ## 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 ###############################################################
participantes (3)
-
Felipe -
Jônatan -
Pedro Campelo