[vsni.co.uk] Contact us
Author Message
Post new topic

  GenStat  ~  REML meta-analysis

Kirsty Hassall
Posted: Fri Mar 15, 2019 8:59 am Reply with quote
Guest
Dear Genstatters,

I am attempting to fit a meta-analysis across 8 different field experiments, each of which has a separate residual structure. The commands to set this up are:

vrmeta [terms=Environment*Genotype; experiments=Environment; random=random] \
'2007','2008R','2008W','2009','2010','2011','2012','2014';\
localterms=!f(block/subblock),!f(block/subblock),!f(block/subblock),\
!f(block/subblock),!f(block/subblock),!f(block/subblock), \
!f(block/(row12F + col12F)), !f(block)
vcomp [Genotype; Experiments=Environment] #random
vresid [experiment='2008R'; term=row08rF.col08rF] model=ar,ar; factor=row08rF,col08rF
vresid [experiment='2010'; term=row10F.col10F] model=ar,ar; factor=row10F,col10F
vresid [experiment='2011'; term=row11F.col11F] model=AR; order=1; factor=col11F
vresid [experiment='2012'; term=block.row12F.col12F] model=AR; order=1; factor=col12F
vresid [experiment='2014'; term=row14F.col14F] model=AR; order=1; factor=col14F
reml [mvincl=y,e; work=500] GY

Running the above, produces output for the random effects and residual terms which all seem to match the structure I would expect (i.e. no additional residual terms are added etc…). However, before reaching the output for the fixed effects I have the error “Genstat server applictaion has stopped working: A problem caused the program to stop working”.

Any advice on how to cope with this error would be greatly appreciated. Note, the above works if I remove the spatial structure over the residuals.

Thanks in advance,
Kirsty

-------------------------------------------------------------------------------------------------------------------------
Dr Kirsty Hassall
Scientific Statistician

Department of Computational and Analytical Sciences
Rothamsted Research
Harpenden, Herts
AL5 2JQ, UK

Tel: +44(0) 1582 763133 Ext: 2080
Email: kirsty.hassall@rothamsted.ac.uk (kirsty.hassall@rothamsted.ac.uk)
http://www.rothamsted.ac.uk


Rothamsted Research is a company limited by guarantee, registered in England at Harpenden, Hertfordshire, AL5 2JQ under the registration number 2393175 and a not for profit charity number 802038.

Post generated using Mail2Forum (http://www.mail2forum.com)

Display posts from previous:  

All times are GMT
Page 1 of 1
Post new topic

Jump to:  

You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum
You can attach files in this forum
You can download files in this forum