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 -
[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.
