Error in "Modeling the technical noise in gene expression" with trendVar()

[marionilab]

Error in "Modeling the technical noise in gene expression" with trendVar()

[vyturing]

Just fixed it! When my local (Mac Majove 10.14.16) was turned off automatically due to low power and I charged the laptop to turn it on again, the RStudio session which was initially open with my half complete workflow became somewhat intrinsically corrupted though everything looked perfectly same as before. It started acting weird. I rebooted the computer and completed the entire workflow without hassle. That is good to know. The laptop does not seem to recover 100% from a low battery shutdown! I think this is worth publishing especially for people like me looking all over for hours to find the bug in a perfectly working pipeline.

[LTLA]

The error suggests that a seq_len is encountering a NULL somewhere. If you are following the code exactly, I don't know how this could be possible. But in any case, I suspect you're looking at an old version of the workflows that don't match up to the latest version of the package. The workflows were absorbed into the OSCA book, note the redirection mentioned here; the updated call is to modelGeneVar(), and you should have gotten a deprecation message to that effect.

[LTLA]

Ah, the good old "turn it off and on again"... literally!

I would doubt that any memory state remains valid after a hard shutdown like that, so it's best to assume that it is not... save frequently and reboot from the last saved state.

[vyturing]

Exactly! Thank you so very much again and again..