Hello,
A couple of questions.
1. I upgraded to asreml 3.0 and I'm running the model:
met.1.asr <- asreml(yld ~ Site, random =~ diag(Site):geno, rcov=~at(Site):ar1(Column):ar1(Row), data=mmData). The sample data is nicely balanced.
1 2 3 4 5 6
SITE1 72 72 72 72 72 72
SITE2 72 72 72 72 72 72
SITE3 72 72 72 72 72 72
However, with this 3.0 version I cannot plot the variogram as I used to with the previous version. I now get an error:
> plot(variogram(met.1.asr))
Error in `$<-.data.frame`(`*tmp*`, "gamma", value = c(0, 0.415962583409219, :
replacement has 1296 rows, data has 0
2. To extract solutions from a simple genetic model fitting pedigree
as "random =~ ped(Pedigree)".
coef(model1.asr)$random
effect
ped(Pedigree)_GARDEL_URUGUAYO 2.23904
I first tried this code:
coef(model1.asr,pattern='ped(Pedigree)'), but it does not work. It returns just the name of the effect, but no solutions.
However, I noticed that the following code does work (note that 'ped(' is missing)
coef(model1.asr,pattern='Pedigree)')
Are there any rules that I'm missing when using "pattern"? It looks like it does not like "("
Thanks for you for your assistance.
Fernando
-------------------- m2f --------------------
Sent using Mail2Forum (http://www.mail2forum.com).
Read this topic online here:
http://www.vsni.co.uk/forum/viewtopic.php?p=1265#1265
-------------------- m2f --------------------
Received on Thu Jan 21 2010 - 21:06:28 EST
This webpage is part of the ASReml-l discussion list archives 2004-2010. More information on ASReml can be found at the VSN website. This discussion list is now deprecated - please use the VSN forum for discussion on ASReml. (These online archives were generated using the hypermail package.)