[Solved] Methodology for prediction equations of phenotypic

Statistical Genomics for dummies and advanced. Discussions, links, usefull information.
Forum rules
Welcome! Please feel free to raise any issue. There is no issues big or small. Let's work on them together.

Please note that the first few posts of newly registered users will be moderated in order to filter out any spammers.
Nicolas Guillemin
Posts: 2
Joined: Wed May 18, 2011 8:41 pm

[Solved] Methodology for prediction equations of phenotypic

Postby Nicolas Guillemin » Fri May 20, 2011 7:35 pm

Hi
I'm working on a methodology for developing prediction equations of phenotypic data, based on SNPs, calculated with the lm function.
I saw that it's possible in the GenABEL tutorial.

I'm wondering if someone has experience in this domain, or if there is some interresting bibliography on it (I don't find publications with GenABEL on prediction model for phenotype based on SNPs) ? My interest is to demonstrate that this type of methodology, using lm and GenABEL, could be reliable for the project.

Thanks
Last edited by Nicolas Guillemin on Mon Jun 13, 2011 9:14 pm, edited 1 time in total.

sumankundu
Posts: 2
Joined: Mon May 23, 2011 1:30 pm

Re: Methodology for prediction equations of phenotypic data

Postby sumankundu » Mon May 23, 2011 1:47 pm

Dear Nicolas,
I am currently associated with comparing different simulation methods which have been used in (geneic) risk prediction studies and working with Cecile Janssens group in Erasmus MC, Rotterdam.
Following links may be usefull in your project.

1. PredictABEL: an R package for the assessment of geentic risk prediction model, developed by us.
weblink: http://www.genabel.org/packages/PredictABEL

2. Strengthening the reporting of genetic risk prediction studies: the GRIPS statement.
weblink: http://www.ncbi.nlm.nih.gov/pubmed/21431409

Best wishes,
Suman

Nicola Pirastu
GenABEL senior expert
GenABEL senior expert
Posts: 151
Joined: Wed Feb 09, 2011 3:24 pm

Re: Methodology for prediction equations of phenotypic data

Postby Nicola Pirastu » Thu May 26, 2011 1:55 am

Hi,

I'm not sure I completely understood what you're asking could you be a little more specific?
I don't think you can find any reference for using GenABEL to do linear regression with lm since lm is a function in the base package of R.

Anyways if you're trying to do something like:

trait~sex+age+SNP1+SNP2+SNP3

so let's say your trait is qt1 from the srdta example data and you want to use rs18 rs29 and rs73 as predictors with lm, you could do something like:

Code: Select all

library(GenABEL)
data(srdta)

df=as.data.frame(as.numeric(srdta[,c("rs18","rs29","rs73")]))
df$qt1=phdata(srdta)$qt1
res=lm(qt1~.,data=df)
summary(res)


Does this answer to your question?

Hope I was of help.

Nicola

Nicolas Guillemin
Posts: 2
Joined: Wed May 18, 2011 8:41 pm

Re: Methodology for prediction equations of phenotypic data

Postby Nicolas Guillemin » Mon Jun 13, 2011 9:13 pm

Thanks for your answers.

According scpecialized people, the GenABEL package is accurate for my project (modeling phenotypic data from SNPs).
So, I will use this package during my work


Return to “Journal Club on Statistical Genomics”

Who is online

Users browsing this forum: No registered users and 2 guests