Displaying 1 result from an estimated 1 matches for "myfun4".
Did you mean:
myfun
2007 Apr 23
1
Dominance in qtl model
Hi,
I'm using R for a QTL analysis of SNP data. I was wondering if anyone
had any advice on fitting a dominance effect into the following
function;
> myfun4
function (x) {
x <- scan(con, nmax=169)
y <- unique(x[which(!is.na(x))])
if(length(y)>1) {
summary(lme(Ad ~ x, random= ~1|sire, na.action="na.omit"))
}
else {print("no.infomation")}
}
Con is the connection to a file of the genotypes for each SNP. It is set
up as a con...