【发布时间】:2021-01-27 05:18:02
【问题描述】:
下午好!
我正在尝试从头开始实现 EM 算法。我开发的代码包含:
k=10
w_k=rep(1,k)/k
n_j=rep(0,k)
print(w_k)
data=iris[1:150,-5]
means=sample(1:dim(data)[1],k,replace=FALSE)
means
mu=iris[means,-5]
sigma=cov(data)
print("mu")
mu
print("sigma")
sigma
print("inv")
solve(sigma)
print("mu[1,] ")
as.vector(mu[1,])
as.vector(mu[1,]) %*% solve(sigma) #line of error
不幸的是,这一行给出了以下错误:
Error in as.vector(mu[1, ]) %*% solve(sigma) :
requires numeric/complex matrix/vector arguments
Execution halted
我多次尝试使用as.matrix() 解决问题,但均未成功。
提前感谢您的帮助!
【问题讨论】:
标签: r