Bookmark and Share

Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down on April 23, and its replacement, statalist.org is already up and running.


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: st: Re: Anovalator output


From   Martyn Sherriff <statams48@gmail.com>
To   statalist@hsphsun2.harvard.edu
Subject   Re: st: Re: Anovalator output
Date   Thu, 8 Dec 2011 14:36:28 +0000

Thanks for the advice. I will check for singletons and see what trace yields.

Martyn

On 8 December 2011 14:17, Joseph Coveney <jcoveney@bigplanet.com> wrote:
> Martyn Sherriff wrote:
>
> I am running xtmixed on a repeated measures design that is unbalanced,
> there are 4 levels of treatment and 7 levels of time. I am then using
> anovalator (from http://www.ats.ucla.edu/stat/stata/ado/analysis) to
> investigate the main and simple effects.
> xtmixed f trt##time || id: , reml
>
> . anovalator time, main
> anovalator main-effect for time
> chi2(1) = 134.57337   p-value = 4.093e-31
>
> . anovalator trt, main
> too few variables specified
> r(102);
>
> Can anybody suggest why I get this error message for 'trt'.
>
> --------------------------------------------------------------------------------
>
> How imbalanced are things?  Singletons?  -anovalator- seems to work with 15%
> missing values (below).  Perhaps you could -set trace on-, see just what command
> inside balks, and look further from there.
>
> Joseph Coveney
>
> . clear *
>
> . set more off
>
> . set seed `=date("2011-12-08", "YMD")'
>
> . set obs 200
> obs was 0, now 200
>
> . generate int pid = _n
>
> . generate double u = rnormal()
>
> .
> . generate byte group = mod(_n, 4)
>
> .
> . forvalues time = 1/7 {
>  2.         generate double response`time' = u + ///
>>                 `time' / 7 * 0.1 + rnormal()
>  3. }
>
> . quietly reshape long response, i(pid) j(time)
>
> .
> . drop if runiform() < 0.15
> (198 observations deleted)
>
> .
> . quietly xtmixed response i.group##i.time || pid: , ///
>>         reml nolrtest nostderr nolog
>
> . anovalator time, main
>
> anovalator main-effect for time
> chi2(6) = 7.3372743   p-value = .29077902
>
> . anovalator group, main
>
> anovalator main-effect for group
> chi2(3) = .20230331   p-value = .97721661
>
> . exit
>
> end of do-file
>
>
> *
> *   For searches and help try:
> *   http://www.stata.com/help.cgi?search
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/


© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index