- In order to install this package:
library(devtools)
install_github("elnazmirzaei/FIESTA")
library(FIESTA)
library(NMF)
library(mixtools)
library(fitdistrplus)
- After installation you need to prepare your gene-expression dataset as a matrix in R such that rows represent genes and columns represent cells. Then, pick whether to use WNMF or sNMF. Default is sNMF. Based on the size of your data this step might take a while.
A = GeneExpressionData
A_imputed = impute( A , method="sNMF" )
- Now it is time for scaling.
A_imputed_scaled = scale( A , A_imputed )
- And at the end run the thresholding step.
A_imputed_scaled_thresholded = thresholding ( A , A_imputed_scaled )