Cellchat log fold change

Question

[miyeok]

Hi,

I'm doing cell-cell interaction anlaysis with Cellchat package in R. To identify over-expressed signaling genes, I want to set the thresh.fc argument of OverExpressedGenes() function.

I have a question about this thresh.fc argument.

Documentation says it's Log fold change threshold. those that mean Log2 fold-change? or Log10 fold-change?

Thank you

Answer 1

[Nitin Narwade]

It seems they are taking the natural log of the mean expression and calculating the difference.

mean.fxn <- function(x) {
    return(log(x = mean(x = expm1(x = x)) + 1)) 
}

data.1 <- apply(X = data.use[features, cell.use1, drop = FALSE],MARGIN = 1,FUN = mean.fxn)
data.2 <- apply(X = data.use[features, cell.use2, drop = FALSE],MARGIN = 1,FUN = mean.fxn)

FC <- (data.1 - data.2)

You can have a look at the code here.

[miyeok]

Thank you so much!!