Question: How to run bioHMM on Affymetrix data?
0
gravatar for Stacey Burrows
10.7 years ago by
Stacey Burrows30 wrote:
I am interested in fitting bioHMM on some non-array CGH data using the function runBioHMM. Say I have have a vector, y, of processed intensity measurements at M SNPs along a chromosome and a vector of corresponding genomic locations, say pos. How can I use this to run bio HMM? Thanks for any help, Stacey __________________________________________________________________ [[elided Yahoo spam]] [[alternative HTML version deleted]]
cgh • 418 views
ADD COMMENTlink modified 10.7 years ago by Oosting, J. PATH550 • written 10.7 years ago by Stacey Burrows30
Answer: How to run bioHMM on Affymetrix data?
0
gravatar for Oosting, J. PATH
10.7 years ago by
Oosting, J. PATH550 wrote:
You will have to create a SegList object in order to do this, for the runBioHMM() function a SegList needs to have the following components M.observed: a matrix with the copynumber values samples in columns, probes in rows genes: A data.frame with at list a Chr and a Position column Make sure that the data is already ordered by chromosome and position # example code starting with vectors of cn.data, chromosome and # base.position Idx<-order(chromosome,base.position) segList<-list() observed<- matrix(cn.data[Idx],ncol=1) colnames(observed)<-"sample1" segList$M.observed <- observed segList$genes <- data.frame(Chr=chromosome[Idx],Position=base.position[Idx]) segList<-new("SegList", segList) result<-runBioHMM(segList) Jan > I am interested in fitting bioHMM on some non-array CGH data using the > function runBioHMM. Say I have have a vector, y, of processed intensity > measurements at M SNPs along a chromosome and a vector of corresponding > genomic locations, say pos. How can I use this to run bio HMM? >
ADD COMMENTlink written 10.7 years ago by Oosting, J. PATH550
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 293 users visited in the last hour