this is for holding javascript data
David LeBauer edited stat conversions.md
almost 10 years ago
Commit id: 94593a18cf315dc4ffd3dd11bb7d25bf42ed6493
deletions | additions
diff --git a/stat conversions.md b/stat conversions.md
index 2cdee77..05bdee6 100644
--- a/stat conversions.md
+++ b/stat conversions.md
...
**Table 8: How to convert statistics from \(P\), \(LSD\), or \(MSD\) to \(SE\)**
Most of these (and others) are implemented in the [transformstats](https://github.com/PecanProject/pecan/blob/master/utils/R/transformstats.R) function within the PEcAn.utils package
| From | To | Conversion | Rcode R code | Notes |
|:-----|:---|:-----------|:------|:------|
| P | SE | SE\( = \frac{\bar{X}_1-\bar{X}_2}{t_{1-P/2,2n-2}\sqrt{2/n}}\) | (x1-x2)/(qt(1-P/2,2*n-2)*sqrt(2/n)) `(x1-x2)/(qt(1-P/2,2*n-2)*sqrt(2/n))` | \(\bar{X}_{1,2}\) are two means being compared. |
| LSD | SE | \(SE = \frac{LSD}{t_{1-\alpha/2,n}*\sqrt{2b}}\) | LSD/(qt(1-P/2,n)*sqrt(2*b)) `LSD/(qt(1-P/2,n)*sqrt(2*b))` | where \(b\) is the number of blocks, \(n\) is the number of replicates, and \(n=b\) in a Randomized Complete Block Design |
| MSD | SE | \(SE = \frac{MSD*n}{t_{1-\alpha, 2n-2}*\sqrt{2}}\) | msd*n/(qt(1-P/2,2*n-2)*sqrt(2)) `msd*n/(qt(1-P/2,2*n-2)*sqrt(2))` | |
See related questions on Stats.SE: http://stats.stackexchange.com/q/2917/1381 and http://stats.stackexchange.com/q/4485/1381