[BioC] SNP Analysis
Scholz Claus-Jürgen
scholz at klin-biochem.uni-wuerzburg.de
Thu Oct 22 22:09:13 CEST 2009
Dear Peter,
PLINK is simply a convenient implementation of analysis tasks usually
performed in association studies. If that's not enough, further
processing of PLINK output is easily done with R - it's designed for
that purpose. PLINK is actually standard software, it's freely available
and also well documented on the website
(http://pngu.mgh.harvard.edu/~purcell/plink/) therefore it might be
worth having at look at it. If you prefer R/Bioconductor, there are the
Packages GeneticsBase, SimHap (my personal favourite), GenABEL,
SNPassoc...
Bests,
Claus-Jürgen
>>> Peter Ganske <peterganske at mac.com> >>>
Dear Claus- Jürgen,
thanks for the reply. In which way you would analyze the genotype
frequency wit PLINK?
And why you use this program instead of any bioconductor- package?
All the best and thanks in advance
Peter
Am 22.10.2009 um 13:13 schrieb Claus-Jürgen Scholz:
>
> Dear Peter,
>
> indeed, Birdseed is a genotyping algorithm and I'd use it for genotype
> calling of SNP6.0 arrays (best suited for this platform). If you have
> the calls, export them into a table (export options and formats should
> be described in the Genotyping Console manual) and analyze the
> genotype
> frequency differences between responders and non-responders (valuable
> free software is e.g. PLINK). However, n=100 is a pretty small sample
> size for a GWAS...
>
> Bests,
> Claus-Jürgen
>
>
> Peter Ganske schrieb:
>> Dear Vincent,
>> thanks for the fast replay. Well, i thought, that the Genotyping
>> console used the Birdseed Algorithm and this algorithm is an
>> Genotyping Algorithm.
>>
>> Its hard to find paper or groups, who worked with this array and for
>> me ( i work as a student for an institue) is hard to find the right
>> workflow without help (nobody worked here with SNP arrays in the
>> past)
>>
>> So, i have 100 Arrays (100 CHP and 100 CEL files) of 100 patients. I
>> want to have a look at the SNPs of the patients. 50 are non-responder
>> and 50 are responder. There should be a difference between the two
>> groups. Since yet, i looked for any papers for getting an "general"
>> workflow for sorting out most of the SNPs of the patients.
>>
>> So you think i have to try this package and create the genotyping
>> calls?
>> Whats about this workflow? So are my following thought right:
>>
>> - The package check every SNP for every Chips and put the result in a
>> table
>> - i can combine the result of the SNPs with a selection of gene i
>> want....
>>
>> My boss talked about a top-list of 50 genes... Maybe this can help me
>> out for the usage of CRLMM.. dont know
>>
>> Thanks a lot and sorry for the questions. First time for me to work
>> with SNP Arrays and the first time to work with Bioconductor/R
>> All the best from Germany
>> Peter
>> Am 21.10.2009 um 16:11 schrieb Vincent Carey:
>>
>>
>>> Briefly, you can perform genotype calling with a confidence measure
>>> using crlmm package, working from the CEL files. The crlmm package
>>> includes a vignette called crlmmDownstream.pdf that illustrates one
>>> approach to GWAS analysis based on 6.0, using snpMatrix package. To
>>> use crlmm you will also need a metadata package called
>>> genomewidesnp6crlmm.
>>>
>>> There are certainly other approaches possible. Our workflow
>>> documentation for this use case probably needs some enhancement.
>>>
>>> On Wed, Oct 21, 2009 at 9:42 AM, Peter Ganske
<Peter.Ganske at hki-jena.de
>>>
>>>> wrote:
>>>>
>>>> Hello,
>>>> first time for me to work with SNP arrays. I got CEL- and CHP-files
>>>> for my Analysis. The CEL are from Affymetrix Human-Wide Genome SNP-
>>>> Array 6.0 and the CHP- files are dealed with the Birdseed-
>>>> Algorithm (part of the Genotyp Console from Affymetrix as well).
>>>> Is there anybody here, who worked with this arrays in the past? I
>>>> am looking for an (general) workflow for my study. I want to
>>>> analyse patients with Rheumatoid Arthritis w
ith regard to SNPs and
>>>> the question "why there are respo>>>> I am looking for an workflow for the arrays. Is it better to work
>>>> with the CHP files or with the CEL- files?
>>>> Would me great, if anybody can help me out.
>>>> Thanks in advance
>>>> Peter
>>>>
>>>>
>>>> The information contained in this email and any attachments is
>>>> confidential and may be subject to copyright or other intellectual
>>>> property protection. If you are not the intended recipient, you are
>>>> not authorized to use or disclose this information, and we request
>>>> that you notify us by reply mail or telephone and delete the
>>>> original message from your mail system.
>>>> [[alternative HTML version deleted]]
>>>>
>>>> _______________________________________________
>>>> Bioconductor mailing list
>>>> Bioconductor at stat.math.ethz.ch
>>>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>>>> Search the archives:
http://news.gmane.org/gmane.science.biology.informatics.conductor
>>>>
>>>>
>>> The information contained in this email and any attachments is
>>> confidential and may be subject to copyright or other intellectual
>>> property protection. If you are not the intended recipient, you are
>>> not authorized to use or disclose this information, and we request
>>> that you notify us by reply mail or telephone and delete the
>>> original message from your mail system.
>>>
>>> The information contained in this email and any attachments is
>>> confidential and may be subject to copyright or other intellectual
>>> property protection. If you are not the intended recipient, you are
>>> not authorized to use or disclose this information, and we request
>>> that you notify us by reply mail or telephone and delete the
>>> original message from your mail system.
>>>
>>
>>
>> [[alternative HTML version deleted]]
>>
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor at stat.math.ethz.ch
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>> Search the archives:
http://news.gmane.org/gmane.science.biology.informatics.conductor
>>
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
http://news.gmane.org/gmane.science.biology.informatics.conductor
>
> The information contained in this email and any attachments is
> confidential and may be subject to copyright or other intellectual
> property protection. If you are not the intended recipient, you are
> not authorized to use or disclose this information, and we request
> that you notify us by reply mail or telephone and delete the
> original message from your mail system.
More information about the Bioconductor
mailing list