powered by
Simulate copy number data for a case-control study.
simul(n, nbSNP, probCas, nbSeg, meanSegmentSize, prob, alpha = 15)
a list containing:
A matrix of size n*nbSeg, containing values of the copy-number signal.
A vector of size n containing the cas/control status.
A vector of size nbSeg containing the center of causal segments.
Number of individuals.
Size of the DNA sequence.
Probability to be a case individual.
Number of causal segments.
The mean size of abnormal segment.
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.
Parameter of the beta(alpha,alpha).
Quentin Grimonprez, Serge Iovleff
data <- simul(50, 10000, 0.4, 10, 150, matrix(c(0.1, 0.8, 0.001, 0.001), nrow = 2))
Run the code above in your browser using DataLab