1 getS3method("print","Anova.mlm")
   2 
   3 get.p.vals <- function(x){
   4     test <- x$test
   5     repeated <- x$repeated
   6     ntests <- length(x$terms)
   7     tests <- matrix(NA, ntests, 4)
   8 
   9     if (!repeated) 
  10         SSPE.qr <- qr(x$SSPE)
  11     
  12     for (term in 1:ntests) {
  13         eigs <- Re(eigen(qr.coef(if (repeated) qr(x$SSPE[[term]]) else SSPE.qr, 
  14             x$SSP[[term]]), symmetric = FALSE)$values)
  15         tests[term, 1:4] <- switch(test, Pillai = stats:::Pillai(eigs, 
  16             x$df[term], x$error.df), Wilks = stats:::Wilks(eigs, x$df[term], 
  17             x$error.df), `Hotelling-Lawley` = stats:::HL(eigs, x$df[term], 
  18             x$error.df), Roy = stats:::Roy(eigs, x$df[term], x$error.df))
  19     }
  20     ok <- tests[, 2] >= 0 & tests[, 3] > 0 & tests[, 4] > 0
  21     ok <- !is.na(ok) & ok
  22     tests <- cbind(x$df, tests, pf(tests[ok, 2], tests[ok, 3], 
  23         tests[ok, 4], lower.tail = FALSE))
  24     rownames(tests) <- x$terms
  25     colnames(tests) <- c("Df", "test stat", "approx F", "num Df", 
  26         "den Df", "Pr(>F)")
  27     return(tests)
  28 }
  29 
  30 
  31 
  32 get.p.vals.all <- function(x){
  33     test <- x$test
  34     repeated <- x$repeated
  35     ntests <- length(x$terms)
  36     tests <- matrix(NA, 1, 4)
  37 
  38     if (!repeated) 
  39         SSPE.qr <- qr(x$SSPE)
  40     
  41     for (term in 1:ntests) {
  42         eigs <- Re(eigen(qr.coef(if (repeated) qr(x$SSPE[[term]]) else SSPE.qr, 
  43             x$SSP[[term]]), symmetric = FALSE)$values)
  44         tests <- rbind(tests,rbind(stats:::Pillai(eigs, x$df[term], x$error.df),
  45                                    stats:::Wilks(eigs, x$df[term], x$error.df),
  46                                    stats:::HL(eigs, x$df[term], x$error.df),
  47                                    stats:::Roy(eigs, x$df[term], x$error.df)))
  48     }
  49 
  50     tests <- tests[-1,]
  51     ok <- tests[, 2] >= 0 & tests[, 3] > 0 & tests[, 4] > 0
  52     ok <- !is.na(ok) & ok
  53     tests <- cbind(rep(x$df,each=4), tests, pf(tests[ok, 2], tests[ok, 3], 
  54         tests[ok, 4], lower.tail = FALSE))
  55     rownames(tests) <- paste(rep(x$terms,each = 4),c("Pillai","Wilks","HL","Roy"))
  56     colnames(tests) <- c("Df", "test stat", "approx F", "num Df", 
  57         "den Df", "Pr(>F)")
  58     return(tests)
  59 }
  60 
  61 #############################################################################
  62 sentences <- read.csv("sentences.csv")
  63 
  64 cols <- names(sentences)[grep("PDWnocc_|PDWhkl_",names(sentences))]
  65 xx <- reshape(sentences[,cols],varying = cols, sep="_",direction = "long")
  66