Skip to content

rg by LCV has different sign with rg by LDSC #8

Description

@WeiCSong

Hi LCV developer,
I'm using LCV for causality estimation on summary statistics:

int=fread("ID94.sumstats.gz",data.table=FALSE) #data from 30837455
d1=fread("ID644.sumstats.gz",data.table=FALSE) # data of neuroticism from pmid30643256
d1=d1[which(!is.na(d1$Z)),]
rownames(d1)=d1[,1]
rownames(int)=int[,1]
l=rownames(ldsc) #ldsc is a data.frame made by ld scores under eur_w_ld_chr folder, and i set its
rownames as rsid
L=intersect(l,rownames(d1))
L=intersect(L,rownames(int))
exp=d1[L,"Z"]
out=int[L,"Z"]
ld=ldsc[L,"L2"]
RunLCV(ld,exp,out)

The output rho.est (0.27, se0.07) had an opposite sign compared with rg (-0.23,0.03) from LDSC:

C:\Users\goubegou\ldsc\ldsc.py --rg ID94.sumstats.gz,ID644.sumstats.gz --ref-ld-chr
F:\selection\eur_w_ld_chr\ --w-ld-chr F:\selection\eur_w_ld_chr\ --out ulcer

I believe I haven't flipped the effect column during these procedures (I made no modification to .sumstats file ever since their generation from munge.sumstats.py), and I'm not sure whether any critical bias has emerged. Should I simply flip gcp.pm and its z-score, assuming this issue just reflected some accidental errors? Thanks for your help.

Best Regards

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions