Question: R: How to test for trend of risk scores
gravatar for m93
23 months ago by
m93170 wrote:

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.


risk scores trend test R • 1.1k views
ADD COMMENTlink modified 23 months ago by Biostar ♦♦ 20 • written 23 months ago by m93170

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


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

ADD REPLYlink modified 23 months ago • written 23 months ago by Nicolas Rosewick8.0k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 855 users visited in the last hour