
22 Jul
2011
22 Jul
'11
18:37
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