what is residualModelFormulaStr in monocle?
2
1
Entering edit mode
3.3 years ago
John ▴ 220

Hi, I want to Cluster cells using marker genes, I use monocle, I get the following error, Here I don't know what to replace with ~Media, can anyone please explain what is residualModelFormulaStr ?what should I give there?

> marker_diff <- markerDiffTable(HSMM[expressed_genes,], cth, residualModelFormulaStr = "~Media + num_genes_expressed", cores = 1)
Error in `[.data.frame`(pd, , i) : undefined columns selected
r R RNA-Seq software error monocle • 2.5k views
ADD COMMENT
0
Entering edit mode
3.1 years ago
Batu ▴ 190

Hello, I also could not find what it is exactly. But I solved it by trial and error, it is solved when i delete Media, leave num_genes_expressed only.

marker_diff <- markerDiffTable(pfcm[expressed_genes,], cth, residualModelFormulaStr = "~num_genes_expressed", cores = 2)

But I also need to learn what it is, since it is used in differentialGeneTest function.

ADD COMMENT
1
Entering edit mode

I think you're right. Here is how they set up the experiment: To initiate myoblast differentiation, we switch media from a high-mitogen growth medium (GM) to a low-mitogen differentiation medium (DM). So they need to exclude the effects of "culture media" from clustering. If you don't have such effects that potentially cause bias then there is no need to add it.

ADD REPLY
0
Entering edit mode

I think you're right. Here is how they set up the experiment: To initiate myoblast differentiation, we switch media from a high-mitogen growth medium (GM) to a low-mitogen differentiation medium (DM). So they need to exclude the effects of "culture media" from clustering. If you don't have such effects that potentially cause bias then there is no need to add it.

ADD REPLY
0
Entering edit mode
3.1 years ago
Batu ▴ 190

In this issue in GitHub, you can find the andwer about formula string: https://github.com/cole-trapnell-lab/monocle-release/issues/151

ADD COMMENT

Login before adding your answer.

Traffic: 2211 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6