sample()and "times" in qtscore()

Questions about GenABEL (aka *ABEL) suite of packages
Forum rules
Please remember not to post any sensitive data on this public forum.
The first few posts of newly registered users will be moderated in order to filter out any spammers.

When get a solution to the problem you posted, please change the topic name (e.g. from "how to ..." to "[SOLVED] how to ..."). This will make it easier for the community to follow the posts yet to be attended.
Posts: 2
Joined: Thu Oct 09, 2014 2:43 am

sample()and "times" in qtscore()

Postby snowy » Tue Nov 03, 2015 2:02 am

I used qtscore() in GenABEL to do GWAS. But my smaple size is small. Someone suggested me to do bootstraps without replacement
and gave me a formula (A)


Code: Select all

qtscore(sample(phdata(b.dat)$disease), replace=F), b.dat, trait ="binomial")

but I think this formula is wrong, because it only sampled phentypes but without relevant genotype information.

I think the right formular is (B)


Code: Select all

qtscore(disease, b.dat,trait ="binomial", idsubset=sample(gtdata(b.dat)@idnames,replace=F),times=1000)

am I right? or anyother good suggestions?

I also have a question:

what's the meaning of "times" in qtscore(), if times=1000, does it mean that in my formula the sample process run 1000 times?

if times=1000, what's the meaning of P1df in results? How can I calculate FDR(false discovery rate)?

Thank you!
Last edited by lckarssen on Tue Nov 03, 2015 10:13 am, edited 2 times in total.
Reason: Some formatting of code

Return to “GenABEL”

Who is online

Users browsing this forum: No registered users and 1 guest