Perhaps I should have added the following: To see that it "works," run the
following:
famfit<-lmer(peg.no~1 + (1|family), na.action=na.exclude, vcdf)
resfam<-residuals(famfit)
for( i in 1:length(colms))
{
print(coef(lm(resfam~colms[,i])))
}
Regards, Mark.
A Singh wrote:
>
>
> Dear All,
>
Hi Aditi,
Parts of _your_ code for the solution offered by Jerome Goudet are wrong;
see my comments.
> famfit<-lmer(peg.no~1 + (1|family), na.action=na.omit, vcdf) ## use:
> na.action=na.exclude
> resfam<-residuals(famfit)
> for( i in 1:length(colms))
+ {
+ print ("Marker", i)
+ regfam<-ab
Dear All,
I am quite new to R and am having a problem trying to run a linear model
with random effects/ a regression- with particular regard to my variable
lengths being different and the models refusing to compute any further.
The codes I have been using are as follows:
vc<-read.table("P:\
3 matches
Mail list logo
