Bonjour !

Je travaille actuellement avec le package 'segmented' de Muggeo pour réaliser une régression linéaire segmentée (deux segments). J'aimerais imposer que la pente de gauche soit nulle. J'ai suivi les instructions du package (http://dssm.unipa.it/vmuggeo/segmentedRnews.pdf) pour écrire le code suivant :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
  fit.glm <- glm(y~x) 
  fit.seg <- segmented(fit.glm, seg.Z=~x,psi=0.3) 
  fit.glm <-update(fit.glm,.~.-x) 
  fit.seg1 <-update(fit.seg)
Les deux premières lignes sont pour la régression segmentée sans imposition, et les 2 dernières sont supposées adapter la pente de gauche. Malheureusement, cela ne fonctionne pas... Je dois avouer que j'ai suivi un peu aveuglément un script modèle parce que je ne comprend pas le .~.-x dans la fonction update.

Merci beaucoup d'avance pour votre aide !