Boa tarde senhores!!

Estou tendo dificuldades ao ajustar uma regressão segmentada com a library 
''segmented'. Não entendo qual a dificuldade que a função está encontrando para 
ajustar uma regressão aos meus dados. Creio que os breaks pontos chutados não 
estão tão fora da realidade. Segue abaixo um CMR real.

Y <- c(18.335, 19.456, 19.383, 19.542, 20.718, 21.531)
X <- c(0,3,6,9,12,15)
plot(X,Y,type='b')
mod <- lm(Y ~ X)
library(segmented)
mod.seg <- segmented(mod, seg.Z=~X, psi=list(X=c(3,8)))

Se alguém puder ajudar, agradeço.

(S,f,P)
Allaman


 
Ivan Bezerra Allaman 
Doutor em Zootecnia/UFLA
email e msn - [email protected] 
_______________________________________________
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.

Responder a