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
Hi LCV developer,
I'm using LCV for causality estimation on summary statistics:
The output rho.est (0.27, se0.07) had an opposite sign compared with rg (-0.23,0.03) from LDSC:
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