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

  GenStat  ~  REML meta-analysis

Posted: Fri Mar 15, 2019 10:13 am Reply with quote
Joined: 30 Jul 2009 Posts: 283
Dear Kirsty,

Have you tried VASERIES/VAMETA – this will try the various models and keep
the best fitting model for each experiment and then combines them.
If your structure fits within this framework, it’s definitely an easy way to go.

Regards, David.
Dr David Baird Statistical Consultant and Genstat Developer
VSN (NZ) Limited (David@VSN.CO.NZ (David@VSN.CO.NZ))
8 Mariposa Crescent, Aidanfield, Christchurch 8025, New Zealand
Ph +64 3 3350588 Cell +64 21 1160803

From: GENSTAT-Request <GENSTAT@jiscmail.ac.uk> On Behalf Of Kirsty Hassall
Sent: 15 March 2019 9:59 PM
Subject: REML meta-analysis

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] \
!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,

Dr Kirsty Hassall
Scientific Statistician

Department of Computational and Analytical Sciences
Rothamsted Research
Harpenden, Herts

Tel: +44(0) 1582 763133 Ext: 2080
Email: kirsty.hassall@rothamsted.ac.uk (kirsty.hassall@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.

To unsubscribe from the GENSTAT list, click the following link:

Post generated using Mail2Forum (http://www.mail2forum.com)
View user's profile Send private message

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