Displaying 1 result from an estimated 1 matches for "snp900".
Did you mean:
snp500
2009 Sep 22
2
glm analysis repeated for 900 variables
...s,
Could you help my with the following problem?
I want to repeat a glm analysis with 2 independent variables for all 900
variables (snps) in my data set. So, I want to check whether snp1 has a
different effect on my outcome variable in patients and
controls(phenotype). And repeat that for snp2 to snp900.
Is there an easy way to get a summary of the data, e.g. a list of P
values of all 900 variables?
I tried something with a loop:
for (i in 1:length(data)) { print (summary (glm
(outcome~data[[i]]*phenotype, data=data))) } # This works, but gives
900 written summaries
for (i in 1:length(data))...