changepoint

[pinson laura]

Bonjour bonjour...

Je reviens encore vers vous car de nouveau un pb

Je souhaite segmenter ma série de données avec le package changepoint mais je n'arrive pas à l'utiliser...

Code : Tout sélectionner

cpt.mean(Moy_inte,method="PELT", penalty="Manual", pen.value="5*log(n)")
plot(seg1,cpt.width=2)

Mais j'ai un msg d'erreur me disant :

Error in checkData(data) : object 'Moy_inte' not found
Moy_inte c'est le nom de mon fichier, mon fichier est composé de 3 colonnes : V1=X , V2=Y et V3=data (mes températures)

Avez vous une idée???



Merci merci

[Dominique Soudant]

L'objet auquel tu références n'existe pas. C'est peut être le nom de ton fichier mais si tu ne l'as pas lu et que tu n'as pas organisé son contenu dans un objet R (e.g. vecteur matrice data frame) la fonction ne peut pas utiliser directement un nom de fichier ...

Ton problème ne porte pas tant sur changepoint que sur le fonctionnement de R.
?

[Serge Rapenne]

Pour compléter la réponse de Dominique :

Code : Tout sélectionner

?cpt.mean

Usage

cpt.mean(data,penalty="MBIC",pen.value=0,method="AMOC",Q=5,test.stat="Normal",class=TRUE,param.estimates=TRUE,minseglen=1)

Arguments

data
A vector, ts object or matrix containing the data within which you wish to find a changepoint. If data is a matrix, each row is considered a separate dataset.

[pinson laura]

je dois donc tracer mon graphique avant?


je suis totalement perdu....

[Dominique Soudant]

Un graphique n'a rien à voir là dedans ...
Mettons que ton fichier contienne plusieurs colonnes dont les noms pourraient être "date" "variable1" "variable2" etc.
Tu charges ton fichier dans R avec un read.delim par exemple, si c'est approprié. read.delim te fait un dataframe que tu vas appeler (avec plein d’à-propos et d'originalité) MonDataFrame (ou alors tu peux aussi l'appeler toto, c'est une autre école de nommage). Et après tu fais kkchose du genre cptmean(MonDataFrame$variable1 ... ) si tu veux appliquer cpt.mean sur la variable1. Exemple perso:

Code : Tout sélectionner

cpt.mean(SansNA$ets
   ,method="AMOC" # at most one change
   ,penalty="None"
   ,class=FALSE)

[pinson laura]

Voila, j'ai fais ma courbe...mais j'ai toujours des souci

Code : Tout sélectionner

x1<-ggplot(grrr,aes(V1, V2))
x1+geom_line()
seg1=cpt.mean(x1,method="PELT", penalty="Manual", pen.value="5*log(n)")
plot(seg1,cpt.width=2)

et voila mon msg d'erreur :

Error in checkData(data) : Only numeric data allowed
> plot(seg1,cpt.width=2)
Error in plot(seg1, cpt.width = 2) :
erreur d'évaluation de l'argument 'x' lors de la sélection d'une méthode pour la fonction 'plot' : Error: object 'seg1' not found

je ne comprends absolument pas.


J'ai essayé ta méthode Serge, mais c'est pas vraiment mieux :'(

[pinson laura]

yyyyyyyyyyyyyyyeeeeeeeeeeeeeessssssssssssssssssssssssssssssssssssssssssssss



1000 mercis dominique!
TU serais en face de moi, je t'aurai pris ds mes bras!!!

Je peux terminer ma thèse, tranquillement!!!!


10000000mercis pour ta patience!!!!!


HEUREUSE!

J'adore ce forum!
Vous etes tous aussi adorable les uns que les autres!