Find centralized, trusted content and collaborate around the technologies you use most.
Teams
Q&A for work
Connect and share knowledge within a single location that is structured and easy to search.
Rには対称であるはずのマトリックスがありますが、マシンの精度により、マトリックスは対称ではありません(値は約10 ^ -16異なります)。マトリックスが対称であることを知っているので、これまでのところ、問題を回避するためにこれを行ってきました。
s.diag = diag(s) s[lower.tri(s,diag=T)] = 0 s = s + t(s) + diag(s.diag,S)
これのためのより良い1行のコマンドはありますか?
s<-matrix(1:25,5) s[lower.tri(s)] = t(s)[lower.tri(s)]