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
[Solved] Methodology for prediction equations of phenotypic
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.
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.
-
- Posts: 2
- Joined: Wed May 18, 2011 8:41 pm
[Solved] Methodology for prediction equations of phenotypic
Last edited by Nicolas Guillemin on Mon Jun 13, 2011 9:14 pm, edited 1 time in total.
-
- Posts: 2
- Joined: Mon May 23, 2011 1:30 pm
Re: Methodology for prediction equations of phenotypic data
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
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
-
- GenABEL senior expert
- Posts: 151
- Joined: Wed Feb 09, 2011 3:24 pm
Re: Methodology for prediction equations of phenotypic data
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:
Does this answer to your question?
Hope I was of help.
Nicola
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
-
- Posts: 2
- Joined: Wed May 18, 2011 8:41 pm
Re: Methodology for prediction equations of phenotypic data
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
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 1 guest