Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/77.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
使用car软件包获得的MANOVA对象的xtable_R_Xtable_R Car_Manova - Fatal编程技术网

使用car软件包获得的MANOVA对象的xtable

使用car软件包获得的MANOVA对象的xtable,r,xtable,r-car,manova,R,Xtable,R Car,Manova,我想知道如何获得xtable的MANOVA对象获得汽车包。这里是MWE: library(xtable) library(car) MANOVA <- Anova(lm(cbind(Al, Fe, Mg, Ca, Na) ~ Site, data=Pottery)) xtable(MANOVA) 库(xtable) 图书馆(汽车) MANOVA问题在于xtable对类anova的对象有一个方法,但是car包的anova函数返回类anova.mlm的对象。我在以下位置找到了修复程序: 库(

我想知道如何获得xtable的MANOVA对象获得汽车包。这里是MWE:

library(xtable)
library(car)
MANOVA <- Anova(lm(cbind(Al, Fe, Mg, Ca, Na) ~ Site, data=Pottery))
xtable(MANOVA)
库(xtable)
图书馆(汽车)

MANOVA问题在于
xtable
对类
anova
的对象有一个方法,但是
car
包的
anova
函数返回类
anova.mlm
的对象。我在以下位置找到了修复程序:

库(xtable)
图书馆(汽车)
#创建一些示例数据

Pottery问题在于
xtable
对类
anova
的对象有一个方法,但是
car
包的
anova
函数返回类
anova.mlm
的对象。我在以下位置找到了修复程序:

库(xtable)
图书馆(汽车)
#创建一些示例数据

我们能想出一个比[汽车]更好的标签名吗?这只是过于通用的一面上的一点,而且会被可怕地滥用……我们能想出一个比[汽车]更好的标签名吗?这只是过度泛化的一点,将会被可怕地误用……这很好,但它删掉了我添加的标题和标签。知道怎么把那些拿回来吗?没关系。看起来像是对函数进行了最后一次xtable调用:
xtable(tests,…)
使它按照我希望的方式工作。这很好,但它删除了我添加的标题和标签。知道怎么把那些拿回来吗?没关系。看起来像是对函数进行了最后一次xtable调用:
xtable(tests,…)
使它按照我希望的方式工作。
library(xtable)
library(car)

# Create some example data
Pottery <- data.frame(
  "Al" = rnorm(10),
  "Fe" = rnorm(10),
  "Mg" = rnorm(10),
  "Ca" = rnorm(10),
  "Na" = rnorm(10),
  "Site" = sample(LETTERS[1:2], 10, replace = TRUE))

# Create a custom function handling `anova.mlm` objects
xtable.Anova.mlm <- function (x, ...) {
  test <- x$test
  repeated <- x$repeated
  ntests <- length(x$terms)
  tests <- matrix(NA, ntests, 4)
  if (!repeated)
    SSPE.qr <- qr(x$SSPE)
  for (term in 1:ntests) {
    eigs <- Re(eigen(qr.coef(if (repeated) qr(x$SSPE[[term]]) else
      SSPE.qr,
      x$SSP[[term]]), symmetric = FALSE)$values)
    tests[term, 1:4] <- switch(test, Pillai = stats:::Pillai(eigs,
      x$df[term], x$error.df), Wilks = stats:::Wilks(eigs,
        x$df[term], x$error.df), `Hotelling-Lawley` = stats:::HL(eigs,
          x$df[term], x$error.df), Roy = stats:::Roy(eigs,
            x$df[term], x$error.df))
  }
  ok <- tests[, 2] >= 0 & tests[, 3] > 0 & tests[, 4] > 0
  ok <- !is.na(ok) & ok
  tests <- cbind(x$df, tests, pf(tests[ok, 2], tests[ok, 3],
    tests[ok, 4], lower.tail = FALSE))
  rownames(tests) <- x$terms
  colnames(tests) <- c("Df", "test stat", "approx F", "num Df",
    "den Df", "Pr(>F)")
  tests <- structure(as.data.frame(tests), heading = paste("\nType ",
    x$type, if (repeated)
      " Repeated Measures", " MANOVA Tests: ", test, " test
statistic",
    sep = ""), class = c("anova", "data.frame"))
  #    print(tests)
  #    invisible(x)
  xtable(tests)
}


MANOVA <- Anova(lm(cbind(Al, Fe, Mg, Ca, Na) ~ Site, data=Pottery))

xtable(MANOVA)

% latex table generated in R 2.15.2 by xtable 1.7-1 package
% Tue Jun 11 20:49:51 2013
\begin{table}[ht]
\centering
\begin{tabular}{lrrrrrr}
\hline
& Df & test stat & approx F & num Df & den Df & Pr($>$F) \\ 
\hline
Site & 1 & 0.61 & 1.24 & 5 & 4 & 0.4288 \\ 
\hline
\end{tabular}
\end{table}