search for: testvcomp10

...th the na.exclude command, because there is a lot of missing data in patches which I have replaced with NA's, but I don't know how to fix this error at all. Its a bit urgent, and any help is hugely appreciated. The data files are at: <http://www.4shared.com/file/131980362/460bdafe/Testvcomp10.html> (excel) http://www.4shared.com/file/131980512/dc7308b/Testvcomp10.html (txt) Cheers, Aditi ---------------------- A Singh Aditi.Singh at bristol.ac.uk School of Biological Sciences University of Bristol

...eg.no.prm<-sample(peg.no, length(peg.no)) try(fit5<-lmer(data=modeldf, peg.no.prm~1 + (1|family/f))) print(summary(fit5)) capture.output(fit5, file="testperm5.txt", append=T) }}} } The data files are at: >> >> >> <http://www.4shared.com/file/131980362/460bdafe/Testvcomp10.ht >> ml> (excel) >> http://www.4shared.com/file/131980512/dc7308b/Testvcomp10.html >> (txt) ------------------------------------------------------------------ ---------------------- A Singh Aditi.Singh at bristol.ac.uk School of Biological Sciences University of Bristol

...; (family wise, which is why this is the blocking factor), explained by each PxLyy, and by 'family' as well. I have managed to run the lmer model for the first marker column (P1L55), and have been able to block it according to 'family' as follows: vc<-read.table("P:\\R\\Testvcomp10.txt",header=T) attach(vc) family<-factor(family) colms<-(vc)[4:13] ## these are the marker columns lapply(colms,factor) library(lme4) for (f in levels(family)) { print("Family") print(f) { t<- vc [ family == f, ] try(lmer<-fit(peg.no~1 + (1|P1L55) + (1|family), data=t,...

..., 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:\\R\\Testvcomp10.txt",header=T) >> attach(vc) > > family<-factor(family) > colms<-(vc)[,4:13] ## this to assign the 10 columns containing marker > data to a new variable, as column names are themselves not in any > recognisable sequence > > vcdf<-data.frame(family,peg.no,...

Dear List, I'm having a curious problem with lapply(). I've used it before to convert a subset of columns in my dataframe, to factors, and its worked. But now, on re-running the identical code as before it just doesn't convert the columns into factors at all. As far as I can see I've done nothing different, and its strange that it shouldn't do the action. Has anybody

....") print(se[i]) { peg.no.prm<-sample(peg.no, length(peg.no)) try(fit5<-lmer(data=modeldf, peg.no.prm~1 + (1|family/f))) print(summary(fit5)) capture.output(fit5, file="testperm5.txt", append=T) }}} } The data files are at: <http://www.4shared.com/file/131980362/460bdafe/Testvcomp10.ht ml> (excel) http://www.4shared.com/file/131980512/dc7308b/Testvcomp10.html (txt) ---------------------- A Singh Aditi.Singh at bristol.ac.uk School of Biological Sciences University of Bristol

...w to fix this >> >> error at all. >> >> >> >> Its a bit urgent, and any help is hugely appreciated. >> >> >> >> The data files are at: >> >> >> >> >> >> <http://www.4shared.com/file/131980362/460bdafe/Testvcomp10.ht >> >> ml> (excel) >> >> http://www.4shared.com/file/131980512/dc7308b/Testvcomp10.html >> >> (txt) >> >> >> >> >> >> Cheers, >> >> Aditi >> >> >> >> ---------------------- >> >...

...marker), but I have been advised to use that as an alternative corroborative approach for later, and stick to the original looped lmer() as my principal approach. For that, the code that I was using, with some fresh modifications, is as follows: > vc<-read.table("E:\\R\\Text Files\\Testvcomp10.txt", header=T, sep="", dec=".", na.strings=NA, strip.white=T) > attach(vc) > colms<-(vc)[4:13] > family<-factor(family) > for (c in 1:levels(family)) + { + print("Family") + print(c) + { + for (i in 1:length(colms)) + { + print(anova(fit&lt...