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 - ivanalaman@yahoo.com.br