<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />Salut tout le monde. J'essaie d'exécuter un Mann-Kendall pour tester si une tendance temporelle est significative. Impossible de le trouver dans le manuel JMP ou quoi que ce soit en ligne. Quelqu'un a des idées?
Merci!
merci MS et Dave ! Je ne savais pas non plus si le Tau de Kendall et le test de Mann-Kendall étaient identiques (et/ou liés). Merci pour l'info sur le test multivarié et ensuite pour le script aussi !
Même si les définitions de la statistique du test de Mann-Kendall et du Tau b de Kendall sont souvent décrites en des termes différents (du point de vue des séries chronologiques et de la corrélation, respectivement), elles sont essentiellement les mêmes pour autant que je sache.
Voici un script que j'ai écrit il y a assez longtemps. Il calcule et affiche la pente de Theil et effectue le test de Mann-Kendall. Une ligne de régression linéaire paramétrique est affichée dans le graphique uniquement à des fins de comparaison. J'ai prévu de généraliser le script en incluant une fonction "By" ainsi que la possibilité de tenir compte des modèles saisonniers. Mais je n'ai pas encore eu le temps de réaliser ces plans. J'espère que cela vous donnera des idées.
dt = Tableau de données actuel () ;
dlg = Boîte de dialogue Colonne (
"Estimer la pente de Theil et Kendall Tau à partir du tableau de données " ||
Carboniser ( dt << obtenir le nom ) ,
Oui = ListeCol ( "Sélectionnez Y" ,
maxcol ( 1 ) ,
mincol ( 1 ) ,
Type de données ( numérique )
) ,
X = ListeCol ( "Sélectionner X" ,
maxcol ( 1 ) ,
mincol ( 1 ) ,
Type de données ( numérique )
) ,
Groupe = ListeCol ( "En (ne travaillant pas)" , maxcol ( 1 ) )
) ;
colY = Colonne ( dt , dlg [ "Oui" ] ) ;
colX = Colonne ( dt , dlg [ "X" ] ) ;
colGroupe = Colonne ( dt , dlg [ "Groupe" ] ) ;
//Calcul de la pente de Theil
n = N lignes ( dt ) ;
pentes = J ( n , n , . ) ;
Pour ( je = 1 , je < n , je ++,
Pour ( j = je + 1 , j <= n , j ++,
pentes [ je , j ] = ( colY [ j ] - colY [ je ]) / ( colX [ j ] - colX [ je ])
)
) ;
theil = quantile ( 0,5 , pentes ) ;
//Tau b de Kendall pour tester si la pente diffère de zéro
m.v. = Multivarié (
Oui ( colY , colX ) ,
Matrice de nuage de points ( 0 ) ,
Corrélations multivariées ( 0 ) ,
Tau de Kendall ( 1 )
) ;
//Tracer les données et ajouter la régression linéaire et la pente de Theil
bivouac = Bivarié (
Oui ( colY ) ,
X ( colX ) ,
Ligne d'ajustement ( { Couleur de la ligne ( "Rouge" )} ) ,
Ajustement spécial ( Pente ( theil ) , { Couleur de la ligne ( "Vert" )} ) ,
invisible
) ;
biv = bivouac << rapport ;
biv [ "Ajustement linéaire" ] << fermer tout comme ça ;
biv [ Boîte de contour ( sept )] << définir le titre ( "Pente de Theil" ) ;
db = Signaler ( m.v. )[ "Multivarié" ] << parent ;
db << ajouter ( biv ) ;
db [ Zone de texte ( 4 )] << Définir le texte ( "Pente de Theil" ) ;
biv << Fermer la fenêtre ;
Je ne sais pas dans quelle mesure (le cas échéant) le test de Mann-Kendall est lié à la statistique Kendall Tau pour la corrélation non paramétrique.Ces derniers sont accessibles sous la plate-forme multivariée : Analyser>Méthodes multivariées, puis dans le rapport multivarié, sélectionnez Corrélations non paramétriques à partir du hotspot.
Si les tests sont complètement différents, je serais intéressé par toutes les références que vous avez pour cette méthode.J'envisage de créer un complément qui implémente un certain nombre de méthodes (par exemple, Hurst Index) qui ont été développées pour le domaine de l'hydrologie
Salutations
David
