Question

我想在xtable输出上使用candisc::cancor函数。但是想不出来。在这方面的任何帮助将受到高度赞赏。感谢

library(xtable)
library(candisc)
data(Rohwer, package="heplots")
X <- as.matrix(Rohwer[,6:10])  # the PA tests
Y <- as.matrix(Rohwer[,3:5])   # the aptitude/ability variables
(cc <- cancor(X, Y, set.names=c("PA", "Ability")))
Canonical correlation analysis of:
         5   PA  variables:  n, s, ns, na, ss 
  with   3   Ability  variables:  SAT, PPVT, Raven 

    CanR  CanRSQ   Eigen percent    cum                          scree
1 0.6703 0.44934 0.81599   77.30  77.30 ******************************
2 0.3837 0.14719 0.17260   16.35  93.65 ******                        
3 0.2506 0.06282 0.06704    6.35 100.00 **                            

Test of H0: The canonical correlations in the 
current row and all that follow are zero

     CanR  WilksL      F df1   df2  p.value
1 0.67033 0.44011 3.8961  15 168.8 0.000006
2 0.38366 0.79923 1.8379   8 124.0 0.076076
3 0.25065 0.93718 1.4078   3  63.0 0.248814
xtable(cc)

Error in UseMethod("xtable") : 
  no applicable method for 'xtable' applied to an object of class "cancor"

Answer 1

如果您查看methods(xtable)，cancor对象没有方法，但我们可以使用xtable.data.frame，因为cancor的输出看起来像{{ 1}}。

但是，data.frames不会返回这些值 - 请查看cancor。它们以str(cc)方法计算。看看print。

因此，创建一个快速提取功能，以candisc:::print.cancor

生成cancor输出

data.frames

运行辅助函数

outfun <- function(X) { 
              CanR = X$cancor  
              CanRSQ = CanR^2  
              Eigen = CanRSQ/(1 - CanRSQ)
              percent = 100 * Eigen/sum(Eigen)
              cum = cumsum(percent)
              list(tab1 = data.frame(CanR, CanRSQ, Eigen, percent, cum) ,
                   tab2 = as.data.frame(candisc:::Wilks.cancor(X)))
        }

然后可以使用(out <- outfun(cc)) # or run it like out <- outfun(cancor(X, Y, set.names=c("PA", "Ability")))生成xtable.data.frame输出

xtable

您可以调整辅助函数以提取标题，然后将其作为参数传递给xtable:::xtable.data.frame(out[[1]]) # % latex table generated in R 3.1.1 by xtable 1.7-3 package # % Wed Sep 17 18:32:07 2014 # \begin{table}[ht] # \centering # \begin{tabular}{rrrrrr} # \hline # & CanR & CanRSQ & Eigen & percent & cum \\ # \hline # 1 & 0.67 & 0.45 & 0.82 & 77.30 & 77.30 \\ # 2 & 0.38 & 0.15 & 0.17 & 16.35 & 93.65 \\ # 3 & 0.25 & 0.06 & 0.07 & 6.35 & 100.00 \\ # \hline # \end{tabular} # \end{table} xtable:::xtable.data.frame(out[[2]])。

可以帮助您入门的快速答案

在candisc :: cancor输出上使用xtable函数

1 个答案: