Simulate copy number data for a case-control study.

1 | ```
simul(n, nbSNP, probCas, nbSeg, meanSegmentSize, prob, alpha = 15)
``` |

`n` |
Number of individuals. |

`nbSNP` |
Size of the DNA sequence. |

`probCas` |
Probability to be a case individual. |

`nbSeg` |
Number of causal segments. |

`meanSegmentSize` |
The mean size of anormal segment. |

`prob` |
A 2*2 matrix containing probabilities: prob[1,1]=probability to have an anomaly to a SNP given the person does not have the disease and the SNP is causal. prob[1,2]=probability to have an anomaly to a SNP given the person does not have the disease and the SNP is not causal. prob[2,1]=probability to have an anomaly to a SNP given the person has the disease and the SNP is causal. prob[2,2]=probability to have an anomaly to a SNP given the person has the disease and the SNP is not causal. |

`alpha` |
Parameter of the beta(alpha,alpha). |

a list containing:

`data` |
A matrix of size n*nbSeg, containing values of the copy-number signal. |

`response` |
A vector of size n containing the cas/control status. |

`causalSNP` |
A vector of size nbSeg containing the center of causal segments. |

Quentin Grimonprez, Serge Iovleff

1 |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.