Question

mydat=iris

#split sample

index <- sample(1:nrow(mydat),round(0.70*nrow(mydat)))
train <- mydat[index,]
test <- mydat[-index,]
str(train)
z <- lda(Species ~ .,data = train)
z

我该如何计算组质心，以便输出为

IRISTYPE    Function    
          1       2
set -   7,608   ,215
vers    1,825   -,728
virg    5,783   ,513

所以我需要获得组质心的此功能。怎么做？

Answer 1

这是您要找的吗？

install.packages("candisc")
library(candisc) 

iris.mod <- lm(cbind(Petal.Length, Sepal.Length, Petal.Width, Sepal.Width) ~ Species, data=iris)

iris.can <- candisc(iris.mod, term="Species")

iris.can$means

如果您输入iris.can $ means，它将显示这些值

使用R计算组质心

1 个答案: