Conditional Random Pattern Algorithm for LOH Inference and Segmentation

Software Usage:

Step 1: Download the package, and decompress it. For example, programs are decompressed at: c:\CRP_LOH

Step 2: Run Matlab (version 7.0 and higher), and change the current directory to the ‘c:\CRP_LOH’.

Step 3: Test the file: “crpSNPLOHInfer.m”.



There is one major function: crpSNPLOHInfer.m

crpSNPLOHInfer.m implements the LOH inference using the CRP model. There are some parameters are adjustable in simSNP_CRP_Infer.m. Please check the Bioinformatics for more details.


 Input Data:

The input data of the CRP_LOH model is the “log2-ratio” signals (which is obtained as: log2(disease snp intensity/normal snp intensity). In this study, we make use of “CNAG” to extract the log2-ratio sequences from the SNP array files (.CEL and .CHP). The users, who want to test the some new SNP array files, please use the CNAG to extract the log2-ratio sequence first. To download CNAG at:


There are two simulated data sets. One is the simulated log2-ratios. The other is the simulated 500K SNP array data based on the HapMap samples public available on the Affymetrix website.

Simulated Log2-ratio data set: simLog2RatioData.mat (load it into matlab)

Simulated SNP array data set: ‘.\HapMap\ *.mat’.


The HapMap sample NA10851 was used to simulated these data at three different SNR levels. At same SNR level, there are two simulated data. One simulated the copy number loss, the other simulated the copy number gain. The original CEL and CHP files are also provided in ‘HapMapOri’ directory. The users can test them on different software packages that designed for Affymetrix array.


Important NOTE about CRF matlab code developed by Dr. Kevin Murthy:

We used the condition random field codes developed by Dr. Kevin Murthy, which could be downloaded at: . We also clarify this in our paper


This work is published: Wu LY, Zhou X, Li FH, Yang XR, Chang CC and Wong STC. Conditional Random Pattern Algorithm for LOH Inference and Segmentation. Bioinformatics, 25(1):61-67, 2009