Caro Ivan com 6 pontos? voce parece estar estimando um modelo com mais parametros que dados, especialmente se estiver estimando tb os pontos de quebra! On Fri, 22 Jul 2011, Ivan Bezerra Allaman wrote:
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