如何绘制多元对应分析的二元（二元）预测变量的置信椭圆？

Question

我试图用分类预测器绘制MCA中的95％置信度elipses时遇到了一些麻烦。以下是我的数据和代码示例。

我不明白问题出在哪里来自和我在（M）CA新：

CONSTRUCTION<-c("in plaats van", "ten opzichte van", "met behulp van", "ten koste van", "op grond van", "onder leiding van", "in tegenstelling tot", "op basis van") 
IDIOSYNC<-c("no", "yes", "yes", "yes", "no", "no", "no", "no") 
ARTICLE<-c("yes", "no", "no", "no", "no", "no", "no", "no") 
PLURAL<-c("no", "no", "no", "no", "no", "no", "no", "no") 
MODIF<-c("no", "no", "no", "no", "no", "yes", "yes", "no") 
EXTRAPOS<-c("yes", "no", "no", "yes", "no", "yes", "yes", "no") 
COMPLEM<-c("no", "no", "no", "no", "no", "yes", "no", "no") 
P2OPTION<-c("no", "no", "no", "no", "no", "yes", "no", "no") 
P1CHBL<-c("no", "no", "no", "no", "no", "no", "no", "no") 
P2CHBL<-c("no", "no", "no", "no", "no", "no", "yes", "no") 

data.comp.prep<-cbind(CONSTRUCTION, IDIOSYNC, ARTICLE, PLURAL, MODIF, EXTRAPOS, COMPLEM, P2OPTION, P1CHBL, P2CHBL) 

ca.comp.prep<-MCA(data.comp.prep[, -c(1, 4, 9)], graph=FALSE) 
plot(ca.comp.prep, cex=0.7, col.var="black", col.ind="grey") 
ca.comp.prep1<-MCA(data.comp.prep[, -c(4, 9)], quali.sup=1, graph=FALSE) 
plot(ca.comp.prep1, invis="ind", col.var="darkgrey", col.quali.sup="black") 

的问题，当我为绘制椭圆运行代码出现。

Answer 1

ca.comp.prep1是包含MCA结果的列表。你试图索引它，就好像它是一个有两个维度的data.frame。这里的[,c(4,9)]是不合情理的。

然后，您试图为第一个变量ca.comp.prep1绘制省略号，它基本上是一个行ID（它对每个人都有不同的值）。没有椭圆是可能的，因为每个类别只有一个点。

我的建议是，你不应该在MCA中添加第一个变量;继续学习ca.comp.prep，忘记ca.comp.prep1。然后绘制椭圆：

plotellipses(ca.comp.prep, keepvar=1, label="quali") 

如果您仍想上ca.comp.prep1完成后，尝试另一个变量：

plotellipses(ca.comp.prep1, keepvar=2, label="quali")