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))])
>
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