Commit c6ffa3b4 by Eric Coissac

Quelques notes...

parent c8cb261e
...@@ -32,6 +32,7 @@ knitr::opts_knit$set( ...@@ -32,6 +32,7 @@ knitr::opts_knit$set(
- connecter avec glasso - connecter avec glasso
- connecter avec path analysis (lavaan) - connecter avec path analysis (lavaan)
- test randomisation anova - test randomisation anova
- regarder les functions : shuffle, how
- test rendomisation coeficients - test rendomisation coeficients
## interrogation : ## interrogation :
...@@ -518,6 +519,7 @@ data = procmod.frame(euk=euk.pco.li,bac=bac.pco.li, ...@@ -518,6 +519,7 @@ data = procmod.frame(euk=euk.pco.li,bac=bac.pco.li,
geo=geo.pco.li,hist=dh.pco.li) geo=geo.pco.li,hist=dh.pco.li)
vars = mvar(data) vars = mvar(data)
``` ```
```{r echo=FALSE} ```{r echo=FALSE}
...@@ -551,7 +553,7 @@ plot(euk.pm) ...@@ -551,7 +553,7 @@ plot(euk.pm)
```{r} ```{r}
W=1/rowSums(euk.pm$residuals^2) W=1/rowSums(euk.pm$residuals^2)
W=W/max(W) W=W/max(W)
euk.pm.w = pm(euk ~ soil + climat + geo + hist ,data=data,weights = W) euk.pm.w = pm(euk ~ soil + climat + geo + hist + climat:hist,data=data,weights = W)
euk.pm.w euk.pm.w
``` ```
...@@ -576,6 +578,7 @@ partition ...@@ -576,6 +578,7 @@ partition
``` ```
```{r} ```{r}
euk.anova.w = anova(euk.pm.w) euk.anova.w = anova(euk.pm.w)
euk.anova.w euk.anova.w
``` ```
```{r} ```{r}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment