[BioC] genotyping output files
Vincent Carey
stvjc at channing.harvard.edu
Sun Jan 15 01:32:26 CET 2012
Dear Rafiq S.
You do not have to reissue your query. You can verify that it was
broadcast by checking the mailing list archive.
There are to my knowledge no immediate translations from crlmm output
to plink input. The "downstream" vignette
http://www.bioconductor.org/packages/2.9/bioc/vignettes/crlmm/inst/doc/crlmmDownstream.pdf
shows some relevant conversions, but the document is obsolete and
needs to be rebuilt. If you substitute "snpStats" for
"snpMatrix" and "SnpMatrix" for "snp.matrix", everything should work
as noted. You can use coercion of the SnpMatrix instance to
"character" to get the calls in "A/B" format.
On Sat, Jan 14, 2012 at 6:42 PM, Rafiq S. <S.Rafiq at soton.ac.uk> wrote:
>
> Dear All,
>
> I have recently attempted a genome wide genotyping experiment on around 500 samples using the Illumina 660quad chip.
> Crlmm genotyping was completed and a snpset object crlmmResult was created.
>
> I have used the calls(crlmmResult) and the configuation command to quickly browse through the genotypes.
>
> I was expecting that at the end of the process crlmm would produce a text file with all the genotypes which could be transported to other softwares like plink for association analysis,
> This hasnt happened and a text file was not generated, however a crlmmResult.Rda file was extracted into SNPchip packages in my rlibrary.
>
> Could you please direct me towards what other prrocesses which I might need to complete before I can obtain an output text file all the with genome wide calls?
>
> Many Thanks,
>
> -S
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor
mailing list