R/sqrtsingular.R

In Hlch1992/Rprojects: Statistical modeling and inference of microbial interaction and stability for microbial dynamical systems

sqrtsingular=function(sqrtmatrix) {
  Xm = eigen(sqrtmatrix)
  Tm = Xm$vectors
  Jm = Diagonal( x=Xm$values )
  Tinv = solve(Tm)
  Tm %*% Jm %*% Tinv
  Jsqrt = Diagonal( x=sqrt( Xm$values ) )
  Jsqrt[is.na(Jsqrt)]=0
  SDsqrt = Tm %*% Jsqrt %*% Tinv
  return(SDsqrt)
}