Description Usage Arguments Value Examples

Calls the robnb function to estimate the coefficients, and then construct the statistical tests for DD and DE. It works for a single gene. y1 and y2 are count vectors for a single gene. diphiseq calls this function to do the calculation for each gene. Normal users often don't need to use this function directly.

1 | ```
robtest(y1, log.depth1, y2, log.depth2, c.tukey.beta = 4, c.tukey.phi = 4)
``` |

`y1` |
counts from group 1. a vector. |

`log.depth1` |
log(sequencing depths) for samples in group 1. a vector. |

`y2` |
counts from group 2. a vector. |

`log.depth2` |
log(sequencing depths) for samples in group 2. a vector. |

`c.tukey.beta` |
The c value for beta in Huber function. The default value, 4, is typically regarded as appropriate and should work for most datasets. |

`c.tukey.phi` |
The c value for phi in Huber function. The default value, 4, is typically regarded as appropriate and should work for most datasets. |

A vector that contains the elements:
`phi1`

: the estimated dispersion of sample group 1.
`phi2`

: the estimated dispersion of sample group 2.
`beta1`

: the estimated (log) expression of sample group 1.
`beta2`

: the estimated (log) expression of sample group 2.
`statistic.phi`

: the z statistic for DD.
`statistic.beta`

: the z statistic for DE.
`p.value.phi`

: the p value for DD.
`p.value.beta`

: the p value for DE.

1 2 3 4 5 |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.