Derniers messages sur Zeste de Savoirhttps://zestedesavoir.com/forums/2022-06-20T07:27:32+02:00Les derniers messages parus sur le forum de Zeste de Savoir.Traiter des données ordinales en fonction du temps dans R, message #2437082022-06-20T07:27:32+02:00Moté/@Mot%C3%A9https://zestedesavoir.com/forums/sujet/16323/traiter-des-donnees-ordinales-en-fonction-du-temps-dans-r/?page=1#p243708<p>Bonjour,</p>
<p>Est-ce que ce ne serait pas le but des médianes et autres quartiles ? Dans ce cas, il faudrait peut-être tracer les boxplot (boîtes à moustaches) de tes données.</p>Traiter des données ordinales en fonction du temps dans R, message #2436982022-06-18T13:00:26+02:00JuRDT/@JuRDThttps://zestedesavoir.com/forums/sujet/16323/traiter-des-donnees-ordinales-en-fonction-du-temps-dans-r/?page=1#p243698<p>Bonjour,</p>
<p>Dans le cadre de mon stage j’essaie de traiter les données ordinales en fonction du temps afin d’analyser s’il existe une différence significative entre deux modalités différentes.</p>
<p>Voici mon ensemble de données</p>
<div class="table-wrapper"><table>
<thead>
<tr>
<th align="center">Nombre</th>
<th align="center">Note</th>
<th align="center">Semaine</th>
<th align="center">Compagnon</th>
<th align="center">Momies</th>
<th align="center">Mycoses</th>
</tr>
</thead>
<tbody>
<tr>
<td align="center">1</td>
<td align="center">2</td>
<td align="center">13</td>
<td align="center">0</td>
<td align="center">1</td>
<td align="center">2</td>
</tr>
<tr>
<td align="center">2</td>
<td align="center">1</td>
<td align="center">13</td>
<td align="center">0</td>
<td align="center">1</td>
<td align="center">1</td>
</tr>
<tr>
<td align="center">3</td>
<td align="center">4</td>
<td align="center">13</td>
<td align="center">0</td>
<td align="center">1</td>
<td align="center">1</td>
</tr>
<tr>
<td align="center">4</td>
<td align="center">2</td>
<td align="center">13</td>
<td align="center">0</td>
<td align="center">1</td>
<td align="center">1</td>
</tr>
<tr>
<td align="center">31</td>
<td align="center">5</td>
<td align="center">13</td>
<td align="center">1</td>
<td align="center">2</td>
<td align="center">1</td>
</tr>
<tr>
<td align="center">1</td>
<td align="center">3</td>
<td align="center">14</td>
<td align="center">0</td>
<td align="center">3</td>
<td align="center">2</td>
</tr>
</tbody>
</table></div>
<p>num est le numéro de la plante, il y a 30 plantes dans chaque modalité, donc un total de 60 plantes.</p>
<p>Note est la note donnée en fonction de la quantité de pucerons trouvés sur la plante :
1 "absence"
2 "fondatrice
3 "fondatrice + larves
4 "colonie
5 "colonie ailée".</p>
<p>Semaine est le niveau de faiblesse dans lequel les déclarations ont été faites.</p>
<p>Companion est la modalité
0 "sans plante compagne
1 "avec plante compagne".</p>
<p>momies est la note donnée en fonction de la quantité de momies trouvées sur la plante :
1 "absence"
2 "peu"
3 "beaucoup".</p>
<p>J’ai effectué une régression linéaire en utilisant la fonction polr mais le temps est pris en compte ici comme un facteur fixe.</p>
<div class="hljs-code-div hljs-code-r"><div class="hljs-line-numbers"><span data-count="1"></span><span data-count="2"></span><span data-count="3"></span><span data-count="4"></span><span data-count="5"></span><span data-count="6"></span></div><pre><code class="hljs language-r">head(puceron, <span class="hljs-number">2</span>)
aphid$Note = factor(aphid$Note)
e = polr(Note ~ num + Semaine + Compagnon, data = aphid, Hess = <span class="hljs-literal">TRUE</span>)
n4 = data.frame(num = <span class="hljs-number">30</span>, Semaine = <span class="hljs-number">14</span>, Compagnon = <span class="hljs-number">0</span>)
predict(e, n4, type = <span class="hljs-string">"probs"</span>)
<span class="hljs-built_in">exp</span>(cbind(OR = coef(e), confint(e)))
</code></pre></div>
<p>Je veux trouver un moyen de voir s’il existe une différence statistique entre le nombre de pucerons, de momies, de mycoses et de semaines des deux modalités. Existe-t-il un modèle qui puisse faire cela ? </p>