Discussion :

[CupOfAppleTea]

Bref aperçu de l'objectif:

Je dispose d'un jeu de données de température corporelle d'organismes en fonction de la température de leur milieu. J'ai calculé la différence entre la température des organismes et la température du milieu dans lequel ils évoluent afin de déterminer si ceux-ci sont plus/moins chauds que leur environnement. S'ils sont plus chauds, le delta de T° calculé (T°organisme - T°milieu) doit-être positif. Au contraire, s'ils sont plus froids, ce delta doit être négatif.

Je souhaiterais donc procéder à un test me permettant de déterminer si, globalement, la plupart de mes individus se situent en dessous de "0", ce qui signifierait qu'il sont plus froids. JE veux donc comparer d'un coté la distribution de mes individus (delta de température) et d'un autre coté une distribution théorique ou ceux-ci auraient pour valeur 0 (ni plus froids ni plus chaud, mais à la même température). C'est pour cela que je voulais utiliser le test de K-S qui permet de préciser l'hypothèse alternative, en l'occurrence pour moi ce serait "less" ici.

Voila, j'espère avoir été clair
Comment puis-je donc procéder dans mon cas? Quelle syntaxe pour le test K-S?

Merci pour votre aide!

[Sengar]

Bonjour,

Alors la tu me sembles un peu embrouillé donc on va essayer de redefinir les choses proprement et ensuite on verra si tu obtiens ce que tu veux

Bonjour,
Je souhaite comparer une distribution observée à une distribution théorique uniforme (y=0).

Déjà , sais tu ce qu'est une loi uniforme?
y=0 ce n'est pas une loi uniforme, ou alors une "loi uniforme" sur [0,0]


Si tu veux vraiment comparer tes observations à la distribution à valeur toujours nulle, alors pas besoin de se prendre la tête avec un test: tu regardes si tes observations sont toujours nulles. Si oui tu ne peux pas rejetter l'hypothèse que ca suit ta loi, et si il y a au moins une observation non nulle alors tu rejettes avec une pvalue de 0 et c'est tout!

Donc question: à quelle loi veut tu comparer ta distribution observée?


Si c'est une loi uniforme standard (dans [0,1]) alors tu fais :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

 ks.test(x, "punif")

et c'est fini.

Si c'est une autre loi, détaille nous laquelle et on te dira comment l'écrire sous R ou la simuler si il n'existe pas directement une fonction de base pour ca dans R.

Edit:
si c'est une loi uniforme sur [29.5,32.5] tu fais:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

 ks.test(x, "punif",29.5,32.5)

De plus, je souhaiterais que la loi uniforme suive une valeur de 0 des valeurs 29.5 à 32.5 pour x.

Cette phrase est incomprehensible.

[CupOfAppleTea]

Bref aperçu de l'objectif:

Je dispose d'un jeu de données de température corporelle d'organismes en fonction de la température de leur milieu. J'ai calculé la différence entre la température des organismes et la température du milieu dans lequel ils évoluent afin de déterminer si ceux-ci sont plus/moins chauds que leur environnement. S'ils sont plus chauds, le delta de T° calculé (T°organisme - T°milieu) doit-être positif. Au contraire, s'ils sont plus froids, ce delta doit être négatif.

Je souhaiterais donc procéder à un test me permettant de déterminer si, globalement, la plupart de mes individus se situent en dessous de "0", ce qui signifierait qu'il sont plus froids. JE veux donc comparer d'un coté la distribution de mes individus (delta de température) et d'un autre coté une distribution théorique ou ceux-ci auraient pour valeur 0 (ni plus froids ni plus chaud, mais à la même température). C'est pour cela que je voulais utiliser le test de K-S qui permet de préciser l'hypothèse alternative, en l'occurrence pour moi ce serait "less" ici.

Voila, j'espère avoir été clair et en effet je pense que je me suis bien égaré avec cette histoire de loi uniforme. Comment puis-je donc procéder dans mon cas?

[Sengar]

Ok j'ai compris.
Donc clairement la comparaison de la distribution des delta à une loi uniforme n'était pas adaptée à la problématique.

Ce que je ferais à ta place, c'est un test de moyenne.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

t.test(x=vecteurdesdelta, alternative ="two.sided")

Te dira si la moyenne est significativement différente de 0

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

t.test(x=vecteurdesdelta, alternative ="less")

Te dira si la moyenne est significativement inférieure à 0

[CupOfAppleTea]

C'est tout à fait ça!
En revanche le seul problème est que dans l'un de mes cas, ni la normalité ni l'homoscédasticité ne sont respectées. C'est pour ça que je m'étais au départ orienté vers K-S qui est non-paramétrique. Y aurait-il une version non-paramétrique du test de student me permettant de faire la même analyse?

[Sengar]

C'est tout à fait ça!
En revanche le seul problème est que dans l'un de mes cas, ni la normalité ni l'homoscédasticité ne sont respectées.

La normalité ce n'est pas un problème dès que tu as plus de ~30 observations avec le TCL.

L'homoscédasticité je vois pas le rapport puisque tu ne modélises rien et que tu n'as donc pas de residus.