This is a test version of Biostars. For the public version, visit https://www.biostars.org.
R: How to test for trend of risk scores

I have some case/control data like so (many more rows)

ID    score     phenotype
1     0.6482772         0
2     0.9519315         1
3     2.7689036         0
4     1.9974466         1
5     1.6287841         1
6     1.9993825         0

I performed a logistic regression on pheno~score to determine if an increase in scores leads to increased odds ratio (a measure of disease risk)

I now want to do a trend test in R to test whether there is trend. This trend test would just give me one odds ratio and a p-value. Does anyone know which R package or function I could use? I've looked around and I'm struggling to find a simple package that can help me do this.

Thanks

r trend test risk scores
• 2,828 views • link
updated by Biostar Biostar • written by m98 m98

As phenotype is either 0 or 1 could you just do a fisher test as :

fisher.test(score,phenotype)

Could also be nice to plot the data as boxplots (x is the data.frame containing your data):

ggplot(x,aes(x=score,y=score))+geom_boxplot()
• 0 views • link
written by Nicolas Rosewick Nicolas Rosewick

0 answers

No answers yet.

Log in to answer this question.

More posts like this