Estou tentando rodar esta rotina no R mas esta dando este erro e tambem ele não esta estimando b2, poderiam me ajudar? lrp.fit<-nls(cv ~ (b0 + b1*x + b2*I(x^2))*(x<=x0p) +(b0 + b1*x0p + b2*x0p^2)*(x>x0p), start=list(b0=10, b1=7, b2=15, x0p=5), trace=F) Erro em nls(cv ~ (b0 + b1 * x + b2 * I(x^2)) * (x <= x0p) + (b0 + b1 * : gradiente singular
Juraci, Seu problema não é reproduzível! Seus chutes iniciais são cartesianamente plausíveis, ou seja, ao plotar uma curva com esses valores iniciais ela passa por entre os dados? Troque os chutes iniciais porque os atuais estão criando uma matriz de derivadas primeiras (F) que em um passo (F'F) não é invertível. Na formula da nls() você não precisa usar I(x^2), só x^2 mesmo. À 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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (2)
-
Juraci Mendes Moreira -
Walmes Zeviani