Replies: 4 comments 5 replies
-
Strange thing is that when I use the same commands, I get both a table of p-values and a table of F statistics. Can you prepare a reproducible example, where 'reproducible' means that we can reproduce the example. Using, say, data in vegan is a good choice so that you do not need to provide your own data (and I usually do not want to download alien data). > m <- with(dune.env, betadisper(vegdist(dune), Management))
> permutest(m, pairwise = TRUE)
Permutation test for homogeneity of multivariate dispersions
Permutation: free
Number of permutations: 999
Response: Distances
Df Sum Sq Mean Sq F N.Perm Pr(>F)
Groups 3 0.13831 0.046104 1.9506 999 0.15
Residuals 16 0.37816 0.023635
Pairwise comparisons:
(Observed p-value below diagonal, permuted p-value above diagonal)
BF HF NM SF
BF 0.570000 0.058000 0.267
HF 0.593549 0.076000 0.351
NM 0.060191 0.064936 0.372
SF 0.281235 0.402333 0.365221 |
Beta Was this translation helpful? Give feedback.
-
Hi @jarioksa, What I meant was that I don't get a F statistic for every pairwise comparision. Using the dune data example that you showed, I would like to know the F statistic of HF x BF, HF x NM, HF x SF.... It only gives you the overall F statisitic for all groups. Thanks |
Beta Was this translation helpful? Give feedback.
-
Why? Pairwise tests use t-statistics. As a first trick, you may try F 1,p = tp2. Nothing of this makes sense in permutation tests that do not use theoretical distributions, but are based on ... permutations. I suggest you consult your statistical text book, and prepare to explain to referees why you used F statistics in pairwise tests instead of t statistics. |
Beta Was this translation helpful? Give feedback.
-
If you want them, the observed statistics (F for the omnibus test, t for each pairwise test) are available in the # load data
data(varespec)
# Hellinger distances between samples
dis <- vegdist(varespec, method = "hellinger")
# First 16 sites grazed, remaining 8 sites ungrazed
groups <- factor(c(rep(1,16), rep(2,8)), labels = c("grazed","ungrazed"))
# Calculate multivariate dispersions
mod <- betadisper(dis, groups)
mod
# Perform test
anova(mod)
# Permutation test for F
pmod <- permutest(mod, permutations = 99, pairwise = TRUE)
# extract the observed test statistics
pmod$statistic which gives
|
Beta Was this translation helpful? Give feedback.
-
I am using the commands:
But the resulting table only give p-values for each pairwise comparison, but not F. Is there a way to see what the F value is for all these pairwise comparisions?
Thanks
Beta Was this translation helpful? Give feedback.
All reactions