Re: [R] Dominance in qtl model

2007-04-23 Thread Charles C. Berry
On Mon, 23 Apr 2007, joseph powell (RI) wrote: > 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))]) >

[R] Dominance in qtl model

2007-04-23 Thread joseph powell \(RI\)
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="n