simulprofiles

Share:

Description

Simulates multiple CGH profiles

Usage

1
simulprofiles(M, n,k.mean,SNR,lambda)

Arguments

M

Number of profiles

n

length of profiles

k.mean

average number of segments per profile

SNR

Signal to Noise Ratio between the mean and the measurement noise. The mean of segments can take 3 values -1;0;0.58.

lambda

Variance ratio between measurement noise and bakground noise

Value

simul

list containing Y (dataframe of simulated CGH signals), Ki (the true number of segments per profile), mu0 (the true signal), theta0 (the true background intensity)

Author(s)

F. Picard, E. Lebarbier, M. Hoebeke, G. Rigaill, B. Thiam, S. Robin

References

Joint segmentation, calling and normalization of multiple CGH profiles, Biostatistics (2011)

Examples

1
2
3
#set.seed(1)
#simul          = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1)
#CGHd           = new("CGHdata",Y=simul$Y)