[Solved] Methodology for prediction equations of phenotypic

[Nicolas Guillemin]

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

[sumankundu]

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]

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]

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