

# Two sample MR report

---

## family.Rikenellaceae.id.967 against BK polyomavirus seropositivity

Date: **22 January, 2026**

---

### Results from two sample MR:


|method                    | nsnp|         b|        se|      pval|
|:-------------------------|----:|---------:|---------:|---------:|
|MR Egger                  |   20| 0.1541870| 0.9789421| 0.8766013|
|Weighted median           |   20| 0.7101236| 0.4245308| 0.0943813|
|Inverse variance weighted |   20| 0.6733523| 0.3151565| 0.0326334|
|Simple mode               |   20| 0.2860323| 0.7082861| 0.6908415|
|Weighted mode             |   20| 0.1982180| 0.7051815| 0.7816811|

---

### Heterogeneity tests


|method                    |        Q| Q_df|    Q_pval|
|:-------------------------|--------:|----:|---------:|
|MR Egger                  | 13.55152|   18| 0.7578241|
|Inverse variance weighted | 13.86529|   19| 0.7914865|

--- 

### Test for directional horizontal pleiotropy


| egger_intercept|       se|     pval|
|---------------:|--------:|--------:|
|       0.0385902| 0.068892| 0.582279|

--- 

### Test that the exposure is upstream of the outcome


| snp_r2.exposure| snp_r2.outcome|correct_causal_direction | steiger_pval|
|---------------:|--------------:|:------------------------|------------:|
|        3.12e-05|          5e-07|TRUE                     |    0.7160603|

Note - R^2 values are approximate

---

### Forest plot of single SNP MR

![plot of chunk chunk5](figure/familyRikenellaceaeid967_against_BK_polyomavirus_seropositivitychunk5-1.png)

---

### Comparison of results using different MR methods

![plot of chunk chunk6](figure/familyRikenellaceaeid967_against_BK_polyomavirus_seropositivitychunk6-1.png)

---

### Funnel plot

![plot of chunk chunk7](figure/familyRikenellaceaeid967_against_BK_polyomavirus_seropositivitychunk7-1.png)

---

### Leave-one-out sensitivity analysis

![plot of chunk chunk8](figure/familyRikenellaceaeid967_against_BK_polyomavirus_seropositivitychunk8-1.png)
