mirror of
https://github.com/Polarolouis/anova-phylogenetique-projet-msv.git
synced 2026-06-17 10:15:25 +02:00
🐛 Add todo + add norm definition
This commit is contained in:
parent
8c56ad8216
commit
ba810efe8c
2 changed files with 11 additions and 8 deletions
19
rapport.Rnw
19
rapport.Rnw
|
|
@ -235,14 +235,17 @@ On testera alors les hypothèses suivantes :
|
||||||
\begin{equation}
|
\begin{equation}
|
||||||
F=\frac{||\hat{Y} - \bar{Y}||^2_{V_\lambda^{-1}}(n-2)}{||Y - \hat{Y}||^2_{V_\lambda^{-1}}} \underset{\mathcal{H}_0}{\sim}\mathcal{F}\text{isher} (1, n-2)
|
F=\frac{||\hat{Y} - \bar{Y}||^2_{V_\lambda^{-1}}(n-2)}{||Y - \hat{Y}||^2_{V_\lambda^{-1}}} \underset{\mathcal{H}_0}{\sim}\mathcal{F}\text{isher} (1, n-2)
|
||||||
\end{equation}
|
\end{equation}
|
||||||
|
\begin{align*}
|
||||||
|
&\text{Où }||Y - \hat{Y}||^2_{V_\lambda^{-1}} = ||Y - X\hat{\beta}||^2_{V_\lambda^{-1}}= Proj_X^{V_\lambda}Y= (Y-\hat{Y})^TV^{-1}_\lambda(Y-\hat{Y})\\
|
||||||
|
&\text{et }||\hat{Y} - \bar{Y}||^2_{V_\lambda^{-1}}=(\hat{Y} - \bar{Y})^TV^{-1}_\lambda(\hat{Y}- \bar{Y})
|
||||||
|
\end{align*}
|
||||||
|
|
||||||
|
|
||||||
TODO: Rajouter def norme
|
|
||||||
|
|
||||||
\subsection{Approximation de Satterthwaite}
|
\subsection{Approximation de Satterthwaite}
|
||||||
|
|
||||||
On va dans notre cas avoir $n-2$ degrés de liberté.
|
On va dans notre cas avoir $n-2$ degrés de liberté.
|
||||||
L'ANOVA, suppose souvent une homoscédasticité des variances entre les groupes ou les échantillons. Cela signifie que les variances des groupes sont égales.
|
L'ANOVA, suppose souvent une homoscédasticité des variances entre les groupes ou les échantillons. Cela signifie que les variances des groupes sont égales.
|
||||||
|
|
||||||
Cependant, lorsque cette condition n'est pas satisfaite, l'approximation de Satterthwaite peut être utilisée pour tenir compte des variances inégales entre les groupes. Elle est particulièrement utile dans le cas des ANOVA à un facteur, mais peut également être appliquée à des ANOVA à plusieurs facteurs.
|
Cependant, lorsque cette condition n'est pas satisfaite, l'approximation de Satterthwaite peut être utilisée pour tenir compte des variances inégales entre les groupes. Elle est particulièrement utile dans le cas des ANOVA à un facteur, mais peut également être appliquée à des ANOVA à plusieurs facteurs.
|
||||||
|
|
||||||
L'approximation de Satterthwaite ajuste les degrés de liberté pour tenir compte de ces différences dans les variances.
|
L'approximation de Satterthwaite ajuste les degrés de liberté pour tenir compte de ces différences dans les variances.
|
||||||
|
|
@ -250,16 +253,19 @@ L'approximation de Satterthwaite ajuste les degrés de liberté pour tenir compt
|
||||||
Cela permet d'obtenir des résultats plus fiables lorsque les conditions d'homoscédasticité ne sont pas respectées.
|
Cela permet d'obtenir des résultats plus fiables lorsque les conditions d'homoscédasticité ne sont pas respectées.
|
||||||
\newline
|
\newline
|
||||||
\newline
|
\newline
|
||||||
On s'est basé sur la documentation du package lmerTest \cite{kuznetsovaLmerTestPackageTests2017} pour calculer les formules explicites de l'approximation dans notre cadre et ensuite les implémenter et voir si cela améliore la fiabilité de la statistique de test.
|
On s'est basé sur la documentation du package lmerTest \cite{kuznetsovaLmerTestPackageTests2017} pour calculer les formules explicites de l'approximation dans notre cadre.
|
||||||
|
En effet il existe des formules explicite dans le cadre du modèle mixte.
|
||||||
|
Cela nous permettra ensuite de les implémenter et voir si cela améliore la fiabilité de la statistique de test.
|
||||||
A partir de \ref{eq:eq2err} on rappelle les valeurs suivantes:
|
A partir de \ref{eq:eq2err} on rappelle les valeurs suivantes:
|
||||||
\[
|
\[
|
||||||
Y \sim \mathcal{N}_n(X\beta, \sigma^2_{phy}K + \sigma^2_{err}I_n) \text{, }\theta=(\sigma^2_{phy}, \sigma^2_{err}) \text{ et } Var_\theta(Y) = V(\theta) = \sigma^2_{phy}K + \sigma^2_{err}I_n
|
Y \sim \mathcal{N}_n(X\beta, \sigma^2_{phy}K + \sigma^2_{err}I_n) \text{, }\theta=(\sigma^2_{phy}, \sigma^2_{err}) \text{ et } Var_\theta(Y) = V(\theta) = \sigma^2_{phy}K + \sigma^2_{err}I_n
|
||||||
\]
|
\]
|
||||||
De là on obtient:
|
|
||||||
|
De la documentation on obtient alors la covariance suivante:
|
||||||
\begin{equation}
|
\begin{equation}
|
||||||
C(\theta) = (Cov(\beta_i , \beta_j))_{i,j} = (X^TV(\theta)^{-1}X)^{-1} = (X^T(\sigma^2_{phy}K + \sigma^2_{err}I_n)^{-1}X)^{-1}
|
C(\theta) = (Cov(\beta_i , \beta_j))_{i,j} = (X^TV(\theta)^{-1}X)^{-1} = (X^T(\sigma^2_{phy}K + \sigma^2_{err}I_n)^{-1}X)^{-1}
|
||||||
\end{equation}
|
\end{equation}
|
||||||
|
TODO: Préciser df degré de liberté de quoi !
|
||||||
Toujours en suivant la documentation \cite{kuznetsovaLmerTestPackageTests2017} on part de l'expression pour les degrés de liberté $df$ et de l'approximation. Ce qui nous donne :
|
Toujours en suivant la documentation \cite{kuznetsovaLmerTestPackageTests2017} on part de l'expression pour les degrés de liberté $df$ et de l'approximation. Ce qui nous donne :
|
||||||
\begin{equation}
|
\begin{equation}
|
||||||
df = \frac{2(l^T\hat{C}l)^2}{[Var(l^T\hat{C}l)]}=\frac{2(f(\hat{\theta}))^2}{[Var(f(\hat{\theta}))]}\approx \frac{2(f(\hat{\theta}))^2}{[\nabla f(\hat{\theta})]^T A[\nabla f(\hat{\theta})]}
|
df = \frac{2(l^T\hat{C}l)^2}{[Var(l^T\hat{C}l)]}=\frac{2(f(\hat{\theta}))^2}{[Var(f(\hat{\theta}))]}\approx \frac{2(f(\hat{\theta}))^2}{[\nabla f(\hat{\theta})]^T A[\nabla f(\hat{\theta})]}
|
||||||
|
|
@ -339,9 +345,6 @@ De là on obtient la Hessienne $\begin{pmatrix}
|
||||||
\end{pmatrix}$ puis A en l'inversant, ce qui peut se faire par des méthodes numériques.
|
\end{pmatrix}$ puis A en l'inversant, ce qui peut se faire par des méthodes numériques.
|
||||||
|
|
||||||
\end{proof}
|
\end{proof}
|
||||||
TODO: trouver erreur de signe
|
|
||||||
% \[\]
|
|
||||||
% Supposons que nous avons une expression $x^2 - 2x + \cancel{3} - 3$. Comme la partie $\cancel{3}$ est nulle, nous pouvons la barrer.
|
|
||||||
|
|
||||||
% TODO REML voir sujet d'exam corrigée
|
% TODO REML voir sujet d'exam corrigée
|
||||||
% Quand estimateur classique on divise par n-p au lieu de diviser par n donc on
|
% Quand estimateur classique on divise par n-p au lieu de diviser par n donc on
|
||||||
|
|
|
||||||
BIN
rapport.pdf
BIN
rapport.pdf
Binary file not shown.
Loading…
Add table
Reference in a new issue