SNPassoc association and weights

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

SNPassoc association and weights

Thierry van de Wetering
Dear Reader,

Please can you tell me how to use weights with the SNPassoc package
function association.

association(NEB~recessive(rs178012)*as.factor(GENDER), weights =
as.factor(WAGA), data=dat2SNP)

SNP: recessive(rs178012 adjusted by: Interaction --------------------- 1
dif lower upper 2 dif lower upper G/G-A/G 1485 1.451 0.03814 0.0000 NA
NA 1702 1.854 0.03879 0.4025 0.2946 0.5105 A/A 11 1.091 0.45636 -0.3603
-1.28 0.5599 20 1.900 0.62786 0.4488 -0.2356 1.1333 p interaction: 0.48707


In this example I want to know if there is an interaction between GENDER
and the SNP in assocation with NEB, weighted for WAGA. Removal of the
weights argument gives exactly the same p.value.

association(NEB~recessive(rs178012)*as.factor(GENDER), data=dat2SNP)

SNP: recessive(rs178012 adjusted by: Interaction --------------------- 1
dif lower upper 2 dif lower upper G/G-A/G 1485 1.451 0.03814 0.0000 NA
NA 1702 1.854 0.03879 0.4025 0.2946 0.5105 A/A 11 1.091 0.45636 -0.3603
-1.28 0.5599 20 1.900 0.62786 0.4488 -0.2356 1.1333 p interaction: 0.48707



Changing the model from recessive into for example dominant gives me
other p interactions and p.values, but still the weighting is not giving
any differences.

Many thanks in advance.
Your sincerely

--
Thierry


        [[alternative HTML version deleted]]

______________________________________________
[hidden email] mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.