cghseg (version 1.0.5)

simulprofiles: simulprofiles

Description

Simulates multiple CGH profiles

Usage

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)

References

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

Examples

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