Подтвердить что ты не робот

Как использовать "hclust" как вызов функции в R

Я попытался построить метод кластеризации как функцию следующими способами:

mydata <- mtcars

# Here I construct hclust as a function
hclustfunc <- function(x) hclust(as.matrix(x),method="complete")

# Define distance metric
distfunc <- function(x) as.dist((1-cor(t(x)))/2)

# Obtain distance
d <- distfunc(mydata)

# Call that hclust function
fit<-hclustfunc(d)

# Later I'd do
# plot(fit)

Но почему это дает следующую ошибку:

Error in if (is.na(n) || n > 65536L) stop("size cannot be NA nor exceed 65536") : 
  missing value where TRUE/FALSE needed

Какой правильный способ сделать это?

4b9b3361

Ответ 1

Прочитайте справку по функциям, которые вы используете. ?hclust довольно ясно, что первый аргумент d является объектом несходства, а не матрицей:

Arguments:

       d: a dissimilarity structure as produced by ‘dist’.

Update

Поскольку OP теперь обновил свой вопрос, в чем нуждаются

hclustfunc <- function(x) hclust(x, method="complete")
distfunc <- function(x) as.dist((1-cor(t(x)))/2)
d <- distfunc(mydata)
fit <- hclustfunc(d)

Оригинальные

Вы хотите

hclustfunc <- function(x, method = "complete", dmeth = "euclidean") {    
    hclust(dist(x, method = dmeth), method = method)
}

а затем

fit <- hclustfunc(mydata)

работает как ожидалось. Обратите внимание, что теперь вы можете передать метод коэффициента несходства как dmeth и метод кластеризации.