我想生成一个非奇异的数据矩阵,并避免在R中出现以下错误:
"Error in solve.default(sigma) :
system is computationally singular: reciprocal condition number = 1.3971e-16"
我在R中使用is.singular.matrix函数。如何在ifelse函数中编写“ discard x”部分的代码?
n=100
p=50
x<-0.5*matrix(rnorm(n*p),nrow=n,ncol=p)
ifelse(is.singular.matrix(x)=T, discard x, x=x)