Boa noite a todosEstou fazendo doutorado na UEL no curso de Ciência animal. Ainda sou novato no uso do R e também não tenho base estatística aprofundada em modelos não lineares.E preciso de um help, estou adaptando uma fórmula de modelo não linear do proc do SAS para o R, mas li o help, busquei dicas no site http://leg.ufpr.br/~paulojus/embrapa/Rembrapa/Rembrapase30.html#x32-20100030 e fui seguindo o modelo explicado por esse autor, mas tem algumas coisas que não entendi e não estou conseguindo rodar o modelo.Primeiramente explicarei resumidamente o que se trata o modelo. Esse modelo é para poder determinar cinco outras variáveis, através da produção de gás pela fermentação de bacteriais do líquido ruminal in vitro, simulando assim o rumen do animal. As variáveis que tenho em mãos são o volume de gás produzido e o tempo da leitura, a partir deles preciso determinar:VcnF = volume de gás produzido pelo carboidratos não fibrosos (mL de gás)kdcnf = taxa de degradação dos carboidratos não fibrosos (%/hora)VcF = volume de gás produzido pelo carboidratos fibrosos (mL de gás)kdcf = taxa de degradação dos carboidratos fibrosos (%/hora)L= tempo de latência (h), este é o tempo que as bactérias gastam para começar a degradar o substrato.Mas ao rodar o modelo abaixo aparece o seguinte erro: fit30=nls(Y~ (VcnF/(1+exp(2-4*kdcnf*I(L-Tempo)))+ VcF/I(1+exp(2-4*kdcf*(L-Tempo)))), + data=cra30, + start=list(VcnF = 10,kdcnf = 0.01,L =0.1,VcF = 50,kdcf =0.01, Tempo=1), + alpha= (0.05))Error in qr.qty(QR, resid) : NA/NaN/Inf in foreign function call (arg 5)In addition: Warning message:In Ops.factor(lhs, rhs) : - not meaningful for factors Aproveitando o ensejo gostaria de perguntar como faço para limitar os parâmetros pois no proc do SAS está assim:VcnF = 10 to 200 by 50 kdcnf = 0.01 to 0.8 by 0.05 L = 0.1 to 5 by 0.2 VcF = 50 to 150 by 30 kdcf = 0.01 to 0.8 by 0.05; bounds L >=0; Caso alguém tenha materiais relativos a modelos não linear para que possa aprofundar melhor, também contribuirá bastante. Agradeço a todos, e espero que alguém possa me ajudar! Eduardo Lucas Terra PeixotoDiscente de doutorado - Ciência AnimalUniversidade Estadual de Londrina Cel. (43) 9653-5574
_______________________________________________ R-br mailing list [email protected] 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.
