我正在运行manova，我的代码是这样的：
fit <- manova(y ~ grp + cov)
summary(fit)
所以汇总的输出是：
Df Pillai approx F num Df den Df Pr(>F)
grp 2 0.185330 5.6511 6 332 1.322e-05 ***
age 1 0.110497 6.8323 3 165 0.0002284 ***
fd 1 0.153049 9.9388 3 165 4.646e-06 ***
scan 1 0.037374 2.1354 3 165 0.0977272 .
Residuals 167
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
我无法找到提取第1列第6个元素（p值）的方法，这样我就可以创建一个条件if循环来对重要模型运行单变量和事后测试。希望有类似的东西：
p <- [method of obtaining p value(1.322e-05) in this case]
if (p < 0.017){
for (i in 1:length(dependentvars){
univariate tests with aov()
p.uni <- univariatetestp(same thing as before with obtaining p value)
if (p.uni <-){
summary(glht(....
非常感谢任何建议或指示。 谢谢 -J